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Abstract. This paper has two parts. The first part surveys Euler's work on the con- 
stant 7 = 0.57721... bearing his name, together with some of his related work on the 
gamma function, values of the zeta function, and divergent series. The second part 
describes various mathematical developments involving Euler's constant, as well as 
another constant, the Euler-Gompertz constant. These developments include connec- 
tions with arithmetic functions and the Riemann hypothesis, and with sieve methods, 
random permutations, and random matrix products. It also includes recent results on 
Diophantine approximation and transcendence related to Euler's constant. 
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1. Introduction 



Euler discovered the constant 7, defined by 

/ " 1 
7 := lim I > — — logn 




This is a fundamental constant, now called Euler 's constant. We give it to 50 decimal 
places as 



In this paper we shall describe Euler's work on it, and its connection with values of 
the gamma function and Riemann zeta function. We consider as well its close relatives 

= 1.7810724... and = 0.561459.... We also inquire as to its possible meaning, by 
describing many situations in mathematics where it occurs. 

The constant 7 is often known as the Euler- M as cheroni constant, after later work of 
Mascheroni discussed in ^2.6[ As explained there, based on their relative contributions 
it seems appropriate to name it after Euler. 

There are many famous unsolved problems about the nature of this constant. The 
most well known one is: 

Conjecture 1.0.1. Euler's constant is irrational. 

This is a long-standing open problem. A recent and stronger version of it is the 
following. 

Conjecture 1.0.2. Euler's constant is not a Kontsevich-Zagier period. In particular, 
Euler's constant is transcendental. 

A period is defined by Kontsevich and Zagier |181) to be a complex constant whose real 
and imaginary parts separately are given as a finite sum of absolutely convergent integrals 
of rational functions in any number of variables with rational coefficients, integrated over 
a domain cut out by a finite set of polynomial equalities and inequalities with rational 
coefficients, see §3.141 Many constants are known to be periods, in particular all zeta 
values {C(^) '■ n > 2}. The set of all periods forms a ring V which includes the field Q 
of all algebraic numbers. It follows that if Euler's constant is not a period, then it must 
be a transcendenta numberl. Conjecture 11.0.21 also implies that 7 would be Q-linearly 
independent of all odd zeta values C(3), C(5), C(7), ••• and of vr. 

This paper also presents results on another constant defined by Euler, the Euler- 
Gompertz constant 



Here we adopt the notation of Aptekarev [13] for this constant (in which 5 follows 7) 
and refer to it by the name given by Finch [1131 Sec. 6.2] (who denotes it C2). Some 
results about S intertwine directly with Euler's constant, see ^3.51 ^3.61 ^3.151 and ^3.161 
Euler's name is also associated with other constants such as e = 2.71828... given by 



He did not discover this constant, but did standardize the symbol "e" for its use. He 
used the notation e in his 1737 essay on continued fractions ([HI p. 120], cf. [95]), where 



7 = 0.57721 56649 01532 86060 65120 90082 40243 10421 59335 93992... 




(1.1) 




n=0 
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he determined its continued fraction expansion 

e = [2;1,2, 1,1,4, 1,1,6, 1,1,8, •••] 



and deduced that it was irrational. Elsewhere he found the famous relation e'^* = — 1. 
The constant e is also conjectured not to be a Kontsevich-Zagier period. It is well known 
to be a transcendental number, consistent with this conjecture. 

In this paper we review Euler's work on the constant 7 and on mathematical topics 
connected to it, and subsequent developments. It is correspondingly divided into two 
parts, which may be read independently. In addition we have aimed to make individual 
subsections separately readable as short essays. 

The first part of this paper considers the work of Euler. Its emphasis is historical, and 
it retains some of the original notation of Euler. It is a travelogue of part of Euler's work 
related to zeta values and various explicit constants. Euler did extensive work related to 
his constant: he returned to its study many times. The basic fact to take away is that 
Euler did an enormous amount of work directly on his constant and related areas, far 
more than is commonly appreciated. 

The second part of this paper addresses mathematical developments made since Eu- 
ler's time concerning Euler's constant. Since Euler's constant is an unusual constant 
that seems unrelated to other known constants, its appearance in different mathemati- 
cal subjects can be a signal of possible connections between these subjects. We present 
many contexts in which Euler's constant appears. These include connections to the Rie- 
mann hypothesis, to random matrix theory, to probability theory and to other subjects. 
There is a set of strong analogies between factorizations of random integers in [l,n] and 
cycle structures of random permutations in the symmetric group Sn, taking N to be 
proportional to logn. In this analogy there appears another constant that might have 
appealed to Euler: the Golomb-Dickman constant, which can be defined as 



where Li{x) := /^^ see ^3.101 The constant A appears in connection with the 
distribution of the longest cycle of a random permutation, while the constant e~'^ appears 
in parallel in connection with the distribution of the shortest cycle. The discovery of 
this analogy is described at the beginning of ^3.111 It was found by Knuth and Trabb- 
Pardo [I75J through noticing a numerical coincidence to 10 decimal places between two 
computed constants, both of which were the Golomb-Dickman constant. 

In passing we present many striking and elegant formulas, which exert a certain fasci- 
nation in themselves. Some of the formulas presented were obtained by specialization of 
more general results in the literature. Other works presenting striking formulas include 
the popular book of J. Havil |155j devoted extensively to mathematics around Euler's 
constant, and Finch jll3l Sects. 1.5 and 4.5]. 

In ^we make concluding remarks and briefly describe other directions of work related 
to Euler's constant. 



Euler (1707-1783) was born in Basel, and entered the University of Basel at age 13, 
to study theology. He also received tutoring in mathematics on the side from Johann 
Bernoulli (1667-1748), who was Professor of Mathematics at the University of Basel. 
Johann recognized his mathematical gift and convinced his family to let him switch to 
mathematics. Euler was excellent at computation, and had an essentially photographic 
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memory. He had great persistence, and returned over and over again to problems that 
challenged him. 

Here we only scratch the surface of Euler's enormous output, which contains over 800 
works. These include at least 20 books, surveyed in Gautschi |122j . We trace a thread 
comprising results corresponding to one of his interests, the determination of specific 
constants given by infinite sums, definite integrals, infinite products, or by divergent 
series, and finding identities giving interrelations among these constants. Ancillary to 
this topic is the related problem of determining numerical values of these constants; 
computing such values allow one to guess and test possible identities. 

In his work Euler made extensive calculations, deriving approximations of numbers to 
many decimal places. He typically reports results of calculations obtained using various 
approximation schemes, and often wants to compare the different values he obtains by 
the different methods. In Euler's numerical results reported below the symbol ~ is 
used to mean "the result of an approximate calculation" . The symbol ~ does not imply 
agreement to the number of digits provided. An underline of a digit in such an expansion 
indicates the first place where it disagrees with the known decimal expansion. 

We refer to Euler's papers by their index numbers in the Enestrom catalogue (|83j). 
These papers may be found online in the Euler archive (|85j). For discussions of Euler's 
early work up to 1750 see Sandifer |261j . [262]. For his work on the Euler (-Maclaurin) 
summation formula see Hofmann |164j and Knoebel et al. ^i73i Chapter 1]. For other 
discussions of his work on the zeta function see Stackel |280j . Ayoub [21], Weil [309^ 
Chap. 3] and Varadarajan |305l Chap. 3]. 

2.1. Background. Euler spent much effort evaluating the sums of reciprocal powers 

oo ^ 

2 

for integer m > 2. He obtained the formula (^(2) = solving the "Basel problem," 
as well as formulas for all even zeta values Ci"^"^)- He repeatedly tried to find analogous 
formulas for the odd zeta values C,{2n + 1), which would explain the nature of these 
numbers. Here he did not completely succeed, and the nature of the odd zeta values 
remains unresolved to this day. Euler's constant naturally arises in this context, in 
attempting to assign a value to C(l)- 

Euler discovered the constant 7 in a study of the harmonic numbers 

n 

Hn:=y^-, (2.1.2) 

i=i ■' 

which are the partial sums of the harmonic series. The series "C(l)" = X^nii ~ diverges, 
and Euler approached it by (successfully) finding a function f{z) (of a real or complex 
variable z) that interpolates the harmonic number values, i.e. it has f{n) = Hn for 
integer n > 1. The resulting function, denoted Hz below, is related to the gamma 
function via (I2.2.2p . In studying it he also studied, for integer m, the related sums 

n ^ 

Hn,m ■■= E (2.1.3) 

i=i 

which we will term m-harmonic numbers; in this notation H^^i = Hn- 

Euler obtained many inter-related results on his constant 7. These are intertwined 
with his work on the gamma function and on values of the (Riemann) zeta function. In 
§2.2-§2.4 we present results that he obtained for Euler's constant, for the gamma function 
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and its logarithmic derivative, and for values of the zeta function at integer arguments, 
concentrating on explicit formulas. In §2.5 we discuss a paper of Euler on how to sum 
divergent series, applied to the example 0! — l! + 2! — 3! + -- - His analysis produced a 
new constant, the Euler-Gompertz constant 5, given above in (jl.ip . In §2.6 we review 
the history of the name Euler-Mascheroni constant, the adoption of the notation 7 for 
it, and summarize Euler 's approach to research. 

Taken together Euler's many papers on his constant 7, on the related values Cl'^) 
for n > 2, and other values show his great interest in individual constants ("notable 
numbers") and on their arithmetic interrelations. 



2.2. Harmonic series and Euler's constant. One strand of Euler's early work in- 
volved finding continuous interpolations of various discrete sequences. For the sequence 
of factorials this led to the gamma function, discussed below. In the 1729 paper \87\ 
E20] he notes that the harmonic numbers Hn = X]fe=i i ^'^^ given by the integrals 

Hn := / dx. 

Jo 1 - a; 

He proposes a function interpolating Hn by treating n as a continuous variable in this 
integral. The resulting interpolating function, valid for real z > 0, is 

H, := / dx. (2.2.1) 

Jo 1 - a: 



Using this definition, he derives the formula 



Hi =2 - 2 log 2. 

2 

The function Hz is related to the digamma function ^p{z) = ^ \ogT{z) by 

= V(2 + 1) + 7, (2.2.2) 

see ^2.31 and also Theorem 13.1.11 

In a paper written in 1731 ("On harmonic progresssions") Euler |90l E43] summed 
the harmonic series in terms of zeta values, as follows, and computed Euler's constant 
to 5 decimal places. 

Theorem 2.2.1. (Euler 1731) The limit 

7 = lim {Hn — log n) 

exists. It is given by the (conditionally) convergent series 

7 = f;(-l)"^. (2.2.3) 

Euler explicitly observes that this series (12.2.3P converges (conditionally) since it is an 
alternating series with decreasing terms. He reports that the constant 7 ~ 0.577218. 

Euler obtains the formula (I2.2.3P using the expansion log(l + x) = 
which was found by Nicholas Mercator (c. 1620-1687) [213| in 1668. Evaluating this 
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formula at x = 1, ^, ^ Euler observes the following in tabular form: 

1 /1\2 1 /l^ 



log2 = 1 , , , 

^ 2 Viy 3 VI 

3 11 /1\2 1 /1n3 

^°^2 = 2-2-12J +3-y +••• 

, 4 1 1 /1\2 1 /1\3 

^°^3 = 3-2-U) +3-U) 

5 11 /1\2 1 /1n3 

^"^4 = i-rUJ +rU] +■• 



Summing by columns the first n terms, he obtains 
One may rewrite this as 



log(n +l) = Hn- ^Hn,2 + ^^n,3 " " " " (2-2.4) 



Hn - log(n + 1) = ^Hn,2 - ^^n,3 + " " " 

Now one observes that for j >2, 

lim Hnj = Hooj = C(j) 

Taking this limit as n — )• oo term by term in (I2.2.4p formally gives the result (I2.2.3p . 

We note that such representations of infinite sums in a square array, to be summed in 
two directions, were already used repeatedly by Johann Bernoulli's older brother Jacob 
Bernoulli (1654-1705) in his 1689 book on summing infinite series (|29j). Another nice 
example of an identity found by array summation is 

(C(2) - 1) + (C(3) - 1) + (C(4) - 1) + . . . = 1, (2.2.5) 

using 

1111 
C(2)-l = _ + _ + _ + _ + ... 

1111 

C(3)-l = ^ + ^ + i3+^ + --- 

1111 
C(4)-l = _ + _ + _ + _ + ... 

since the column sums yield the telescoping series 

1 13 14 111 

^x2 + — X- + — x-H = \ \ \ = 1. 

22 32 2 42 3 1-2 2-3 3-4 

Compare (j2.2.5p to Euler's expression ()2.2.1ip below for Euler's constant. 

In 1732 Euler |88l E25] stated his first form of his summation formula, without a 

proof. He then gives it in detail in his 1735 paper [911 E47], published 1741. Letting s 

denote the sum of the function f{n) over some set of integer values (say n = 1 to A^), he 

writes 

(■ ^dt dH dH 

s = tdn + at + p- — h + ^~r^ + 6tc.. 

J dn dn"^ dn-^ 

where ct = ^)/3 = f — g,7 = f — f + j5,'5 = ^ — § + ^ ~ T^- These linear equations 
for the coefficients are solvable to give a = |,/3 = ^,7 = 0, = —755 etc. Thus 

t 1 dt 1 d^t 

tdn H 1 ; —T + etc. 

2 12 dn 720 dn^ 
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The terms on the right must be evaluated at both endpoints of this interval being 
summed. In modern terms, 

n=l •^'J 

where the i?2fc are Bernoulli numbers, defined in (12.2.6P below, and the sum must be cut 
off at a finite point with a remainder term (cf. [288, Sec. 1.0]). Euler's formulas omit any 
remainder term, and he applies them by cutting off the right side sum at a finite place, 
and then using the computed value of the truncated right side sum as an approximation 
to the sum on the left, in effect making the assumption that the associated remainder 
term is small. 

A similar summation formula was obtained independently by Colin Maclaurin (1698- 
1746), who presented it in his 1742 two- volume work on Fluxions [203]. A comparison 
of their formulas is made by Mills |215j . who observes they are substantially identical, 
and that neither author keeps track of the remainder term. It is appropriate that this 
formula, now with the remainder term included, is called the Euler- Maclaurin summation 
formula. 

Jacob Bernoulli introduced the numbers now bearing his name (with no subscripts), 
in Ars Conjectandi [30j, published posthumously in 1713. They arose in evaluating (in 
modern notation) the sums of n-th powers 



m-i ^ / / I 1 \ 

1 — 1 \ I — n \ / 



k=l \k=0 

We follow the convention that the Bernoulli numbers Bn are those given by the expansion 

t ^B 



1 ^ ri\ 



n=0 

so that 

B^ = , Bo = — , B-^ = 0, Ba = , B^ = 0, Bp. = — , ... 

2' 6 30 42' 

Euler later named the Bernoulli numbers in his 1768 paper (in [104, E393]) ("On 
the sum of series involving the Bernoulli numbers"), and there listed the first few as 
5' M' 5S' M' definition led to the conventioru that the notation B^ corresponds 
to \B2n\ in (|2.2.6p . This convention is followed in the classic text of Whittaker and 
Watson pT, p. 125]. 

In his 1755 book on differential calculus [1001 E212, Part II, Chap. 5 and 6] Euler 
presented his summation formula in detail, with many worked examples. (See Pengelley 
[239j, |240j for a translation and for a detailed discussion of this work.) In Sect. 122 
Euler computes the Bernoulli numbers \B2n\ through n = 15. In Sect. 129 he remarks 
that: 

The [Bernoulli numbers] form a highly divergent sequence, which grows 
more strongly than any geometric series of growing terms, u 



However in [104] and elsewhere Euler does not use subscripts, writing 21, 23, £ etc. 
^"pariter ac BernouUiani 21, ©,£,2) &c. serium maxime divergentem, quae etiam magis increscat, 
quam uUa series geometrica terminis crescentibus procedens." [Translation: David Pengelley [239], 
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In Sect. 142 he considers the series for the harmonic numbers s = Hx = X]n=i ~' where 
his summation formula reads (in modern notation) 

where the constant C to be determined is Euler's constant. He remarks that this series 
is divergent. He substitutes x = 1, stating that s = 1 in ()2.2.7p . obtaining formaUy 

^ _ 1 I-B2I I-B4I l-Bel l-Bsl 

although the right side diverges. Substituting this expression in (I2.2.7P yields 

f I \B2\ ^ \Bi\ \B6\ \Bs 



(1 -Bo \Ba\ \Br\ \Bq\ 
2 2 4 6 8 

In Sect. 143 he sets x = 10 in (|2.2.7|) so that s = Hiq, and gives the numerical value s = 
2.92896 82539 68253 968. The right side contains C as an unknown, and by evaluating 
the other terms (truncating at a suitable term) he finds 

C ^ 0.57721 56649 01532 5 

a result accurate to 15 places. We note that the divergent series on the right side of (|2.2.7p 
is the asymptotic expansion of the digamma function -0(2; + 1) + 7, where 0(x) = 
see Theorem 13.1.31 and in effect Euler truncates the asymptotic expansion at a suitable 
term to obtain a good numerical approximation. 

Over his lifetime Euler obtained many different numerical approximations to Euler's 
constant. He calculated approximations using truncated divergent series, and his papers 
report the value obtained, sometimes with a comment on accuracy, sometimes not. In 
contrast, some of the convergent series he obtained for Euler's constant converge slowly 
and proved of little value for numerics. He also explored various series acceleration 
methods, which he applied to both convergent and divergent series. He compared and 
cross-checked his work with various different numerical methods. In this way he obtained 
confidence as to how many digits of accuracy he had obtained using the various methods, 
and also obtained information on the reliability of his calculations using divergent series. 

In 1765 in a paper studying the gamma function \W3\ E368] (which will be discussed 
in more detail in §2.3) he derives the formula r'(l) = —7. It is given in Sec. 14 of his 
paper as equation (I2.3.6P in §2.3, which is equivalent to the integral formula 

/•oo 

-7 = / e-'^'logxdx. (2.2.8) 
Jo 

This result follows from differentiation under the integral sign of Euler's integral ()2.3.4p . 

In a 1768 paper mainly devoted to relating Bernoulli numbers to values of C(^)5 Euler 
[1041 E393] obtains more formulas for 7 (denoting it O). In Sect. 24 he reports the 
evaluation of 7 found earlier and remarks: 

= 0, 57721 56649 01532 5. This number seems also the more notewor- 
thy because even though I have spent much effort in investigating it, I 
have not been able to reduce it to a known kind of quantity. □ 



O — 0.5772156649015325 qui numerus eo maiori attentione dignus videtur, quod eum, cum olim in 
hac investigatione multum studii consumsissem, nuUo modo ad cognitum quantitatum genus reducere 
valui." [Translation: Jordan Bell] 
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He derives in Sect. 25 the integral formula 



7=/ (r^^-v'^^- ^^-^-^^ 



Here we may note that the integrand 

I -e-y V v\ey-l J ^ V " ^ V n! ' 

^ " 71=1 

In Sect. 26 by change of variables he obtains 

7= C (T^ + T^^dz (2.2.10) 



/o V 1 - ^ log ^ 
In Sect. 27 he gives the following new formula for Euler's constant 

7 = ^(C(2) - 1) + ^(C(3) - 1) + ^(C(4) - 1) + • • • (2.2.11) 
and in Sect. 28 he obtains a formula in terms of log 2 and the odd zeta values: 

7=4-2'°s2 + E(l-2^)(<P'+l)-l 

fc = l 

He concludes in Sect. 29: 

Therefore the question remains of great moment, of what character the 
number O is and among what species of quantities it can be classified^ 

In a 1776 paper, posthumously published in 1785, Euler [1071 E583] studied the con- 
stant 7 in its own right, in a paper titled: "On a memorable number naturally occurring 
in the summation of the harmonic series 0." In Sect. 2 he speculates that the number 
N = e^ should be a notable number: 

And therefore, if x is taken to be an infinitely large number it will then 
be 

. 11 1^1 

1 + - + - + --- + - = C + logx. 

2 3 X 

One may suspect from this that the number C is the hyperbolic [= nat- 
ural] logarithm of some notable number, which we put = iV, so that 
C = log N and the sum of the infinite series is equal to the logarithm of 
the number N ■ x. Thus it will be worthwhile to inquire into the value 
of this number N, which indeed it suffices to have defined to five or six 
decimal figures, since then one will be able to judge without difficulty 
whether this agrees with any known number or not|j 



Manet ergo quaestio magni momenti, cujusdam indolis sit numerus iste 0[:= 7] et ad quodnam 
genus quantitatum sit referendus." [Translation: Jordan Bell] 

^De numero memorabili in summatione progressionis harmonicae naturalis occurrente 
^"Quod si ergo numerus x accipiatur magnus, turn erit 

11 1 ^ , 

l + - + - + --- + -^C + lx 
2 3 z 

istum numerum C[~ 7] esse logarithmum hyperbolicum cuiuspiam numeri notabilis, quem statuamus 
= A'^, ita ut sit C = log A'^, et summa illius serial infinitae aequeture logarithmo numeri Nx, under operae 
pretium erit in valorem huius numeri A*' inquirerre, quem quidem sufiiciet ad quinque vel fex figuras 
decimales definiuisse, quoniam hinc non difiiculter iduicari poterit, num cum quopiam numero cogito 
conueniat nec ne." [Translation: Jordan Bell] 
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Euler evaluates N = 1.78106 and says he cannot connect it to any number he knows. 
He then gives several new formulae for 7, and tests their usefulness for calculating its 
value, comparing the values obtained numerically with the value 0.57721 56649 01532 5 
that he previously obtained. His first new formula in Sect. 6 is 

1 - 7 = lim - 1) + lim - 1) + lim - 1) + • • • 

Using the first 16 terms of this expansion, he finds 7 ~ 0, 57721 69. He finds several 
other formulas, including, in Sect. 15, the formula 

l_log^-7 = _L^(C(3)-l) + ^(C(5)-l) + ^(C(7)-l) + --- , 

from which he finds 7 ~ 0.57721 56649 01791 3, a result accurate to 12 places. He 
concludes the paper with a list of eight formulas relating Euler's constant, log 2, and 
even and odd zeta- values. One striking formula (number VII), involving the odd zeta 
values is 

^ C(2fc + 1) 

In another 1776 paper, titled "T/ie expansion of the integral formula J dx{j^ + j^) 
with the term extended from x = to x = V\ not published until 1789, Euler [109, E629] 
further investigated the integral formula (j2.2.10p for his constant. He denotes its value 
by n, and says: 

the number n, [...] which is found to be approximately 0, 57721 56649 01532 5, 
whose value I have been able in no way to reduce to already known tran- 
scendental measures; therefore, it will hardly be useless to try to resolve 
the formula in many different waysQ 

He then finds several expansions for related integrals. One of these starts from the 

identity, valid for m,n >0, 



dx = log 



log X \ n + 1 

and from this he deduces, for n > 1, that 



'1 -x)" , 
■dx 



logx 

This formula gives the moments of the measure iog(i_^) dx on [0, 1]. 

2.3. The gamma function. Euler was not the originator of work on the factorial func- 
tion. This function was initially studied by John Wallis (1616-1703), followed by Abra- 
ham de Moivre (1667-1754) and James Stirling (1692-1770). Their work, together with 
work of Euler, Daniel Bernoulli (1700-1782) and Lagrange (1736-1813) is reviewed in 
Dutka [76]. 

In 1729, in correspondence with Christian Goldbach, Euler discussed the problem 
of obtaining continuous interpolations of discrete series. He described a (convergent) 
infinite product interpolating the factorial function, ml, 



. t^in 2^"™' • 3™ 3^""^ • 4"^ 41—^" . 5*^ 



1 + m 2 + m 3 + m 4 + m 



etc. 



'''"numerus n [...], et quern per approximationem olim inveni esse = 0,5772156649015325, cuius val- 
orem nuUo adhuc modo ad mensuras transcendentes iam cogitas redigere potui; unde haud inutile erit 
resolutionem huius formulae propositae pluribus modis tentare." [Translation: Jordan Bell] 
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This infinite product, grouped as shown, converges absolutely on the region C\{— 1, —2, ...}, 
to the gamma function, which we denote T{m + 1), following the later notation of Le- 
gendre. John Wallis (1616-1703) had earlier called this function the "hyper geometric 
series" . Euler substitutes m = ^ and finds 




by recognizing a relation of this infinite product to Wallis's infinite product for ^. In 
this letter he also introduced the notion of fractional integration, see Sandifer |262j . 

Also in 1729 Euler [55.1 E19] wrote a paper presenting these results, but it was not pub- 
lished until 1738. This paper describes the infinite products, suitable for interpolation. 
It then considers Euler 's first integral 

Jo 

which in modern terms defines the Beta function B[e + l,n + 1), derives a recurrence 
relation from it, and in Sect. 9 derives B{e + 1, n + 1) = j^^^^^p^y - More generally, in the 
later notation of Legendre this extends to the Beta function identity 

i?(r, .) := f x^-\l - xf-'dx = ^^p^. (2.3.1) 
Jo r(r + s) 

relating it to the gamma function. Euler \86\ Sect. 14] goes on to obtain a continuous 
interpolation of the factorial function f{n) = n!, given by the integral 

/(n)= [\-logxrdx, (2.3.2) 
Jo 

where n is to be interpreted as a continuous variable, n > —1. He then computes formulas 
for r (l + ^) at half-integer points as above [Ml Sect. 20]. 

In 1765, but not published until 1769, Euler \103\ E368] ("On a hypergeometric curve 
expressed by the equation y = l*2*3*---*x") studied the "hypergeometric curve" 
given (using Legendre's notation) by y = T{x + 1), where 

/•OD 

r(x) := / f'-^e-^dt. 
Jo 

Euler interpolates the factorial function by an infinite product. In Sect. 5 he gives the 
product formula 




with any fixed a with 1 < a < x, and with a = — ^ obtains 

^ \ 2 J ^}^\k + x\2k-l + xj J 
In Sect. 9 he derives a formula, which in modern notation reads 

logy = -jx + y;(-l)«^x", (2.3.3) 

n=2 



Here the letter e denotes an integer. Only in a later paper, [E71], written in 1734, does Euler use 
the symbol e to mean the constant 2.71828.... 
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where 7 is Euler's constant. In Sect. 11 he makes an exponential change of variables 
from his formula (j2.3.2p to obtain the integral formula 

/•oo 

y= / e-Vdi;, (2.3.4) 
Jo 

which is (essentially) the standard integral representation for the gamma function. This 
formula shows that the "hypergeometric curve" is indeed given aqj 

y = T{x+l). 

In Sect. 12 he obtains a formula for the digamma function ip{x) = j^{x), which is 

r'(x + 1) dy X 
r(x + 1) '~ yllx ~ ~^ ^ k{x + k)' 

In formula VII in this section he gives the asymptotic expansion of logr(x+l), (Stirling's 
series), expressed as 

1 , 1 , A 1-2B 1-2-3-4C l-2-3-4-5-6L> 

logy = -log27r+(x+-)logx-x+ — -^-3+ ^+^*^- 

(2.3.5) 

with A = ^, B = C = gig, D = ... Euler also gives here the value of Euler's 

constant as 7 w 0, 57721 56649 01422 5. In Sect. 14 he finds the explicit value 

r'(l) := I L=o = -7- (2.3.6) 
In Sect. 15 he uses r(|) = ^\/tt to derive the formula 

^Q)=-7 + 2-21og2. 

In Sect. 26 he observes that, for x = n a positive integer, and with a treated as a 
variable, 

00 

r(x + l) = J](-l)'=(^)(a-fc)^ (2.3.7) 

the series terminates at A; = n, and the resulting function is independent of a. This leads 
him to study in the remainder of the paper the series 

m(m — 1) , ,„ m(m — l)(m — 2) , 
s = - m(x - 1)" + \ x - 2)" + ^ '-{x - 3)" + &c 

In this formula he sets n = m + X, and he investigates the expansion s = s{x,X) as a 
function of two variables {\,x) which corresponds to the variables (x — n,a) in ()2.3.7p . 
These expansions involve Bernoulli numbers. 

2.4. Zeta values. Already as a young student of Johann Bernoulli (1667-1748) in the 
1720's Euler was shown problems on the summation of infinite series of arithmetical func- 
tions, including in particular the sums of inverse fc-th powers of integers. The problem of 
finding a closed form expression for this value, C(2), became a celebrated problem, called 
"The Basel Problem." The Basel problem had been raised by Pietro Mengoli (1628-1686) 
in his 1650 book on the arithmetical series arising in geometrical quadratures (see [210] . 



'^If one instead used the notation Tl{x) for this integral used by Gauss, cf. Edwards [781 p. 8], then 
the hypergeometric curve is y = 11 (x). Gauss made an extensive study of hypergeometric functions. 
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|123j ). In that book Mengoli [2121 Prop. 17] summed in closed form the series of (twice 
the) inverse triangular numbers 

^ 1 _ 1 1 ^ _^ 
^ n(n + 1) ~ 2 ^ 6 ^ 12 "' ~ ' 

n=l 

which is a telescoping sum. The problem became well known because it was also raised by 
Johann Bernoulli's older brother Jacob (also called James or Jacques) Bernoulli (1654- 
1705), earlier holder of the professorship of mathematics at the University of Basel. He 
wrote a book on infinite series ([29]) in 1689; and several additional works on them. A 
collection of this work appeared posthumously in 1713 [30j. He could tell that the sum 
of inverse squares converged to a value between 1 and 2, using the triangular numbers, 
via 

' 1 




n + l 



2. 



He was unable to relate this value to known constants, and wrote ([30l art. XVII, p. 
254]): 

And thus by this Proposition, the sums of series can be found when the 
denominators are either triangular numbers less a particular triangular 
number or square numbers less a particular square number; by [the results 
of art.] XV, this happens for pure triangular numbers, as in the series 
T~'^l~^k~^To~^TE ^^i^^ •^'^ other hand, when the numbers are 
pure squares, as in the series j + \ + ^ + jq + -^ &c., it is more difficult 
than one would have expected, which is noteworthy. If someone should 
succeed in finding what till now withstood our efforts and communicate 
it to us, we would be much obliged to them0 

Here he says that if a > 1 is an integer and the denominators range over n(n + 1) — 
a(a + 1) or alternatively — a^, with n > a + 1, then he can sum the series, but he 
cannot sum the second series when a = 0. 

Euler solved the Basel problem by 1735 with his famous result that 

Euler's first contributions on the Basel problem involved getting numerical approxima- 
tions of C(2)- In ^87i, E20], written in 1731, integrating his representation for the harmonic 
numbers led to a representation (in modern terms) 

1 r 1 - t" , ^1 
By a series acceleration method he obtained the formula 

oo ^ 

C(2) = (log2)2 + j;^^, 

n=l 







-'^'^"Atque ita per hanc Propositionem, inveniri possunt summae serierum, cum denominatores sunt vel 
numcri Trigonales minuti alio Trigonali, vel Quadrati minuti alio Quadrato; ut & per XV. quando sunt 
puri Trigonales, ut in serie y + | + | + ^ + ^ &c. at, quod notatu dignum, quando sunt puri Quadrati, 
ut in serie j + j + ^ + jn + 'i^ ifec difRcilior est, quam quis expectaverit, summae pervestigatio, quam 
tamen finitam esse, ex altera, qua manifesto minor est, coUigimus: Si quis inveniat nobisque communicet, 
quod industriam nostram elusit hactenus, magnas de nobis gratias feret." [Translation: Jordan Bell] 
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in which the final sum on the right is recognizable as the dilogarithm value 2 Li2(^). Euler 
used this formula to calculate C(2) ~ 1.64492 4, cf. Sandifer [262] . He next obtained 
a much better approximation using Euler summation (i.e. Euler-Maclaurin expansion). 
He announced his summation method in \88\ E25] in 1732, and gave details in |9H E47], 
a paper probably written in 1734. In this paper he computed Bernoulli numbers and 
polynomials, and also applied his summation formula to obtain ("(2) and C(3) to 20 
decimal places. 

In December 1735, in |89l E41], published in 1740, Euler gave three proofs that C{2) = 

2 

The most well known of these treated sinx as though it were a polynomial with 
infinitely many roots at ivr, ±2tt, ib37r, • • • . He combines these roots in pairs, and asserts 
the infinite product formula 

He then equates certain algebraic combinations of coefficients of the Taylor series expan- 
sion at X = on both sides of the equation. The right side coefficients give the power 
sums of the roots of the polynomials. Writing these Taylor coefficients as 1 — ax^ + /3x^ — 
7x^ + . . -jhe creates algebraic combinations giving the power sums of the roots, in effect 
obtaining^ 

= « 

^C(4) = a' -2/3 

1 Q 

—(^(6) = - 3al3 + 37, 
vr" 

and so on. Now the coefficients a,/3,7 may be evaluated as rational numbers from the 
power series terms of the expansion of which Euler computes directly. In this way 
he determined formulas for C{2n) for 1 < n < 6. He also evaluated by a similar approach 
alternating sums of odd powers, obtaining in Sect. 10 the result of Leibniz (1646-1716) 
that 



2A; + 1 4 

fc=0 

He views this as confirming his method, saying: 

And indeed this is the same series discovered some time ago by Leibniz, 
by which he defined the quadrature of the circle. From this, if our method 
should appear to some as not reliable enough, a great confirmation comes 
to light here; thus there should not be any doubt about the rest that will 
be derived from this method. 

In Sect. 12 he obtained a parallel result for cubes, 

Vl)^ I = ^ 

' (2k + 1)3 32 

k=o ^ ' 

These proofs are not strictly rigorous, because the infinite product expansion is not 
completely justified, but Euler makes some numerical cross-checks, and his formulas 
work. Euler was sure his answer was right from his previous numerical computations. In 



-'^^Euler only writes down the right sides of these equations. 

"'^■^"Atque haec est ipsa series a Leibnitio iam pridem prolata, qua circuh quadraturam definiuit. Ex 
quo magnum huius methodi, si cui forte ea non satis certa videatur, firmamentum elucet; ita ut de 
rehquis, quae ex hac methodo deriuabantur, omnino non hceat dubitari" [Trans: Jordan BeU[ 
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this same paper he gives two more derivations of the result, stiU susceptible to criticism, 
based on the use of infinite products. 

We note that directly equating power series coefficients on both sides of (j2.4.ip . for 
the the 2n-th coefficient would give the identity 

(-l)n _ (_l)n 1 

(2n + f)! vr^" ^ m^m^ • • • mi' 

^ ' l<mi<m2<---<mn ^ ^ 

The sum on the right hand side is the multiple zeta value C(2, 2, 2) (using n copies of 
2, see (|2.4.5p below), and the case n = 1 gives C(2)- Relations like this may have inspired 
Euler's later study of multiple zeta values. 

Euler communicated his original proof to his correspondents Johann Bernoulli, Daniel 
Bernoulli (1700-1782), Nicolaus Bernoulli (1687-1759) and others. Daniel Bernoulli 
criticized it on the grounds that sirnrx may have other complex roots, cf. [21] Sec. 4], 
[309\ p. 264]. In response to this criticism Euler eventually justified his infinite product 
expansion for sinvra;, obtaining a version that can be made rigorous in modern terms, cf. 
Sandifer [262] . 

In a paper [921 E61] published in 1743, Euler obtained more results using infinite 
products, indicating how to derive formulas for C(2?^), for all n > 1, as well as for the 
sums 

oo 

[L(2n + l,X-4) H E(-l) \2fc + 11^^+1 ' 
fc=o ^ ^ 

for n > 0. In this paper he studied the function = X^^q 7J- recalls his formula 
for sins, as 

^3 5^ 5^ 

^-rT2:3 + 1.2.3.4.5---- = '^'-^)('-4^)^'-9^)('-l6^)--- 

Then he obtains the formula sinz = — , introducing the modern notation for e, 
saying: 

For, this expression is equal to , where e denotes that number 

whose logarithm = 1, and 

with n being an infinitely large numbeiF^. 

Now he introduces a new principle. He evaluates the following integrals, where p and q 
are integers with < p < q: 

-ax 



and 

-1 ^p-l _ ^.g-p-l TTCOS ^ 

-dx " 



l-xi g sin 2^ 



Expanding the first integral as an indefinite integral in power series in x, evaluating term 
by term, and setting x = 1 he obtains the formula 

vr 1 1 1 1 1 1 1 

- + ^ + + ^ ^ + ••• 



q sin ^ P Q — P Q + P 2q — p ' 2q + p ' 3q — p 3q + p 



"'^'^"Haec enim expressio aequivalet isti ^ '^^^^ denotante e numerorum, cujus logarithmus est 
= 1, & cum sit = (l + ■!■)" existente n numero infinito, " [Translation; Jordan Bell] 
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and similarly for the second integral. Setting ^ = s, and rescaling, he obtains 

vr 1 1 1 1 1 1 

+ z rr-. + + ^ 



sin ITS s 1 — s 1 + s 2 — s 2 + s 
and 

VrCOSTTS 1 1^1 1^1 1^ 

sin ITS s 1 — s 1 + s 2 — s 2 + s 3 — s 
By differentiating with respect to s, he obtains 

vr^ cos vrs 1 1 1 



(sin7rs)2 s2 (l-s)2 (l + s)2 (2-s)2 (2 + s)2 {3-sf 

and 

7r2 1 1 1 1 1 

+ ... , ^o + 77^ ^ + ,^ . xo + 



(sin7rs)2 s2 (l-s)2 (l + s)2 (2-s)2 (2 + s)2 
Now he substitutes s = ^ with < p < q integers and obtains many identities. Thus 

1 1 1 1 1 



g2(sm^)2 p2 {q-py {q + pf {2q - p)^ {2q + p)^ 
For g = 4 and p = 1 he obtains in this way 

7r2 ^1 1 



8^2 32 52 ' 72 ' 92 112 132 

and 

7r2 _ 1 1 1 1 1 

T~ 32+52 + 72 + 92+Tl2+'" 

2 

He observes that it is not difficult to derive on like principles the value C{2) = (It 
follows for example from the observation that the series above is (1 - 2^)C(2).) He 
concludes that one may continue differentiating in s to obtain formulas for various series 
in n-th powers, and gives several formulas for the differentials. In particular, one can 
obtain formulas for C(2f^) by this method, which show that it will be a rational multiple 
of 7r2". However this approach does not make evident the closed formula for C(2ra) 
expressed in terms of Bernoulli numbers (given below as ()2.4.2p ). 

2 

In a paper [931 E63], published in 1743, Euler gave a new derivation of C(2) = 
using standard methods in calculus involving trigonometric integrals, which is easily 
made rigorous by today's standards. He takes s to be arclength on a circle of radius 1 
and sets x = sins so that s = arcsinx. He observes x = 1 corresponds to s = ^ and here 
he uses the symbol vr with its contemporary meaning : 

It is clear that I employ here the letter vr to denote the number of Ludolf 

van Ceit/en 3. 14159265.. E3 
Here Euler refers to Ludolf van Ceulen (1540-1610), who was a professor at Leiden 
University and who computed vr to 35 decimal places. In his calculations Euler works 
with differentials and writes ds = —^==, so that s = J —j==. Multiplying these 
expressions gives 

dx f dx 



sds 



yjl — XX J \/T 



XX 



"'^^"11 est clair que j'emplois la lettre vr pour marquer le nombre de LUDOLF a KEULEN 3,14159265 



etc." 
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He integrates both sides from a; = to x = 1. The left side is Jq^ sds = For the right 
side, he uses the binomial theorem 



1 



1 • 3 



:[=(!- xx)-2] = l + -x' + —X^ + 



1 • 3 • 5 



-x" + etc. 



^/l-xx' ' ' ' ' 2 '2-4 ' 2-4-6" 

In fact this expansion converges for |x| < 1. Euler integrates it term by term to get 

1-3-5 



/ 



dx 



X + 



1 



XX 



2-3 



X + 



1-3 
2-4-5 



X + 



2-4-6-7 



x' + etc. 



He then obtains 

xdx 



sds 



+ 



x^dx 



TTTT 



+ 



1 • 3 x^dx 



\/l — XX 2 • 3 ^1 — XX 2 • 4 • 5 yT^^"xx 
Integrating from x = to x = 1 would give 

x^dx 



+ 



1 • 3 • 5 x^dx 
2-4-6-7vT^r^ 



+ etc. 



+ 



1 



XX 



2 • 3 



+ 



XX 



1 • 3 
2-4-5 



x^dx 



etc. 



XX 



The individual terms on the right can be integrated by parts 



r x"+^dx n + 1 f x"dx x"+^ ^ 
y VI - a;x ~ n + 2 7 ^/l-xx n + 2 



XX. 



When integrating from x = to x = 1, the second term on the right vanishes and Euler 
gives a table 



f xdx 
Jo Vl - XX 

x^dx 
Jo Vl — XX 

x^dx 
Jo VT^ 

L 



= 1 — Vl — iCX = 1 



XX 



x^c/x 

XX 

x^dx 



5 70 
6 

7 Jo 



xdx 


2 






y/l — XX 


~ 3 






x^dx 


2 


4 




\/l — XX 


3 


5 




x^dx 


2 


4 


6 


Vl — XX 


3 


5 


7 


x'^dx 


2 


4 


6-8 


\/l — XX 


3 


5 


7-9 



\/l — XX 

Substituting these values in the integration above yields 



TTTT 

^ = 1 + 



1111 

+ + 1^ + TT^ + etc. 



3-3 5-5 7-7 9-9 



This gives the sum of reciprocals of odd squares, and Euler observes that multiplying by 



i gives 



from which follows 



1111 

-H 1 1 h etc. 

4 16 36 64 



TTTT , 1 1 1 



This proof is rigorous and unarguable. 

He goes on to evaluate sums of larger even powers by the same methods, and gives a 
table of explicit evaluations for C{2n) for 1 < ra < 13, in the form 



22n-l 
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In this expression, he says part of the formula is known and explainable, and that the 
remaining difficulty is to explain the nature of the fractions [ C2n] which take values of 
a different character: 

111 3 5 691 35 

2 6 6 10 6 210 Y 
He says that he has two other methods for finding these numbers, whose nature he does 
not specify. As we know, these involve Bernoulli numbers. 

In this same period, around 1737, not published until 1744, Euler |96, E72] obtained 
the "Euler product" expansion of the zeta function. In [951 Theorem 8] he proved: 

The expression formed from the sequence of prime numbers 

2" . 3" • 5" • 7^^ • 11" • etc. 



(2" - 1)(3" - 1)(5" - 1)(7" - 1)(11" - 1) etc. 
has the same value as the sum of the series 

1 1 1 1 1 1 . [T^ 

1 + 2^ + ^ + 4^7 + ^ + ^ + 7^7 + c 113 

This infinite product formula for Cl'^) supplies in principle another way to approximate 
the values Cin) for n > 2, using finite products. 

Euler made repeated attempts throughout his life to find closed forms for the odd 
zeta values, especially C(3)- In a 1739 paper [981 E130], not published until 1750, he 
developed partial fraction decompositions and a method equivalent to Abel summation. 
In particular, letting 9{s) = I^^q (2n+T)^ ~ Snii ^"ll^ — ' obtained the 

formula 

^(l-2„,)= '-"°-':,f"'-^" ^(2..) 

for m = 1, 2, 3, 4, in which the left side is determined as the Abel sum 

^ (-l)"-ix" 
(f>{m) := hm > ; 

n=l 

in fact the function (p{s) is Abel summable for all s S C. 

In a paper written in 1749 [102^ E352], not presented until 1761, and published in 
1768, Euler derived the functional equation for the zeta function at integer points, and 
he conjectured 

0(1 -s) _ -r(s)(2^ - l)cosf 
(/>(s) ~ (2'^-! - l)7r^ ■ 
He apparently hoped to use this functional relation to deduce a closed form for the odd 
zeta values, including C(3), and says: 

As far as the [alternating] sum of reciprocals of powers is concerned, 

1 1 1 1 1 „ 

1 \ \ h&c 

2" 3" 4" 5" 6" 

I have already observed the sum [</>(n)] can be assigned a value only when 

n is even and that when n is odd, all my efforts have been useless up to 

now 13 



"'^^"Si ex serie numerorum primorum sequens formetur espressio ( 2" - 1 ) ( 3"^- 1 x 5 " - 1) ( 7"^- i^l'ii'" - 1 ) etc 
erit eius aequalis summae huius serial l + ^ + -^ + ^ + -^ + -(pr + Y7^ + etc." [Translation: David 
Pengelley [239] . |240] .] Here n is an integer. 

^''"A regard des series reciproques des puissances 

1 _L 1_ J 1_ 

2" ^3" 4" + 5n 6" ^ 
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In his 1755 book on differential calculus Euler [100, E212, Part II, Cfiap. 6, Arts. 141- 
153] gave closed forms for CC^^) explicitly expressed in terms of the Bernoulli numbers. 
These are equivalent to the modern form 

C(2n) = (-l)"+^ ^^;gp . (2.4.2) 

He also used his summation formula to numerically evaluate the sums using divergent 
series, as follows. To estimate C(2) he considers the finite sums 



1 1 1 

4 9 



where x is fixed, and obtains the summation 



s = C-- + — -^ + ^-^ + ---, (2.4.3) 

X 2x^ x^ x^ X'' ' 

where the constant C is to be determined by finding all the other terms at one value 
of X. Choosing x = oo we have s = C{2) and the summation formula formally gives 
s = C = C(2). On setting x = 1, the finite sum gives s = 1, therefore one has, formally, 

C = 1 + 1- i + l^sl- 1^41 + 1^61 

The right side of this expression is a divergent series and Euler remarks: 

But this series alone does not give the value of C, since it diverges 
strongly. Above [§125] we demonstrated that the sum of the series to 

2 

infinity is and therefore setting x = oo and s = we have C = 
because then all other terms vanish. Thus it follows that 

l + l-l: + \B2\-\B,\ + \B,\---- = ^E 
2 b 

Euler's idea is to use these divergent series to obtain good numerical approximations to 

C by substituting a finite value of x and to add up the terms of the series "until it begins 

to diverge" , i.e. to truncate it according to some recipe. We do not address the question 

of what rule Euler used to determine where to truncate. This procedure has received 

modern justification in the theory of asymptotic series, cf. Hardy [153J . Hildebrand 

Sec. 5.8], [781 Sec 6.4]. By taking large x one can get arbitrarily accurate estimates 

(for some functions). If the series is alternating, one may sometimes estimate the error 

as being no larger than the smallest term. Euler demonstrates in this particular case 

how to find the value C = ^ numerically using the value x = 10, with s = Yli=i by 

truncating the expansion ()2.4.3p after power obtaining C 1, 64493 40668 48226 430. 

He writes: 

2 

And this number is at once the value of the expression as will be 
apparent by doing the calculation with the known value of vr. Whence 
it is at once understood that even though the series [I-B2I, I-B4I, [-Bel • • • ] 
diverges, still the true sum is produced]^. 



j'ai deja observe, que leurs sommes ne sauroient etre assignees que lorsque I'exposant n est un nombre 
entier par, & que pour les cas ou n est un nombre entier impair, tons mes soins ont ete jusqu'ici inutiles." 

^^"quae series autem cum sit maxime divergens, valorem constantis C[:— C(2)] non ostendit. Quia 
autem supra demonstavimus summam huius seriei in infinitum continuatac esse — facto a:: = 00, si 
ponatur s = fiet C = ob reliquos terminos omnes evanescentes. Erit ergo 1 + 1— i+il— 23 + 
e:-2) + £-&c=^." [Translation: David Pengelley [233 Sec. 148], [240] ] 

Hicque numerus simul est valor expressionis quemadmodem ex valore ipsius tv cogito calculum 
institutenti patebit. Unde simul intelligitur, etiamsi series 21, 58, £ & c. divergat, tamen hoc modo veram 
prodire summan." [Translation: Jordan Bell] 
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Euler also applies this method to the odd zeta values, obtaining numerical approxima- 

3 

tions for all zeta values up to C(15), in particular ^(3) 25 79436 • -^^ ^^^o (implicitly) 
introduces in Sec. 153 a "Bernoulli function" b{z) which interpolates at integer values 
the variable z his version of the even Bernoulli numbers, i.e. 

b{n) = \B2n\ = {-ir+^B2n, u e Z>i. 

Euler remarks: 

The series of Bernoulli numbers [I-B2I) I-S4I, l-Bel) " " " ] however irregular, 
seem to be interpolated from this source; that is, terms constituted in 
the middle of two of them can be defined: for if the middle term lying 
between the first A and the second B, corresponding to the index 1^ were 
= p, then it would certainly be 

1 1 2^p 



and hence 



3 11 
P = ^(1 + ^ + ^ + &c) = 0, 05815 2270 



That is, Euler shows the interpolating value p := 6(|) = p is related to C(3), by the 
relation 

3C(3) 

and he finds the numerical value p ~ 0.05815227. He similarly obtains a value for 6(|) 
in terms of C(5)- 

In a 1771 paper, published in 1776, Euler [1061 E477] ( "Meditations about a singular 
type of series ") introduced multiple zeta values and with them obtained another formula 
for C(3)- He had considered this topic much earlier. In Dec. 1742 Goldbach wrote Euler 
a letter about multiple zeta values, and they exchanged several further letters on their 
properties. The 1771 paper recalls this correspondence with Goldbach, and considers 
the sums, for m,n>l involving the n-harmonic numbers 

Sm,n ■= ^ ^ ( X] 1 ^ X^X] [tH^J ' (2-4.4) 

k=l k>l \l<l<k I l>\ k>l ^ ' 

Euler lets J denote C,(rn) and lets / ps:(^)" denote the sum Sm,n- He considers all 
values 1 < m, n < 8. This includes divergent series, since all sums are divergent. 
The modern definition of multiple zeta values, valid for positive integers (ai, a^) with 
oi > 2, is 

C(«l:«2,-,afc) = ^ ] a, - (2.4.5) 

ni>n2>--->nfc>0 ^ ^ 

In this notation ()2.4.4p expresses the identities (valid for m > 2,n > 1) 

Sm,n = Cin + m)+ C{m, n). 



"Ex hoc fonte series numerorum BernouUiarum y <s j [• • ■ ] quantumvis irregularis videatur inter- 
polari, seu termini in medio binorum quorumcunque constituti assignari poterunt; se enim terminus 
medium interiacens inter primum 21 and secundum 23, sue indici li, respondens fuerit = p; erit utique 

^ + W + W + = T^^^ ideoque p = ^{1 + ^ + ^ + kc) ^ 0, 05815227." [Translation: Jordan 
Bell] 
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Euler presents three methods for obtaining identities, and estabhshes a number of addi- 
tional identities for each range m + n = j for 2 < j < 11. For j = 3 he obtains 

52,1 = 2 C(3). 

This is equivalent to 

C(3) = C(2,l) = i(Ef). 

fc=l 

Being thorough, Euler also treats cases m = 1 where the series diverge, and observes 
that for j = 2 he cannot derive useful identities. 

In a 1772 paper [1051 E432] {^'Analytical Exercises"), published in 1773, Euler obtained 
new formulas for values Ci2n + 1) via divergent series related to the functional equation 
for the zeta function. As an example, in Sect. 6 he writes, treating n as a positive integer 
and oj as an infinitely small quantity, the formula 

11 —1 ^n+u) , , 

^ + gn+c + 5„+.. + 2cos(^^)l-2---(n-l+a;)l ^ 7' 

Here the right side contains a divergent series. In Sect. 7, he in effect takes a derivative 
in ijj and obtains, for example in the case n = 3, the formula 

1 + ^ + ^ + ^ + etc. = - y(2^ log2 - 3^ log3 + 4^ log4 - 5^ log5 + etc, 

The right side is another divergent series, which he names 

Z := 2^ log 2 - 3^ log 3 + 4^ log 4 - 5^ log 5 + 6^ log 6 - 7^ log 7 + etc. 

He applies various series acceleration techniques on this divergent series, compare the 
discussion of divergent series in §2.5 below. In Sect. 20 he transforms Z into a rapidly 
convergent series, 



^ _ 1 avr^ /^TT^ 77r6 ^vr^ e7r^° 

~ 4 ~ 3 • 4 • 22 ~ 5 • 6 • 24 ~ 7 • 8 • 2^ ~ 9 • 10 • 28 ~ 11 • 12 • 2^0 ~ ^ ^' 



In his notation ("(2) = avr^, = etc, so he thus obtains the valid, rapidly conver- 
gent formula 

7^,^, TT^n ^ C(2n) 



8-C(3) 2(4 (2n + l)(2n + 2)22^ 

After further changes of variable, he deduces in Sect. 21 the formula 

111 TT^ rf 

1 + ^ + ^ + ^ + = Y log 2 + 2 y ^ log(sin ^)d<l>. 

This may be rewritten 

7 vr^ [2 

-C(3) = — log 2 -|- 2 / X log(sinx)(ix. 

8 4 Jq 

Euler notes also the "near miss" that a variant of the last integral is exactly evaluable: 



/ log(sinx)dx = ——log 2. 
Jo 2 



Euler continued to consider zeta values. In the 1775 paper [1081 E597] ("A new and 
most easy method for summing series of reciprocals of powers" ) , published posthumously 
in 1785, he gave a new derivation of his formulas for the even zeta values Ci2n). In the 
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1779 paper \110\ E736], published in 1811, he studied the dilogarithm function Li2{x) 
Yl^=i f^- Sect. 4 he establishes the functional equation 



n 



vr 

Li2(x) + Li2(\ ~ ^) = Y ~ (loga;)(log(l - x)). 

2 

The special case x = of this functional equation gives again ^(2) = \. Euler deduces 
several other functional equations for the dilogarithm, and by specialization, a number 
of new series identities. 

2.5. Summing divergent series: the Euler-Gompertz constant. In the course 
of his investigations Euler encountered many divergent series. His work on divergent 
series was a by-product of his aim to obtain accurate numerical estimates for the sum 
of various convergent series, and he learned and developed summability methods which 
could be applied to both convergent and divergent series. By cross-checking predictions 
of various summability methods he developed useful methods to sum divergent series. 
Of particular interest here, Euler's summation of a particular divergent series, the Wal- 
lis hypergeometric series, produced an interesting new constant, the Euler-Gompertz 
constant. 

Euler's near contemporary James Stirling (1692-1770) developed methods of accelera- 
tion of convergence of series by differencing, starting in 1719, influenced by Newton's work 
([2291 PP- 92-101]) on interpolation, see Tweddle [302]. Euler obtained Stirling's 1730 
book ( ^^Differential methods and summation of series") ( |286j . translated in [303J) which 
discussed infinite series, summation, interpolation and quadrature, in which Proposition 
14 gives three examples of a method of accelerating convergence of series of hyperge- 
ometric type. On 8 June 1736 Euler wrote to Stirling, saying ( \303\ p. 229], [301', p. 
141]): 

I searched all over for your excellent book on the method of differences, a 
review of which I had seen in the Acta Lipenses, until I have achieved my 
desire. Now that I have read through it diligently, I am truly astonished 
at the great abundance of excellent methods in such a small volume, by 
which you show to sum slowly convergent series with ease and how to in- 
terpolate progressions which are very difficult to deal with. But especially 
pleasing to me was Prop. 14 of Part I, in which you present a method 
for summing so easily series whose law of progression is not even estab- 
lished, using only the relation of the last terms; certainly this method 
applies very widely and has the greatest use. But the demonstration of 
this proposition, which you seemed to have concealed from study, caused 
me immense difficulty, until at last with the greatest pleasure I obtained 
it from the things which had gone before... 

A key point is that Proposition 14 applies when the members a„ of the series satisfy linear 
recurrence relations whose coefficients are rational functions of the parameter n. This 
point was not made explicit by Stirling. It is known that there exist counterexamples to 
Proposition 14 otherwise (cf. Tweddle [3031 PP- 229-233]). The convergence acceleration 
method of Stirling is really a suite of methods of great flexibility, as is explained by Gosper 
Euler incorporated these ideas in his subsequent work! I. 



^'-'Gosper [1331 p. 122] says: "We will be taking up almost exactly where James Stirling left off, and 
had he been granted use of a symbolic mathematics system such as MIT MACSYMA, Stirling would 
probably have done most of this work before 1750." 

We do not attempt to untangle the exact influence of Stirling's work on Euler. In the rest of his 
1736 letter Euler explains his own Euler summation method for summing series. In a reply made 16 
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In a 1760 paper [101, E247] {"On divergent series") Euler formulated his views on 
the use and philosophy of assigning a finite value to some divergent series. This paper 
includes work Euler had done as early as 1746. In it Euler presents four different methods 
for summing the particular divergent series 

1 - 1 + 2 - 6 + 24 - 120 + • • • [= 0! - 1! + 2! - 3! + 4! - 5! + • • • ], (2.5.1) 

which he calls "Wallis' hypergeometric series," and shows they each give (approximately) 
the same answer, which he denotes A. Barbeau and Leah |26) give a translation of the 
initial, philosophical part of this paper, followed by a summary of the rest of its contents, 
and Barbeau [27] gives a detailed exposition of Euler's treatment of this example. 

The series (j2.5.ip is termed "hypergeometric" because each term in the series is mul- 
tiplied by a ratio which varies from term to term, rather than by a constant ratio as in 
a geometric series. Series of this type appeared in the Scholium to Proposition 190 in 
Wallis's 1656 book Arithmetica Infinitorum |308p^ . 

Euler's four methods were: 

(1) [lOl_, Sect. 13-16] The iteration of a series acceleration method involving a 
differencing operation applied to this series. Euler iterates four times, truncates 
the resulting divergent series and obtains an approximation to its "sum" A of 
A w 38015/65536 = 0.58006...; 

(2) |10H Sect. 17-18] The interpolation of solutions to the difference equation Pn+i = 
nPn + 1, with Pi = 1. Euler obtains the closed formula 

P„ = 1 + (n - 1) + (n - l)(n - 2) + (n - l)(n - 2)(n - 3) + • • • , 

which is a finite sum for each n > 1, and observes that substituting n = 
formally gives the Wallis series. He considers interpolation methods that rescale 
Pn by an invertible function f{x) to obtain values f{Pn), and then attempts to 
interpolate the value f{Po)- For the choice /(x) = ^ he obtains A ~ 0.6 and for 
f{x) = logx he obtains A 0.59966; 

(3) |10H Sect. 19-20] The construction and solution of an associated differential 
equation. Euler views the Wallis series as the special case x = 1 of the power 
series 

s{x) = 1 - 1! X + 2! - 3! + 4! x"^ - 5! + • • • . 

Although this series diverges for all nonzero x S C, he observes that s(x) formally 
satisfies a linear differential equation. This differential equation has a convergent 
solution given by an integral with variable endpoint x, and by specializing it to 
X = 1 he obtains an exact answer A. We discuss this method in detail below. 



April 1738 Stirling tells Euler that Colin Maclaurin will be publishing a book on Fluxions that has a 
similar summation formula that he found some time ago. Euler replied on 27 July 1738 that "since Mr. 
Maclaurin probably came upon his summation formula before me, and consequently deserves to be named 
as its first discoverer. For I found that theorem about four years ago, at which time I also described its 
proof and application in greater detail to our Academy." Euler also states that he independently found 
some summability methods in Stirling's book, before he knew of its existence f [301l Chap. 6]). 

^^This reference was noted by Carl Boehm, an editor of Euler's Collected Works, cf. [261 p. 157]. 
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(4) [101', Sec. 21-25] The determination of a continued fraction expansion for the 
series s{x) given by the integral in (3). This is 

s(x) = i (2.5.2) 

X 

1 + 



1 + 



2x 

1 + 



1 

1 + - 



3x 

1 + 



1 + 

This continued fraction is convergent at a; = 1. Using it, he computes a good 
numerical approximation, given in ()2.5.8p below, to the answer A. 

In the third of four summation methods, Euler considers the series to be the "value" 
at X = 1 of the asymptotic series 

oo 

s{x) r^Y,^-l)''n\x''+\ (2.5.3) 

n=0 

He notes this series satisfies the inhomogeneous hnear differential equation 

s'{x) + ^s{x) = -. (2.5.4) 

He writes down an explicit solution to the differential equation (j2.5.4p that vanishes at 
X = 0, namely 

1 r 1 1 

s{x) := J -e~~^dt. (2.5.5) 

The expansion of this integral near x = has the asymptotic expansion (j2.5.3p . Euler 
transforms the integral (|2.5.5p into other forms. Using the change of variable v = « 
he obtains 

s(x) = e(^i) r ' (2.5.6) 
Jo 1-logt; 

He then obtains as his proposed summation for the divergent series the value 

S := s(l) = e -e-'-dt = C . (2.5.7) 

^ ' Jo t Jo 1-log?; ^ ' 

He numerically estimates this integral using the trapezoidal rule and obtains 5 ~ 0.59637255. 
Using his fourth method, in Sec. 25 Euler obtained a more precise numerical value for 
this constant, finding 

6 ^ 0.59634 73621 237. (2.5.8) 

The constant 5 defined by (j2.5.7p is of some interest and we call it the Euler- Gompertz 
constant, following Finch [113^ Sect. 6.2]. This constant was earlier named the Gompertz 
constant hy Le Lionnais [198]. El Benjamin Gompertz (1779-1865) was a mathematician 
and actuary who in 1825 proposed a functional model for approximating mortality tables, 
cf. Gompertz |129| . [T30] . The Gompertz function and Gompertz distribution are named 
in his honor. I have not located the Euler-Gompertz constant in his writings. However 
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Le Lionnais [198) gives no explanation or reference for his choice of name "Gompertz" 
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Finch [113J Sect. 6.2.4] points out a connection of the Euler-Gompertz constant with a 
Gompertz-hke functional model. 

The Euler-Gompertz constant has some properties analogous to Euler's constant. The 
integral representations for the Euler-Gompertz constant given by (|2.5.5p and (j2.5.6p 
evaluated at s = 1 resemble the integral representations (|2.2.9p and (j2.2.10p for Euler's 
constant. An elegant parallel was noted by Aptekarev [13j, 

/"OO /'OO 

7 = - / logxe"^dx, ^= log{x + l)e~''dx. (2.5.9) 
Jo Jo 

There are also integral formulae in which both constants appear, see (|3.5.1ip and ()3.6.8p . 

In 1949 G. H. Hardy [153', p. 26] gave a third integral representation for the function 

s(x), obtained by using the change of variable t = jq^, namely 



s{x) = X dw. 

Jq l + xw 

Evaluation at x = 1 yields the integral representation 

dw. (2.5.10) 



1 + w 

From this representation Hardy obtained a relation of S to Euler's constant given by 

« = -«(7-l + ^-ji^ + j^ + -). (2^5,11) 

One can also deduce from (12.5. lOp that 

S = -eEi{-l) (2.5.12) 
where Ei{x) is the exponential integral 

Ei{x) := I —dt, (2.5.13) 

J —oo ^ 

which is defined as a function of a complex variable x on the complex plane cut along 
the nonnegative real axis0 

Here we note that Euler's integral (j2.5.5p is transformed by the change of variable 
t = i to another simple form, which on replacing x by 1/x becomes 

1 f'°° p~" 

s(-) = eM —du = e''Ei{x). (2.5.14) 
X u 

Here Ei{z) is the principal exponential integral, defined for |arg2;| < tt by a contour 
integral 

roo -u 

Ei{z) := / du, (2.5.15) 

Jz U 

that eventually goes to +oo along the positive real axis, see [T, Chap. 5]. Choosing x = 1 
yields the identity 

Ei{l) = / dt = -. (2.5.16) 

./i t e 



^^This is sufficient to define Ei{—1). For applications the function Ei{x) is usually extended to positive 
real values by using the Cauchy principal value at the point x = where the integrand has a singularity. 
The real- valued definition using the Cauchy principal value is related to the complex-analytic extension 
by 

Ei{x) := lim - ( Ei{x + ei) + Ei{x - ei)) . 

e^0+ 2 V / 
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Further formulas involving the Euler-Gompertz constant are given in ^3.51 ^3.61 ^3.151 
and Em 

Regarding method (4) above, Euler had previously done much work on continued 
fractions in the late 1730's in [9H E71] and in particular in his sequel paper |97i E123], 
in which he expands some analytic functions and integrals in continued fractions. This 
work of Euler is discussed in the book of Khrushchev [171| . which includes an English 
translation of [97J. In 1879 E. Laguerre |188j found an analytic continued fraction expan- 
sion for the principal exponential integral, which coincides with e~^s(l/x)by ()2.5.14p . 
obtaining 



— dt = . 2.5.17 

t 1 ^ ' 

x + 1 

x + 3 



1 

x+5 4 

4 1 

x±7 9_ 

9 x±9 _ 
16 



Laguerre noted that this expansion converges for positive real x. His expansion simplifies 

OO g-t g— a; 

dt = ^ , (2.5.18) 

^ + 1 ^2 

X + 3 -75 



x + 5 



x + 7 ■ 



x + 9 

which is the form in which it appears in Wall \307\ (92.7)]. On taking x = 1 it yields 
another continued fraction for the Euler-Gompertz constant, 

OO -t I 

dt = . . (2.5.19) 

t l2 ^ ^ 



6 



42 



10 

In his great 1894 work T. J. Stieltjes ([284j. |285j ) developed a general theory for 
analytic continued fractions representing functions F{x) given by '^x+u ' '^hich takes 
the form 

d^{u) 1 



f 

Jo 



X + u 1 
aix H 



(2.5.20) 



1 

a2 H 

asx H 

04 H 



a^x + 



a6 + 



^^This form of fraction with minus signs between convergents is often called a J -continued fraction, 
see Wall [3011 P- 103]. 
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for positive real an- He formulated the (Stieltjes) moment problem for measures on the 
nonnegative real axis, and related his continued fractions to the moment problem. In 
the introduction of his paper he specifically gives Laguerre's continued fraction (j2.5.18p 
as an example, and observes that it may be converted to a continued fraction in his form 
()2.5.20p with a2n-i = 1 and a2n = On taking x = 1 his converted fraction is closely 
related to Euler's continued fraction (j2.5.2p . Stieltjes notes that his theory extends the 
provable convergence region for Laguerre's continued fraction (|2.5.18p to all complex x, 
excluding the negative real axis. In [284^ Sect. 57] he remarks that associated measure 
in the moment problem for this continued fraction must be unbounded, and in |284^ 
Sect. 62] he obtains a generalization of (his form of) Laguerre's continued fraction to 
add extra parameters; cf. Wall |307l eqn. (92.17)]. 

It is not known whether the Euler-Gompertz constant is rational or irrational. How- 
ever it has recently been established that at least one of 7 and d is transcendental, as 
presented in §3.161 

2.6. Euler-Mascheroni constant; Euler's approach to research. Euler's constant 
is often called the Euler-Mascheroni constant, after later work of Lorenzo Mascheroni 
[208j, [209j . In 1790, in a book proposing to answer some problems raised by Euler, 
Mascheroni [208] introduced the notation A for this constant. In 1792 Mascheroni [ 2091 
p. 11] also reconsidered integrals in Euler's paper [109, E629], using both the notations A 
and a for Euler's constant. In his book, using an expansion of the logarithmic integral, 
he gave a numerical expression for it to 32 decimal places ([2081 P- 23]), which stated 

A 0.57721 56649 01532 86061 81120 90082 39 

This computation attached his name to the problem; Euler had previously computed it 
accurately to 15 places. At that time constants were often named after a person who 
had done the labor of computing them to the most digits. However in 1809 Johann 
von Soldner [270| . while preparing the first book on the logarithmic integral function 
li(^x) = J found it necessary to recompute Euler's constant, which he called H 
(perhaps from the harmonic series). He computed it to 22 decimal places by a similar 
method, and obtained ( [270i p. 13]) 

H w 0.57721 56649 01532 86060 65 

Von Soldner's answer disagreed with Mascheroni's in the 20-th decimal place. To settle 
the conflict von Soldner asked the aid of C. F. Gauss, and Gauss engaged the 19-year old 
calculating prodigy F. G. B. Nicolai (1793 — 1846) to recheck the computation. Nicolai 
used the Euler-Maclaurin summation formula to compute 7 to 40 places, finding agree- 
ment with von Soldner's calculations to 22 places (see [1551 p. 89]). Thus the conflict was 
settled in von Soldner's favor, and Mascheroni's value is accurate to only 19 places. In 
the current era, which credits the mathematical contributions, it seems most appropriate 
to name the constant after Euler alone. 

The now-standard notation 7 for Euler's constant arose subsequent to the work of both 
Euler and Mascheroni. Euler himself used various notations for his constant, including 
C and O and n; Mascheroni used A and a; von Soldner used H and C. The notation 
7 appears in an 1837 paper on the logarithmic integral of Bretschneider [441 p. 260]. 
It also appears in a calculus textbook of Augustus de Morgan |22H p. 578], published 
about the same time. This choice of symbol 7 is perhaps based on the association of 
Euler's constant with the gamma function. 

To conclude our discussion of Euler's work, it seems useful to contemplate the approach 
to research used by Euler in his long successful career. C. Truesdell [298^ Essay 10], [299^ 
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pp. 91-92], makes the following observations about the methods used by Euler, his 
teacher Johann Bernoulli, and Johann's teacher and brother Jacob Bernoulli. 

(1) Always attack a special problem. If possible solve the special problem in a way 
that leads to a general method. 

(2) Read and digest every earlier attempt at a theory of the phenomenon in question. 

(3) Let a key problem solved be a father to a key problem posed. The new problem 
finds its place on the structure provided by the solution of the old; its solution 
in turn will provide further structure. 

(4) If two special problems solved seem cognate, try to unite them in a general 
scheme. To do so, set aside the differences, and try to build a structure on the 
common features. 

(5) Never rest content with an imperfect or incomplete argument. If you cannot 
complete and perfect it yourself, lay bare its flaws for others to see. 

(6) Never abandon a problem you have solved. There are always better ways. Keep 
searching for them, for they lead to a fuller understanding. While broadening, 
deepen and simplify. 

Truesdell speaks here about Euler's work in the area of foundations of mechanics, but his 
observations apply equally well to Euler's work in number theory and analysis discussed 
above. 

3. Mathematical Developments 
We now consider a collection of mathematical topics in which Euler's constant appears. 

3.1. Euler's constant and the gamma function. Euler's constant appears in several 
places in connection with the gamma function. Already Euler noted a basic occurrence 
of Euler's constant with the gamma function, 

r'(i) = -7. (3.1.1) 

Thus 7 appears in the Taylor series expansion of T{z) around the point z = 1. Further- 
more Euler's constant also appears (in the form e^) in the Hadamard product expansion 

for the entire function which states: 

r(2) ' 

Euler's constant seems however more tightly associated to the logarithm of the gamma 
function, which was studied by Euler [1031 E368] and many subsequent workers. The 
function logr(2;) is multivalued under analytic continuation, and many results are more 
elegantly expressed in terms of the digamma function, defined by 

«='^=TW = 5;'°«rW^ (3,L2) 

The theory of the digamma function is simpler in some respects because it is a mero- 
morphic function. Nielsen |234l Chap. 1] presents results and some early history, and 
Ferraro [112j describes the work of Gauss on these functions. 

We recall a few basic facts. The digamma function satisfies the difference equation 



(3.1.3) 
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which is inherited from the functional equation zT{z) = T{z -\- 1). It has an integral 
representation, valid for Re{z) > 0, 

oo /g— i ^—zt 



found in 1813 by Gauss [12H Art. 30 ff] as part of his investigation of the hyper geometric 
function 2-^1 (a, /3; 7; a;) (see |112j ). The difference equation gives 

n—l ^ 

i^{x) = il^ix + n) - Y] — — 

from which follows 

n— 1 ^ 

V'(x) = lim ( logn — > 

This limit is valid for all x G C \ Z<o, using the fact that ^(x) — log 2; goes uniformly to 
in the right half-plane as Re{x) — t- 00 ( |234l p. 85, 93]). This formula gives '0(1) = ~7; 
a fact which also follows directly from (j3.1.ip . 

Euler's constant appears in the values of the digamma function ^(x) at all positive 
integers and half- integers. In 1765 Euler [Ml E368] found the values r(i),r(l) and 
r'(i),r'(l). These values determines ■0(1/2), V'(l)) ™d then, by using the difference 
equation ()3.1.3p . the values of ip{x) for all half- integers. 

Theorem 3.1.1. (Euler 1765) The digamma function ip{x) := satisfies, for inte- 

gers n > 1, 

7p{n) = -7 + Hn~i 

using the convention that Hq = 0. It also satisfies, for half integers n + ^, and n> 0, 

i^{n + ^) = -7 - 21og2 + 2H2n-i - Hn-i, 
using the additional convention H-i = 0. 

Proof. The first formula follows from 0(1) = —7 using ()3.1.3p . The second formula will 
follow from the identity 

V^(^) = -7-21og2, (3.1.4) 

because the difference equation (|3.1.3p then gives 

Ix .Ix 1 1 

V'(n + -) = ^'(tt) + 2 1 + - + ••• + 



2' ""2' \ 3 2n-l, 
To obtain (j3.1.4p we start from the duplication formula for the gamma function 

T{2z) = ^2^^-'T{z)nz + \). (3.1.5) 
Jtt 2 



We logarithmically differentiate this formula at z = |, and solve for ^(^). □ 

The next two results present expansions of the digamma function that connect it with 
integer zeta values. The first one concerns its Taylor series expansion around the point 
z = l. 



30 



JEFFREY C. LAG ARIAS 



Theorem 3.1.2. (Euler 1765) The digamma function has Taylor series expansion 
around z = 1 given by 

oo 

^l^{z + l) = -j + J2{-l)''+\{k + 1) z'' (3.1.6) 

k=l 

This expansion converges absolutely for \z\ < 1. 

Remark. The pattern of Taylor coefficients in ()3.1.6p has the k = term formally 
assigning to "C(l)" the value 7. 

Proof. Euler [103^ E368, Sect. 9] obtained a Taylor series expansion around z = 1 for 
logr(2:), given as (j2.3.3p (with x = z — 1). Differentiating it term by term with respect 
to X yields (f3X6]l . □ 

The second relation to zeta values concerns the asymptotic expansion of the digamma 
function around z = +00, valid in the right half-plane. In this formula the non-positive 
integer values of the zeta function appear. 

Theorem 3.1.3. The digamma function ip{z) has the asymptotic expansion 

^(z + l)~logz + J^(-l)^C(l-A;)(^-J , (3.1.7) 

which is valid on the sector —it + e < aigz < vr — e, in the half-plane Re{s) > 0, for any 
given e > 0. This asymptotic expansion does not converge for any value of z. 

The statement that (j3.1.7p is an asymptotic expansion means: For each finite n > 1, 
the estimate 

k 



i;iz + 1) = log z + - ^) [-) + 



holds on the sector — vr + e < arg z < vr — e, with e > where the implied constant in the 
O-symbol depends on both n and e. 

Theorem 13.1.31 does not seem to have been previously stated in this form, with coef- 
ficients containing ({1 — k). However it is quite an old result, stated in the alternate 
form 



^(. + l)~log.-^-f;||-L. (3.1.8) 

fc=2 

The conversion from (j3.1.8p to (j3.1.7p is simply a matter of substituting the formulas for 
the zeta values at negative integers: ({0) = — ^ = Bi and ({1 — k) = (— l)'^~^^^||r, noting 
B2k+i = for A; > 1, see [H Sec. 1.18], [B (6.3.18)]. As a formal expansion (l3X8]) 
follows either from Stirling's expansion of sums of logarithms obtained in 1730 ( [286^ 
Prop. 28], c.f. Tweddle |301l Chap. 1]), or from Euler's expansion of logr(z) given 
in (I2.3.5p . by differentiation term-by-term. Its rigorous interpretation as an asymptotic 
expansion having a remainder term estimate valid on the large region given above is due 
to Stieltjes. In his 1886 French doctoral thesis Stieltjes [282J treated a general notion 
of an asymptotic expansion with remainder term, including the case of an asymptotic 
expansion for logr(a) ([2821 Sec. 14-19]). In 1889 Stieltjes [283] presented a detailed 
derivation of an asymptotic series for logr(a), whose methods can be used to derive the 
result (|3.1.8p above. 
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Proof. (1) The 1889 formula of Stieltjes [283J for logr(a) embodys Stirling's expansion 
with an error estimate. Stieltjes started with the function J{a) defined by 

logr(a) = (a - ^)loga - a + ^log(27r) + J(a). 

He notes that J(a) has an integral representation 

/■°° P(x) 

J(a)= / ^^-^dx, (3.1.9) 
^ ^ Jo x + a ' ^ ^ 

with periodic numerator 

P{x) = i - X + [xj = -Bi{x - lx\), 

in which Bi{x) is the first Bernoulli polynomial, and observes that this integral formula 
has the merit that it converges on the complex a-plane aside from the negative real axis, 
unlike two earlier integral formulas for J (a) of Binet, 

roo ]^ ]^ -ax 

and 

1 f°° a , / 1 



•^(") = ~/ — 2^°^[-, ^2^1 )^^' 

vr Jo a + \ 1 — e ^™ / 

which are valid on the half-plane Re{a) > 0. He now repeatedly integrates (|3.1.9p by 

parts to obtain a formula ri for logr(a), valid on the entire complex plane aside from 

the negative real axis, 

logr(a) = (a--)loga-a + -log(2vr) + ^-^ + ---+ ^^^_^^^^jti +'^fe(«)- 

(3.1.10) 

in which Jk{a) = j^^^^^^+rdx where Pk{x) is periodic and given by its Fourier series 

00 ^ 

P,{x) = {-l)H2k)\ sin(2vrnx), 

and is related to the {2k + l)-st Bernoulli polynomial. He obtains an error estimate for 
Jfe(a) ([2831 eqn. (34)]), valid for -vr + e < 6* < vr - e, that 

(2A: + 2)(2fc + l)^'''2^^ ^ • 

We obtain a formula for ^{a) by differentiating (j3.1.10p . and using i?2fc = (~1)^~^|-B2fc| 
gives 

V^(a) = loga---^-^ 2^ + '^^(")- 

This formula gives the required asymptotic expansion for '(/'(o)) ^ soon as one obtains 
a suitable error estimate for J).{a), which may be done on the same lines as for Jk{a)- 
Formula (j3.1.8p follows after using ^/^(a + 1) = ip{a) + ^, which has the effect to reverse 
the sign of the coefficient of -. 

Finally, the known growth rate for even Bernoulli numbers, 

\B2n\ ~ (^(2^)! ~ 4V^(7re)-2"n2". 



^^In Stieltjes's formula Bj denotes the Bernoulli number |i32j| in (|2.2.6|l . We have altered his notation 
to match (|2.2.6p . 
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implies that the asymptotic series (j3.1.7p diverges for all finite values of z. 

(2) The formula ()3.1.7p follows by direct substitution of the zeta values C(0) = —\ = B\ 

and C(l -k) = (-1)''+^^ for fc > 2, noting fiafc+i = for /c > 1, see pi Sec. 1.18], [U 
(6.3.18)]. □ 

A consequence of the asymptotic expansion (j3.1.7p for ^(x) is an asymptotic expansion 
for the harmonic numbers, given by 

oo ^ 

i/n ~ logn + 7 + 5](-l)'C(l - A;)^. (3.1.11) 

k=\ 

This is obtained on taking z = n in (|3.1.6p and substituting ■ip{n + 1) = Hn — 7. Sub- 
stituting the zeta values allows this expansion to be expressed in the more well-known 
form 

k=2 

J ^_ 1 _ 1 

log n + 7 + 2^ 12n2 + 120n4 252^6 + ' ' ' 
Note that although Hn has jumps of size ^ at each integer value, nevertheless the as- 
ymptotic expansion is valid to all orders (^)'^ for A; > 1, although the higher order terms 
are of magnitude much smaller than the jumps. 

Theorem 13.1.21 and Theorem 13.1.31 (2) have the feature that these two expansions 
exhibit between them the full set of integer values C{k) of the Riemann zeta function, 
with Euler's constant appearing as a substitute for "(^(1)". This is a remarkable fact. 

3.2. Euler's constant and the zeta function. The function C(s) was studied at 
length by Euler, but is named after Riemann |254j . commemorating his epoch-making 
1859 paper. Euler's constant appears in the Laurent series expansion of Cis) around the 
point s = 1, as follows. 

Theorem 3.2.1. (Stieltjes 1885) The Riemann zeta function C(s) is meromorphic in 
the plane with only singularity a simple pole at s = 1. Its Laurent expansion around 
s = 1 has the form 

C(s) = 7 + 70 + E -^^n{s - IT (3.2.1) 

s — 1 ^-^ n! 

n=l 

in which 70 = 7 is Euler 's constant, and the coefficients 7„ for n > are defined by 

(logfc)" (logmr+i\ 

■= JiToo I^A. —k ^Tl- j ^'-'-'^ 

Remark. Theorem 13.2.11 gives the coefficients 7„ of the zeta function in a form resem- 
bling the original definition of Euler's constant. These constants arose in a correspon- 
dence that Stieltjes carried on with his mentor Hermite (cf. [23] ) from 1882 to 1894. In 
Letter 71, written in June 1885, Stieltjes stated that he proposed to calculate the first 
few terms in the Laurent expansiorF^ of C(s) around s = 1, 

({z + 1) = - + ^0 + ^12: + ••• 
z 

noting that Aq = 0.577215665... is known and that he finds Ai = -0.072815520... In 
Letters 73 and 74, Hermite formulated (j3.2.2p . which he called "very interesting," and 



In Letter 71 Stieltjes uses the notation Co,Ci, but in Letter 75 he uses the same notation with 
a different meaning. We substitute Ao,Ai,... for the Letter 71 coefficients following [36] . 
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inquired about rigorously establishing itj^ In Letter 75 Stieltjes [23, pp. 151-155] gave 
a complete derivation of (|3.2.2p . writing 

The Laurent series coefficients An are now named the Stieltjes constants by some authors 
who follow Letter 71 ([36], [39]), while other authors follow Letter 75 and call the scaled 
values 7„ = C„ the Stieltjes constants ( jl69j . [57j). 

Proof. We follow a proof of Bohman and Froberg . It starts from the identity 

(s — l)C(s) = YlT=ii^ ~ which is valid for Re{s) > 1. For real s > 1, one has the 

telescoping sum 

oo 

k=l 

Using it, one obtains 

oo 

(s - l)C(s) = 1 + ^{(A; + ly-' - k"-' + {s- l)k-']. 

k=l 

Noting that at s = 1 the telescoping sum above is 0, one has 

oo 

(s-l)C(s) = l + ^{exp(-(s-l)log(A: + l))-exp(-(s-l)logA:) 
k=l 

+(s-l)iexp(-(s-l)log k)} 



= i+E{E^^7^ — ^[(iog(fc + i))"-(iogfe)"] 

k=l n=0 

s-l^ (-l)"(s - l)"(log A;)" -1 
^~A^^ ^. /■ 

n=0 

Dividing by (s — 1) now yields 



s — 1 ^-^ n\ 



n=0 

with 7o = 7, and with 

_ f (log kY (log(fc + 1))"+^ - (log fc)" + l -| 

k=\ 

This last formula is equivalent to ()3.2.2p . since the first m terms in the sum add up to 
the m-th term in (|3.2.2p . □ 

Table [J presents values of 7^ for small n as computed by Keiper jl69 ] • A good deal 
is known about the size of these constants. There are infinitely many 7„ of each sign 
(Briggs [45j). Because C(s) — is an entire function, the Laurent coefficients must 
rapidly approach in absolute value as n — )• 00. However despite the sizes of the initial 
values in Table [H the scaled values 7^ are known to be unbounded (and to sometimes 
be exponentially large) as a function of n, cf. Cohen ^9} Sec. 10.3.5]. Keiper estimated 
that 7150 ~ 8.02 x 10^^. More detailed asymptotic bounds on the 7^ are given in Coffey 



^^The result (13.2.21) is misstated in Erdelyi et al [H (1.12.17)], making the right side equal 
a mistake tracing back to a paper of G. H. Hardy [151| . 
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n 


7n 





+ 0.57721 56649 01532 


1 


- 0.07281 58454 83676 


2 


- 0.00969 03631 28723 


3 


+ 0.00205 38344 20303 


4 


+ 0.00232 53700 65467 


5 


+ 0.00079 33238 17301 


6 


- 0.00023 87693 45430 


7 


- 0.00052 72895 67057 


8 


- 0.00035 21233 53803 


9 


- 0.00003 43947 74418 


10 


+ 0.00020 53328 14909 



Table 1. Values for Stieltjes constants 7, 



|57j and Knessl and Coffey |172j. At present nothing is known about the irrationahty or 
transcendence of any of the Stieltjes constants. 

Various generalizations of the Riemann zeta function play a role in theoretical physics, 
especially in quantum field theory, in connection with a method of "zeta function regu- 
larization" for assigning finite values in certain perturbation calculations. This method 
was suggested by S. Hawking |156j in 1977, and it has now been widely extended to other 
questions in physics, see the book of Elizalde [79] (See also [80], [81].) Here one forms 
a generalized "zeta function" from spectral data and evaluates it an an integer point, 
often s = 0. If the function is meromorphic at this point, then the constant term in the 
Laurent expansion at this value is taken as the "regularized" value. This procedure can 
be viewed as a quite special kind of "dimensional regularization," in which perturbation 
calculations are made viewing the dimension of space-time as a complex variable s, and 
at the end of the calculation one specializes the complex variable to the integer value of 
the dimension, often s = 4, see Leibbrandt |197j . Under this convention Euler's constant 
7 functions as a "regularized" value "C(l)" of the zeta function at s = 1, because by 
Theorem 13.2.11 it is the constant term in the Laurent expansion of the Riemann zeta 
function at s = 1. This interpretation as "C(l)" also matches the constant term in the 
Taylor expansion of the digamma function ^|J{z) around z = 1 in Theorem 13. 1.2[ 

Euler's constant occurs in many deeper ways in connection with the behavior of the 
Riemann zeta function. The most tantalizing open problem about the Riemann zeta 
function is the Riemann hypothesis, which asserts that the non-real complex zeros of 
the Riemann zeta function all fall on the line Re(s) = ^. Sections 13.71 and 13.81 present 
formulations of the Riemann hypothesis in terms of the relation of limiting asymptotic 
behaviors of arithmetic functions to e"'. Section [3.91 presents results relating the extreme 
value distribution of the Riemann zeta function on the line Re(s) = 1 to the constant 

In 1885 Stieltjes announced a proof of the Riemann hypothesis ( |281] . English trans- 
lation in \285\ p. 561]). In Letter 71 to Hermite discussed above, he enclosed an an- 
nouncement of a proof of the Riemann hypothesis, requesting that it be communicated 
to Comptes Rendus, if Hermite approved, and Hermite did so. In Letter 79, written 
later in 1885, Stieltjes sketched his approach, saying that his proof was very complicated 
and that he hoped to simplify it. Stieltjes considered the Mertens function 

n 

M{n) = ^f,{j), 
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in which is the Mobius function (defined in ^3.3p and asserted that he could prove 
that \M{n)\/^/n is a bounded function, a result which would certainly imply the Riemann 
hypothesis. Stieltjes had further correspondence in 1887 with Mittag-Leffler, also making 
the assertion that he could prove ^J^^ is bounded between two limits, see [231 Tome II, 

Appendix], |253[ Sect. 2]. He additionally made numerical calculation^^ of the Mertens 
function M{n) for ranges 1 < n < 1200, 2000 < n < 2100, and 6000 < n < 7000, 
verifying it was small over these ranges. However he never produced a written proof 
of the Riemann hypothesis related to this announcement. Modern work relating the 
Riemann zeta zeros and random matrix theory now indicates that \M{n)\/^/n should 
be an unbounded function, although this is unproved, see te Riele [2531 Sect. 4] and 
Edwards [THl Sect. 12.1]. 

3.3. Euler's constant and prime numbers. Euler's constant appears in various sta- 
tistics associated to prime numbers. In 1874 Franz Mertens |214) determined the growth 
behavior of the sum of reciprocals of prime numbers, as follows. 

Theorem 3.3.1. (Mertens's Sum Theorem 1874) For x >2, 

log log X + B + R{x) , 



p<x 



with B being a constant 



B = j+y" /i(m)^^^^^ « 0.26149 72129- •• (3.3.1) 

' m. 



m 

m=2 



and with remainder term bounded by 



log(P + l) PlogP 
where P denotes the largest prime smaller than x. 

Here /u(n) is the Mobius function, which takes the values /x(l) = 1 and /u(n) = (— l)'^ 
if n is a product of k distinct prime factors, and /x(n) = otherwise. This result asserts 
the existence of the constant 



iim I — 

\p<x ^ 



B := lim I > log log x 



which is then given explicitly by (13.3.ip . Mertens's proof deduced that B = j — H with 
H :=- J]]{log(l --) + -} = 0.31571 84519 .... (3.3.2) 



p p 
p 



For modern derivations which determine the constant B see Hardy and Wright |1541 
Theorem 428] or Tenenbaum and Mendes Prance [2891 Sect. 1.9]. 

Mertens |214[ Sect. 3] deduced from this result the following product theorem, involv- 
ing Euler's constant alone. 

Theorem 3.3.2. (Mertens's Product Theorem 1874) One has 

lim (logj;) TT f 1 - - I = e"^. (3.3.3) 

p<x ^ ^ ' 



^^Such calculations were found among his posthumous possessions, according to te Riele [2531 p. 69]. 
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More precisely, for x > 2, 



in which 



n 



1 



p J log X 



4 2 1 

l*^^^^' - log(P + 1) + Pb^ + 2P' 
where P denotes the largest prime smaller than x. Thus S{x) — )• as x —)• cxd. 

In this formula the constant e'"^ encodes a property of the entire ensemble of primes 
in the following sense: if one could vary a single prime in the product above, holding the 
remainder fixed, then the limit on the right side of (j3.3.3p would change. 

In Mertens's product theorem the finite product 

^(-) - n (l - ^) 

is the inverse of the partial Euler product for the Riemann zeta function, taken over 
all primes below x, evaluated at the point s = 1. The number D{x) gives the limiting 
density of integers having no prime factor smaller than x, taken on the interval 
and letting T — t- oo while holding x fixed. Now the set of integers below T = x"^ having 
no prime smaller than x consists of exactly the primes between x and x^. Letting 7r(x) 
count the number of primes below x, the prime number theorem gives 7r(x) ~ jj^; and 
this yields the asymptotic formula. 

2. / N a;^ X 1 x^ 

7r(x j — 7r(x) ~ ^ — ~ — . 

log x^ log X 2 log X 

On the other hand, if the density D{x) in Merten's theorem were asymptotically correct 
already at T = x^ then it would predict the number of such primes to be about e~'^T^ — . 

log X 

The fact that e~"' = 0.561145... does not equal ^, but is slightly larger, exhibits the 
failure of the inclusion-exclusion formula to give an asymptotic formula in this range. 
This is a basic difficulty studied by sieve methods. 

A subtle consequence of this failure is that there must be occasional unexpectedly 
large fluctuations in the number of primes in short intervals away from their expected 
number. This phenomenon was discovered by H. Maier |205j in 1985. It is interesting 
to note that Maier's argument used a "double counting" argument in a tabular array of 
integers that resembles the array summations of Euler and Bernoulli exhibited in section 
12.21 His argument also used properties of the Buchstab function discussed in section [3^ 
See Granville |136| for a detailed discussion of this large fluctuation phenomenon. 

3.4. Euler's constant and arithmetic functions. Euler's constant also appears in 
the behavior of three basic arithmetic functions, the divisor function d{n), Euler's totient 
function <^(n) and the sum of divisors function a{n). For the divisor function d{n) it 
arises in its average behavior, while for (j){n) and cr(n) it concerns extremal asymptotic 
behavior and there appears in the form . These extremal behaviors can be deduced 
starting from Mertens's results. 

A basic arithmetic function is the divisor function 

d{n) := #{d : d\n, l<d<n}. 

which counts the number of divisors of n. In 1849 Dirichlet [72j introduced a method to 
estimate the average size of the divisor function and other arithmetic functions. In the 
statement of Dirichlet's result {x} := x — [x] denotes the fractional part x (mod 1). 
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Theorem 3.4.1. (Dirichlet 1849) (1) The partial sums of the divisor function satisfy 

n 

^d{k) = nlogn + {2-i-l)n + 0{y/n) (3.4.1) 

k=l 

for \ < n < oo, where 7 is Euler's constant. 

(2) The averages of the fractional parts of ^ satisfy 

n 

^{^} = (l-7)n + 0(V^) (3.4.2) 

for 1 < n < 00. 

Proof. Dirichlet 's proof starts from the identity 

n n 

E'^(^) = EL^J' (3.4.3) 

k=i k=i 

which relates the divisor function to integer parts of a scaled harmonic series. The right 
side of this identity is approximated using the harmonic sum estimate 

n 

— = nlogn + 7n + 0(1), 

k=l 

in which Euler's constant appears. The difference term is the sum of fractional parts. 
In Sect. 2 Dirichlet formulates the "hyperbola method" to count the sum of the lattice 
points on the right side of ()3.4.3p . and in Sect. 3 he applies the method to get (1). 
In Sect. 4 he compares the answer ()3.4.ip with ()3.4.3p and obtains the formula ()3.4.2p 
giving (2). A modern treatment of the hyperbola method appears in Tenenbaum [288^ 
1.3.2]. □ 

The problem of obtaining bounds for the size of the remainder term in (I3.4.ip . 

n 

A(n) ■.= "^d{k) - (nlogn+ (27 - l)n), (3.4.4) 

k=l 

is now called the Dirichlet divisor problem. This problem seeks the best possible bound 
for the exponent 6 of the form | A(n)| = 0(n^^'^), valid for any given e > for 1 < n < 00, 
with O-constant depending on e. This problem appears to be very difficult. Dirichlet's 
argument shows that the size of the remainder term A(n) is dominated by the remainder 
term in the fractional part sum (I3.4.2p . In 1916 Hardy [152j established the lower bound 
9 > J. The exponent 6 = 1/4 is conjectured to be the correct answer, and would 
correspond to an expected square-root cancellation. The best fi-result on fluctuations 
of A(n), which does not change the lower bound exponent, is that of Soundararajan 
[2 78 1 in 2003. Concerning upper bounds, a long sequence of improvements in exponent 
bounds has arrived at the current record < ^ 0.31490, obtained by Huxley ( |165j . 
[166]) in 2005. Tsang |297] gives a recent survey on recent work on A(n) and related 
problems. 

Dirichlet noted that (|3.4.2p implies that the fractional parts {{^};1 < k < n} are 
far from being uniformly distributed (modi). In Sect. 4 of [72] he proved that the 
set of fractional parts k having < {j} < ^ has a limiting frequency (2 — log4)n = 
(0.61370. . . )n. and those with i < {^} < 1 have the complementary density (log 4 — 
l)n = (0.38629 . . . )n, as n ^ 00. 
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In 1898 de La Vallee Poussin [184] generalized the fractional part sum estimate (|3.4.2p 
above to summation over k restricted to an arithmetic progression. He showed for each 
integer a > 1 and integer < 6 < a, 

(n-b)/a 
fc=0 

for 1 < n < CO, with 0-constant depending on a. He comments that it is very remarkable 
that these fractional parts approach the same average value 1 — 7 no matter which 
arithmetic progression we restrict them to. Very recently Pillichshammer [241J presented 
further work in this direction. 

We next consider Euler's totient function 

(j)[n) := #{j : 1 < i < n, with gcd(j,n) = 1}, 

which has (/>(ra) = |(Z/nZ)*|. Euler's constant appears in the asymptotics of the minimal 
growth rate (minimal order) of (t){n). In 1903 Edmund Landau [190j showed the following 
result, which he also presented later in his textbook [191' pp. 216-219]. 

Theorem 3.4.2. (Landau 1903) Euler's totient function (j){n) has minimal order jj^p^^- 
Explicitly, 

0(n)loglogn 
iim mf = e ' . 

n— >oo n 

Since (/'(n) = ?^rTp|n(-'-~p)i extremal behavior is attained by the primorial numbers, 
which are the numbers that are products of the first k primes, i.e. A^^ := piP2 ■ ■ ■ Pki for 
example A'4 = 2 • 3 • 5 • 7 = 210. The result then follows using Mertens's product theorem. 

The sum of divisors function is given by 

d\n 

SO that e.g. a{6) = 12. Euler's constant appears in the asymptotics of the maximal 
growth rate (maximal order) of a{n). In 1913 Gronwall [142] obtained the following 
result. 

Theorem 3.4.3. (Gronwall 1913) The sum of divisors function cr{n) has maximal order 
n log log n. Explicitly, 

lim sup = e' . 

n^oo n log log n 

For cr{n) the extremal numbers form a very thin set, with a more complicated descrip- 
tion, which was analyzed by Ramanujan |248j in 1915 (see also Ramanujan |250j) and 
Alaoglu and Erdos [3] in 1944. They have the property that their divisibility by each 
prime has exponent ep(n) — )• 00 as n — )• 00, see §3.7 and formula (|3.7.3p . 

The Euler totient function (j){n) and the sum of divisors function cr(n) are related by 
the inequalities 

—^n^ < 4>(n)a(n) < r? , 

31 



which are easily derived from the identit; 



^{n)a{n) = U (l - ^) 

p<!||n \ -t^ / 

e>l 



30 



p'^Wn means p'^\n and p"^^^ \ n. 
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We note that the ratio of extremal behaviors of the "smallest" (j){n), resp. the "largest" 
(T(n), given in Theorem 13.4.21 and Theorem 13.4.31 have a product asymptotically growing 
like n^. The extremal numbers N^; = piP2 • • - Pk in Landau's theorem are easily shown 
to satisfy (f>{Njf)a{Nf^) ~ '^^k^ ^k ~^ oo. The extremal numbers n in Gronwall's 
theorem can be shown to satisfy the other bound 

4>{n)a{n) ~ 

as n — )• oo. The latter result implies that the Gronwall extremal numbers satisfy 

(/)(n) log log n _ 
lim = e 

n^oo n 

showing they are asymptotically extremal in Landau's theorem as well. 

In §3.71 and §3.81 we discuss more recent results showing that the Riemann hypothesis 
is encoded as an assertion as to how the extremal limits are approached in Landau's 
theorem and Gronwall's theorem. 

3.5. Euler's constant and sieve methods: the Dickman function. Euler's con- 
stant appears in connection with statistics on the sizes of prime factors of general integers. 
Let ^{x,y) count the number of integers n < x having largest prime factor no larger 
than y. These are the numbers remaining if one sieves out all integers divisible by some 
prime exceeding y. The Dirichlet series associated to the set of integers having no prime 
factor > y is the partial zeta function 



p" 

1 



p<y 

and we have 

Diyy 

where D{y) is given in (I3.3.4p . It was shown by Dickman [70] in 1930 that, for each 
u > 0, a positive fraction of all integers n < x have all prime factors less than x^^". 
More precisely, one has the asymptotic formula 

^'(x, x^) ~ p{u)x, 

for a certain function p{u), now called the Dickman function. The notation p{u) was 
introduced by de Bruijn [39], [50].) This function is determined for u > by the prop- 
erties: 

(1) (Initial condition) For < u < 1, it satisfies 

p{u) = 1. 

(2) (Difference-differential equation) For u > 1, 

up{u) = —p{u — 1). 

The Dickman function is alternatively characterized as the solution on the real line to 
the integral equation 

up{u) = p{u- t)dt, (3.5.2) 



JO 

(for up{u)) with initial condition as above, extended to require that p{u) = for all 
u < 0, as shown in de Bruijn [49] (2.1)]. Another explicit form for p{u) is the iterated 
integral form 

' t^kl ■+*fe<" ti t2 tk' ^ ^ 

k=l ti,...,tf.>i 
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which also occurs in connection with random permutations, see ()3.10.2p ff. 

It is known that for u> 1 the Dickman function p{u) is a strictly decreasing positive 
function that decreases rapidly, satisfying 

p{u) < — — — , 
r(n + 1) 

see Norton |235t Lemma 4.7]. The key relation of the Dickman function to Euler's 
constant is that it has total mass e^, see Theorem 13 . 5 . 1 1 b elow . 

In 1951 N. G. de Bruijn |49j gave an exact expression for the Dickman function as a 
contour integral, 

p(^u) = — exp(7+/ ^-^—ds]e'''''dz{u>0). (3.5.4) 
271"^ y_joo \ Jo s J 

One may also consider the one-sided Laplace transform of the Dickman function which 
we denote 

poo 

p{s) := / p{u)e~'^''du, 
Jo 

following Tenenbaum |288l III.5.4]. This integral converges absolutely for all s G C 
and defines p{s) as an entire function. To evaluate it, recall that the complementary 
exponential integral 

Ein(.) := f ^-^dt = Ei-ir-'^z-, (3.5.5) 

n=l 



is an entire function of z. 



Theorem 3.5.1. (de Bruijn 1951, van Lint and Richert 1964) The one-sided Laplace 
transform p{s) of the Dickman function p{u) is an entire function, given by 

p(s) = eT-Ei>^(^), (3.5.6) 

in which Ein(s) is the complementary exponential integral. In particular, the total mass 
of the Dickman function is 



p{0) = / piu)du = e^. (3.5.7) 
Jo 

Proof. Inside de Bruijn's formula (|3.5.4p there appears the integral 

f'' e* - 1 
/(s) := / dt = -Ein(-s). 

Jo t 

This formula determines the Fourier transform of the Dickman function to be e"^^^*-"**^ . 
The change of variable s = —it yields the the one-sided Laplace transform of p{u) above. 
The total mass identity (|3.5.7p was noted in 1964 in van Lint and Richert |199j . Moreover 
they noted that as u — )• oo, 

/ p{t)dt = + O (e"") . 
Jo 

A direct approach to the determination of the Laplace transform given in ()3.5.6p was 
given by Tenenbaum \288\ III. 5. 4, Theorem 7] in 1995. His proof uses a lemma showing 
that for all s G C \ (— oo, 0] one has 

- I{-s) = Ein(s) = 7 + log s + J{s), (3.5.8) 

with 



oo ^—s—t 



J{s) := I -J^dt. (3.5.9) 
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Tenenbaum sketches a direct arithmetic proof of the identity (j3.5.7p in [288, III. 5 Exercise 
2, p. 392]; in this argument the quantity e< is derived using Mertens's product theorem. 
We note the identity 

J(s) = Ei(s), 

in which Ei(s) := J°° ^'^^ is the principal exponential integral. The function Ei(s) 
appeared in connection with Euler's divergent series in §2.51 see ()2.5.14p . □ 

In 1973 Chamayou [54J found a probabilistic interpretation of the Dickman function 
p{x) which is suitable for computing it by Monte Carlo simulation. 

Theorem 3.5.2. (Chamayou 1973) Let Xx^X^^Xz^ ■■■ a sequence of independent iden- 
tically distributed random variables with uniform distribution on [0,1]. Set 

P{u) := Prob[Xi + X1X2 + + • • • < n.] 

This function is well defined for u > and its derivative P'{u) satisfies 

p{u) = e^P'{u). 

A discrete identity relating the Dickman function and e^' was noted by Knuth and 
Trabb-Pardo |175j . which complements the continuous identity ()3.5.7p . 

Theorem 3.5.3. (Knuth and Trabb-Pardo 1976) For < x < 1 the Dickman function 
satisfies the identity 

00 

x + ^{x + n)p{x + n) = e^. (3.5.10) 

n=l 

In particular, taking x = 0, 

p{l) + 2p(2) + 3p(3) + • • • = eT. 

Proof. The identity ()3.5.10|) immediately follows from (I3.5.7P , using the integral equation 
([33:2]) for up{u). □ 

We conclude this section by using formulas above to deduce a curious identity relating 
a special value of this Laplace transform to the Euler-Gompertz constant 5 given in §2.4; 
compare also Theorem 13.16.21 

Theorem 3.5.4. The Laplace transform p{s) of the Dickman function p{u) is given at 
s = 1 by 

/5(1) := / p{u)e~^du = e~'^ . 
Jo 

in which 6 := i-^ogv Euler-Gompertz constant. 

Proof. By definition p{l) = p{u)e~'^du. We start from Hardy's formula (|2.5.1ip for 
the Euler-Gompertz constant: 

5 = s{l)= / -r—dw. 
Jo l + w 

Comparison with ()3.5.9p then yields 



roo -i-t 

Now ()3.5.8p evaluated at s = 1 gives 



e 



Ein(l) =7 + Ei(l) = 7 + -, (3.5.11) 

e 
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which can be rewritten using ()3.5.5p as 



M - e"* 5 

dt = j + -. (3.5.12) 







t ' e 

We conclude from Theorem 13.5.11 that 

p(l) = e^-^-« =e-i, 
which is the assertion. □ 

Hildebrand and Tenenbaum |162| and Granvihe |137j give detailed surveys including 
estimates for a field sometimes called psixyology, see Moree [219] . It has re- 

cently been uncovered that Ramanujan had results on the Dickman function in his 'Lost 
Notebook', prior to the work of Dickman, see Moree [220^ Sect. 2.4]. 

3.6. Euler's constant and sieve methods: the Buchstab function. Euler's con- 
stant also arises in sieve methods in number theory in the complementary problem where 
one removes integers having some small prime factor. Let ^{x,y) count the number of 
integers n < x having no prime factor p < y. In 1937 Buchstab [5T| (see also |52] ) 
established for n > 1 an asymptotic formula 

1 X 

^{x, X") ~ uuj{u)- 



logx 

for a certain function u}{u) named by de Bruijn 08] the Buchstab function. 
The Buchstab function is defined for n > 1 by the properties: 

(1) (Initial conditions) For 1 < u < 2, it satisfies 

Ijj(u) = —. 

u 

(2) (Difference-differential equation) For u > 2, 

{uuj{u)) = uj{u — 1). 

This function is alternatively characterized as the solution on the real line to the integral 
equation 

uj{t)dt, (3.6.1) 

as shown in de Bruijn [48^ (2.1)], who determined properties of this function. Another 
explicit formula for this function, valid for n > 2, is 

1 f dyidy2 ■ ■ ■ dyk-i 



UUJ[U) 



1+ E 



J,, , iA'<^«<i yiy2 ■ ■ ■ 2/fc-i(l - (yi + y2 + h yk~i)) 

It is known that the function oj{u) is oscillatory and satisfies 

lim oj[u) = e~"' . 

The convergence to is extremely rapid, and satisfies the estimate 

uj{u) = e"^ + 0(n-"/2) for u G [2, oo). (3.6.2) 

Proofs of these results can be found in Montgomery and Vaughan [2181 Sec. 7.2] and in 
Tenenbaum [Ml Chap. III.6]. 

The convergence behavior of the Buchstab function to e'*' as x — )• cxd played a crucial 
role in the groundbreaking 1985 work of Helmut Maier |205j on the existence of large 
fluctuations in the density of primes in short intervals that was mentioned at the end 
of §3.31 Maier [2(15] Lemma 4] showed that a;(n) — e~'^ changes sign at least once on 
each interval [a — l,a] for any a > 2. In 1990 Cheer and Goldston [55] showed there 
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are at most two sign changes and at most two critical points on an interval [a — l,a], 
and Hildebrand [l_61j showed that the spacing between such sign changes approaches 1 
as ti — oo, 

A consequence of the analysis of de Bruijn [l8] is that the constant e~'^ appears as a 
universal limiting constant for sieve methods, for a wide range of sieve cutoff functions 
y = y{x) growing neither too slow nor too fast compared to x. 

Theorem 3.6.1. (de Bruijn 1950) Suppose that y = y{x) depends on x in such a way 
that hath y — )• oo and j^f^ — ?■ hold as x ^ oo. Under these conditions we have 

iOg X 

cl>(x,y)~e-^-^. (3.6.3) 
logy 

Proof. This result can be deduced from de Bruijn |481 (1-7)]. This asymptotic formula 
also follows as a corollary of very general estimate ( [2881 III. 6. 2, Theorem 3]) valid 
uniformly on the region x > y > 2, which states 

^ , , , log X , X y ^ / X 

logy logy logy v(logy)^ 
The formula ()3.6.3p follows using the estimate (I3.6.2P for uj{u), which applies since the 
hypotheses x/y — )• oo when x — )• oo imply — oo. □ 

The one-sided Laplace transform of the Buchstab function is given as 

POO 

u;{s) := / e-'"6j(n)dn, (3.6.4) 
Jo 

where we make the convention that uj{u) = 0, for < n < 1. Under this convention, 
this Laplace transform has a simple relation with the Laplace transform of the Dickman 
function. 

Theorem 3.6.2. (Tenenbaum 1995) The one-sided Laplace transform uj{s) defined for 
Re{s) > by ( fg. 6.4^ , extends to a meromorphic function on C, given explicitly by the 
form 

l + u{s) = ^, (s/0). (3.6.5) 
sp{s) 

When s is not real and negative, 

l + ujis) = e-^'-'\ (3.6.6) 

with J{s) = '^glt dt = Ei{s). 

Proof. The Laplace transforms individually have been evaluated many times in the sieve 
method literature, see Wheeler [3101 Theorem 2, ff], who also gives references back to 
the 1960's. One has explicitly that 

w(s) = -exp(-Ein(s)), (3.6.7) 
s 

in terms of the complementary exponential integral. The identity (j3.6.5p connecting 
them seems to be first explicitly stated in Tenenbaum \288\ III. 6, Theorem 5]. □ 

Combining this result with Theorem 13.5.11 gives 

e"^ 

i2}{s) ~ , as s — )• 0. 

s 

We also obtain using ()3.6.6p . (I3.6.7P along with (I2.5.16P that 

/•oo 

cj(l) := / oj{u)e-''du = e"^ exp(-Ei(l)) = e'^e'^/^, (3.6.8) 
Jo 
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where 5 is the Euler-Gompertz constant. 

Wheeler [310\ Theorem 2] also observes that the solutions to the differential-difference 
equations for the Dickman function and the Buchstab function are distinguished from 
those for general initial conditions by the special property that their Laplace transforms 
analytically continue to entire (resp. meromorphic) functions of s G C. 

For a retrospective look at the work of N. G. de Bruijn on p(x) and '^{x,y) in ^3.6^ 
and of u}{x) and ^{x,y) in this section, see Moree |220j . 

3.7. Euler's constant and the Riemann hypothesis. Several formulations of the 
Riemann hypothesis can be given in terms of Euler's constant. Here we describe some 
that involve the approach towards the extremal behaviors of Euler's totient function 
0(n) and the sum of divisors function cr{n) given in ^3.4[ The Riemann hypothesis is 
also related, in another way, to generalized Euler's constants considered in §3.81 

In 1981 J.-L. Nicolas ( [231j . [232j ) proved that the Riemann hypothesis is encoded 
in the property that the Euler totient function values (t){nk) for the extremal numbers 
nk = P1P2 ' ' ' Pk given in Landau's theorem approach their limiting value from one side 
only. Nicolas stated his result using the inverse of Landau's quantity, as follows. 

Theorem 3.7.1. (Nicolas 1981) The Riemann hypothesis holds if and only if all the 
primorial numbers Nk = P1P2 ' ' ' Pk with k >2 satisfy 

Nk ^ 
(j){Nk) loglogNk 

If the Riemann hypothesis is false, then these inequalities will be true for infinitely many 
primorial numbers Nk and false for infinitely many Nk- 

This inequality is equivalent to the statement that the Riemann hypothesis implies 
that all the primorial numbers Nk for k >2 undershoot the asymptotic lower bound e~"' 
in Landau's Theorem 13.4.21 As an example, (j){2 • 3 • 5 • 7) = 0(210) = 48 and 

(^(714) log log n4 (A(210) log log 210 

^ ^ ^ = ' ^ ^ = 0.38321 . . . < e"^ = 0.56145 .... 

n4 210 

For most n one will have ^"g" > e~'^ , and only a very thin subset of n will be close 

to this bound. For example taking n = pk, a single prime, one sees that as A; — )• 00 this 
ratio goes to +00. 

Very recently Nicolas [233j obtained a refined encoding of the Riemann hypothesis in 
terms of a sharper asymptotic for small values of the Euler i;^)-function, in which both 7 
and e'^ appear. 

Theorem 3.7.2. (Nicolas 2012) For each integer n >2 set 

n 



c(n) := ( — e'~' log log n ] \/log n. 
\(p[n) ) 

Then the Riemann hypothesis is equivalent to the statement that 

limsup c(n) = e'^(4 + 7 - log(47r)). (3.7.1) 

Proof. The assertion that the Riemann hypothesis implies (j3.7.ip holds is \233\ Theorem 
1.1], and the converse assertion is [233\ Corollary 1.1]. Here the constant 

(4 + 7 - log(4^)) = e^(2 + /3) = 3.64441 50964..., 

in which 

/3 = y — = 2 + 7 - log(47r) = 0.04619 14179..., 
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where in the sum p runs over the nonreal zeros of the zeta function, counted with 
multiplicity. For the converse direction, he shows that if the Riemann hypothesis fails, 
then 

limsupc(n) = +00. 

n— >oo 

Note that liminf„_^oo c(n) = —00 holds unconditionally. □ 



In 1984 G. Robin [260] showed that for the sum of divisors function (j{n) the Riemann 
hypothesis is also encoded as a one-sided approach to the limit in Gronwall's theorem 
(Theorem [3X3|). 

Theorem 3.7.3. (Robin 1984) The Riemann hypothesis holds if and only if the in- 
equalities 

a(n) 



<e^ (3.7.2) 

n log log n 

are valid for all n > 5041. If the Riemann hypothesis is false, then this inequality will he 
true for infinitely many n and false for infinitely many n. 

Robin's result says that Riemann hypothesis is equivalent to the limiting value e^ 
being approached from below, for all sufficiently large n. The inequality (j3.7.2p fails to 
hold for a few small n, the largest known exception being n = 5040. Here the extremal 
numbers giving record values for f{n) = will have a different form than that in 
Nicolas's theorem. It is known that infinitely many of the extremal numbers will be 
colossally abundant numbers, as defined by Alaoglu and Erdos [3] in 1944. These are 
numbers n such that there is some e > such that 

ain) a(k) ^ 

^ > -i-^, for 1< A; < n. 

Alaoglu and Erdos showed that the "generic" colossally abundant number is a product 
of powers of small primes having a sequence of exponents by a parameter e > as 

n = n(e) := JJp"f(^), 

p 

in which 

, log(p'^'-l)-l°g(p'-l) _ 1 ,3 , 3) 

logp 

For fixed e the sequence of exponents ap(e) is non- increasing as p increases and become 
for large p. However for fixed p the exponent ap(e) increases as e decreases towards 
0, so that n(e) — )• cxd as e — )• 0"*". Colossally abundant and related numbers had actually 
been studied by Ramanujan |248j in 1915, but the relevant part of this paper was sup- 
pressed by the London Mathematical Society to save expense. The suppressed part of 
the manuscript was recovered in his Lost Notebook, and later published in 1997 ( [250]). 

In 2002 this author ( [186] ). starting from Robin's result, obtained the following el- 
ementary criterion for the Riemann hypothesis, involving the sum of divisors function 
cr(n) and the harmonic numbers Hn- 

Theorem 3.7.4. The Riemann hypothesis is equivalent to the assertion that for each 
n > 1 the inequality 

a{n) <e^- log Hn + Hn (3.7.4) 
is valid. Assuming the Riemann hypothesis, equality holds if and only if n = l. 
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The additive term Hn is included in this formula for elegance, as it makes the result 
hold for all n > 1 (assuming RH), rather than be valid for n > 5041, as in Robin's 
theorem 13.7.31 The converse direction of this result requires additional proof, which 
is accomplished in |186) using certain asymptotic estimates obtained in Robin's paper 
[260]. It shows that if the Riemann hypothesis is false then the inequalities ()3.7.4p will 
fail to hold for infinitely many positive integers n. 

3.8. Generalized Euler constants and the Riemann hypothesis. In 1961 W. 

Briggs [461 introduced the notion of an Euler constant associated to an arithmetic pro- 
gression of integers. This notion was studied in detail by D. H. Lehmer |196j in 1975. 
For the arithmetic progression h (mod k) with < h < k we set 

1 log X ^ 



j{h,k) := lim f ^^"j 



0<n<x 
n=h ( mod k) 

These constants were later termed Euler-Lehmer constants by Murty and Saradha |223j . 
Here one has 7(0,1) = 7, 7(1,2) = ^7 and 7(2,4) = ^7, where 7 is Euler's constant. 
It suffices to study the constants where gcd(/i, k) = 1, other cases reduce to these by 
dividing out the greatest common factor. 

We next define generalized Euler constants 7(fi) associated to a finite set of primes 
= {pii, possibly empty. These constants were introduced by Diamond and 

Ford [68j in 2008. We define a zeta function Zq{s) by 



where we set 



in which 



00 

as) = ^ln{n)n-^= ^ n 

n=l gcd{n,Pn)=l 

fl if(n,Pn) = l, 
I otherwise, 



Pn =IIp 

pen 

The function Zq(s) defines a meromorphic function on the entire plane which has a 
simple pole at s = 1 with residue 



^(^) ^= n (1 - ^) • 



pen 

The generalized Euler constant associated to 0, is the constant term in the Laurent 
expansion of Zq(s) around s = 1, namely 

Zn{s) = ^ +7(f^) + f2lkms - I)', 
s — 1 ^-^ 

i=i 

In the case that = is the empty set, this function Zq(s) is exactly the Riemann 
zeta function, and 7(0) = 7 by Theorem 13.2.11 These numbers 7(f^) generalize the 
characterization of Euler's constant in (j3.2.2p in the sense that 

7(1^)= lim y" Mil _ z)(i^) log 
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The numbers 7(^2) are easily shown to be finite sums of Euler-Lehmer constants 

7(0)= Yl ^(h,Pn). (3.8.1) 

l<h<P.p 
gcd(h,P^) = l 

Diamond and Ford [681 Theorem 1] show that these generahzed Euler constants are 
related to Euler's constant in a way which involves the constant e~'^ rather than 7. Let 
f^r denote the set of the first r primes and set '■= 7(^r). 

Theorem 3.8.1. (Diamond and Ford 2008) Let T = inf{7(P) : all finite V}. Then: 

(1) The values of ^{V) are dense in [r,oo). 

(2) The constant F satisfies 

F < e"^. (3.8.2) 
IfT< e~'^ then the value F is attained at some 7,.. Furthemore F > 0.56. 

Proof. (1) This is shown as [68l Prop. 2]. 

(2) Diamond and Ford show that each ^iV) > 7-r for some 1 < r < \V\. This implies 

F = inf{7r : r > 1}. (3.8.3) 

Merten's formulas for sums and products yield the asymptotic formula 

7,. ~ e~^ as r — )• 00. (3.8.4) 

which by (j3.8.3p implies that F < e~'^ . If F < e~'^ then it must be attained by one of the 
7r, For there must be at least one 7^^ < e^"^' , and by (|3.8.4p there are only finitely many 
7r ^ 7ri. and one of these must attain the infimum by the bound above. Diamond and 
Ford established the lower bound ^ 0.56 for all r > 1. □ 

Diamond and Ford [681 Theorem 2] also show that the behavior of the quantities 7(ri) 
is complicated as the set Q, increases, in the sense that 7^ is not a monotone function of 
r, with 7r+i > 7r and 7r+i < 7r each occurring infinitely often. Nevertheless Diamond 
and Ford [684 Theorems 3, 4] obtain the following elegant reformulation of the Riemann 
hypothesis. 

Theorem 3.8.2. (Diamond and Ford 2008) The Riemann hypothesis is equivalent to 
either one of the following assertions. 

(1) The infimum F of ^(V) satisfies 

F = e"^. 

(2) For every finite set VL of primes 

7(ri) > e"^. (3.8.5) 

Here (2) implies (1), and (2) says that if the Riemann hypothesis is true, then (j3.8.5p 
holds for all and the infimum F = is not attained, while if it does not hold, then 
by Theorem 13. 8. II (1) the reverse inequality holds for infinitely many 0, and the infimum 
is attained. Theorem 13.8.21 seems remarkable: on taking "P = to be the empty set, it 
already requires, for the truth of the Riemann hypothesis, that 

7 = 0.57721... > e"^ = 0.56145... 

In consequence the unique real root xq w 0.56714 of the equation xq = e~^° necessarily 
satisfies 7 > j;o > 

Results on the transcendence of Euler-Lehmer constants 7(/i, k) and generalized Euler 
constants 7(f^) have recently been established, see §3.161 
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Finally we note that one may define more generally higher order Euler-Lehmer con- 
stants. For j > 1 Dilcher [71j defines for h> 1, and < k < h, the constants 

0<n<x ' 
n=h{ mod k) 

Here j is the order, with ji{h,k) = j{h,k) being the Euler-Lehmer constants, and 
7j(l,l) = 7j is the j-th Stieltjes constant (see Theorem 13.2. ip . Dilcher relates these 
constants to values of derivatives of the digamma function at rational points, and also 
to derivatives of Dirichlet L-functions evaluated at s = 1. 

3.9. Euler's constant and extreme values of ^(1 + it) and L{l,X-d)- Some of the 
deepest questions in number theory concern the distribution of values of the Riemann 
zeta function and its generalizations, Dirichlet L-functions. In this section we consider 
such value distributions at the special point s = 1 and on the vertical line Re(s) = 1. 
Here Euler's constant appears in connection with the size of extreme values of (^{1 + it) 
as t — 7- oo. It also appears in the same guise in connection with extreme values of the 
Dirichlet L-function L(s, X-d) at the point s = 1 where X-d is the real primitive Dirichlet 
character associated to the quadratic field Q{V—d), and we let d — )• oo. (The quantity —d 
denotes a fundamental discriminant of an imaginary quadratic field, which is necessarily 
squarefree away from the prime 2.) In both case these connections were first made by 
J. E. Littlewood. Some of his results are conditional on the Riemann hypothesis or the 
generalized Riemann hypothesis, while others are unconditional. 

In these results Euler's constant occurs by way of Merten's Theorem 13.3.21 The basic 
idea, in the context of large values of |C(1 + it)\, is that large values will occur (only) 
when a suitable finite truncation of its Euler product representation is large, and this in 
turn will occur when the phases of the individual terms in the product line up properly. 
Mertens's theorem is then used in estimating the size of this Euler product. 

We first consider extreme values of ("(1 + it). To place results in context, it is known 
unconditionally that |C(1 + it)\ = 0{log\t\) and that = 0(log|i|) as \t\ — )■ oo. 

In 1926 Littlewood determined, assuming the Riemann hypothesis, the correct maximal 
order of growth of the Riemann zeta function on the line Re(s) = 1, and in 1928 he 
obtained a corresponding bound for ^^jq:^. 

Theorem 3.9.1. (Littlewood 1926, 1928) Assume the Riemann hypothesis. Then 

(1) The maximal order of \C,{1 + it)\ is at most log log with 

lim sup — — < 2e ' . 

t-^oo log log t 

(2) The maximal order o/ j^^j^^ is at most log log t, with 

^'f:^^ WTW^^t - cM 

Proof. The bound (1) is [2001 Theorem 7] and the bound (2) is [201, Theorem 1]. □. 

Littlewood also established unconditionally a lower bound for the quantity in (1) and, 
conditionally on the Riemann hypothesis, a lower bound for that in (2), both of which 
matched the bounds above up to a factor of 2. In 1949 a method of Chowla [56j made 
this lower bound for the quantity in (2) unconditional. These combined results are as 
follows. 
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Theorem 3.9.2. (Littlewood 1926, Titchmarsh 1933) The following bounds hold uncon- 
ditionally. 

(1) The maximal order of |C(1 + it)\ is at least log log with 

lim sup > e ' . 

t^oo log log t 

(2) The maximal order of jg^^^^^^^j is at least log log t, with 

1 1 „ 

lim sup J— — — \ > 



i^oo"|C(l + »t)|loglogt - C(2) 

Proof. The bound (1) is Theorem 8 of [200J. The lower bound (2) was proved, assuming 
the Riemann hypothesis, in 1928 by Littlewood [201]. An unconditional proof was later 
given by Titchmarsh |294j . see also |295l Theorem 8.9(B)]. □ 

Littlewood's proof of (1) used Diophantine approximation properties of the values 
\ogp for primes p, namely the fact that they are linearly independent over the rationals. 
This guarantees that there exist suitable values of t where the all terms in the Euler 
product np<x(^ ~ prit )~^ have most phases tlogp near (mod 2it) for p < X. 

Littlewood was struck by the fact that the unconditional lower bound for 
limsup(_j.oo ™ Theorem 13.9.2( 1) and the conditional upper bound for it given in 



oo log log t 

Theorem 13.9.1( 1) differ by the simple multiplicative factor 2. He discusses this fact at 
length in [20T|, [202] . The Riemann hypothesis is not sufficient to predict the exact value! 
It seems that more subtle properties of the Riemann zeta function, perhaps related to 
Diophantine approximation properties of prime numbers or of the imaginary parts of 
zeta zeros, play a role in determining the exact constant. Littlewood favored the right 
answer as being the lower bound in Theorem 13.9.21 saving ( |202^ p. 359]): 

The results involving c [= e"'] are evidently final except for a certain 
factor 2. I showed also that on a certain further hypothesis (which there 
is, perhaps, no good reason for believing) this factor 2 disappears. 

In the fullness of time this statement has been elevated to a conjecture attributed to 
Littlewood. 

There is now relatively strong evidence favoring the conjecture that the lower bound 
(1) in Theorem 13.9.21 should be equality. In 2006 Granville and Soundararajan [140] ob- 
tained a bound for the frequency of occurrence of extreme values of |C(l+it)|, considering 
the quantity 

1>T(r) := ^meas{t G [T,2T] : \C{1 + it)\ > eV}, 

where meas denotes Lebesgue measure. They established a result [140t Theorem 1] 
showing that for all sufficiently large T, uniformly in the range 1 <C r < log log T + 20, 
there holds 

$T(r) = exp f - ^^^^ (l + 0{-^ + ir^)'/')) ) . (3.9.1) 



r V \/r logT 

where C is a positive constant. An extension of these ideas established the following 
result ([nni Theorem 2]). 

Theorem 3.9.3. (Granville and Soundararajan 2006) For all sufficiently large T , the 
set of points having 

^max^ |C(1 + it) I > e'^ ( log log T + log log log T - log log log log T - log A + 0(1)) 
has measure at least T^~a^ uniformly for A > 10. 
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They observe that if the estimate (|3.9.ip remained vahd without restriction on the 
range of r, this would correspond to the foUowing stronger assertion. 

Conjecture 3.9.4. (Granville and Soundararajan 2006) There is a constant Ci such that 
for T > 10, 

^max^|C(l + = e'^(loglogr + logloglogT + Ci +o(l)) 

This conjecture implies that the lower bound (1) in Theorem 13.9.21 would be an equality. 

We next consider the distribution of extreme values of the Dirichlet L-functions 
L{s, x~d) at s = 1, for real primitive characters X-diiT-) = (ir)' '^^ere —d is the discrim- 
inant of an imaginary quadratic field Q{\/—d). The Dirichlet L- function 



L{s,x-d) = ^' ^ 

n=l 



n 



n 



is a relative of the Riemann zeta function, having an Euler product and a functional 
equation. Dirichlet's class number formula states that for a fundamental discriminant 
—d with d > 0, 

2t: h{-d) 

L{l,X-d) = r=-, 

w^dyd 

where h{—d) is the order of the ideal class group in ^{y/—d) and W-d is the number of 
units in Q(-v/— d) so that W-d = 2 for d > 4, w-^ = 6 and w_4 = 4. Thus the size of 
-L(l, X-d) encodes information on the size of the class group, and the formula shows that 
L(l,x-d) >0. 

Theorem 3.9.5. (Littlewood 1928) Assume the Generalized Riemann hypothesis for 
all real primitive characters L{s,X-d) where —d is the discriminant of an imaginary 
quadratic field Q{\/—d). Then: 

(1) The maximal order of L{l,x-d) is at most log log with 

L{^,X-d) ^ r, 7 

hmsup- — < 2e'. 

d^oo log log d 

(2) The maximal order of jjj^^-^ is at most log log with 

1 2 
hm sup — — < 



■ L(l,x-d)loglogd - C(2) 

Proof. The bounds (1) and (2) are the content of \202\ Theorem 1]. □ 

There are corresponding unconditional lower bounds that differ from these by a factor 
of 2, as noted by Littlewood in 1928 for the bound (1), and completed by a result of 
Chowla [56j for the bound (2). 

Theorem 3.9.6. (Littlewood 1928, Chowla 1949) The foUowing bounds hold uncondi- 
tionally as d —>• oo with —d a fundamental discriminant. 

(1) The maximal order of L{l,X-d) is at least log log with 

iim sup - — > e ' . 

d^oo log log d 



(2) The maximal order of jjjj^-^ is at least log log d, with 



lim sup — — > — -— 

d^oo ^(1, X-d) log log d C(2) 
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Proof. The bound (1) is due to Littlewood |202j. The bound (2) is due to Chowla 1 56\ 
Theorem 2]. □ 

There is again a factor of 2 difference between the unconditional lower bound and the 
conditional upper bound. In 1999 Montgomery and Vaughan [217j formulated a model 
for the general distribution of sizes of Dirichlet character values at s = 1, and based on 
it they advanced very precise conjectures in favor of the lower bounds being the correct 
answer, as follows. 

Conjecture 3.9.7. (Montgomery and Vaughan 1999) 

(1) For each e > 0, for all D > D{e) there holds 

maxL(l,X-d) < e'^ log log D + (1 + e) log log log D, 

d<D 

(2) For each e > 0, for ah D > D{e) there holds 

S L(^) - Wf + ^((loglogI.)(logloglogI.)) • 

In 2003 Granville and Soundararajan jl39j obtained detailed probabilistic estimates for 
the number of characters having extreme values e'^r. Their results imply unconditionally 
that there are infinitely many d with 

-^(li X-d) > e"^ ( log log d + log log log d — log log log log d—lQ). 

In their later paper |140t Theorem 3] they sketched a proof that for any fixed A > 10 
for all sufficiently large primes q there are at least q^~^ Dirichlet characters x (modg) 
such that 

1-^(1, x)\ > e'^(loglogd + log log log d — log log log log d — log^ + 0(1)) . 

Finally we remark on a related problem, concerning the distribution of the phase 
arg((^(l + it)). In 1972 Pavlov and Faddeev j238j observed that the phase of C(l + it) 
appears in the scattering matrix data for the hyperbolic Laplacian operator acting on 
the modular surface ^(1) = -P5'L(2, Z)\E[, where M. = {z = x + iy : y = Im{z) > 0} 
denotes the upper half plane. Then in 1980 Lax and Phillips [194] Sect. Theorem 7.19] 
treated this case as an example of their version of scattering theory for automorphic 
functions (given in |193j ). Recently Y. Lamzouri |189j obtained interesting results on 
the joint distribution of the modulus and phase (|C(1 + arg((^(l + it))). 

3.10. Euler's constant and random permutations: cycle structure. Let Sn de- 
note the symmetric group of all permutations on [1, A^] := {1, 2, 3, • • • , A''}. By a random 
permutation we mean an element a G Sm drawn with the uniform distribution, picked 
with probability We view o" as a product of distinct cycles, and let Cj{a) count the 
number of cycles of length j in a. Three interesting statistics on the cycle structure of 
a permutation a G are its total number of cycles 

n 

i=i 

the length of its longest cycle 

M{a) := max{j : Cj > 0}, 
and the length of its shortest cycle 



m{a) := min{j : Cj > 0}. 
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Euler's constant 7 appears in connection with the the distributions of each of these 
statistics, viewed as random variables on Sn- 

In 1939 J. Touchard [296^ p. 247] expressed the distribution of cycle lengths using 
exponential generating function methods. In 1944 Goncharov [1321 Trans, pp. 31- 
34] (announced 1942 |131j ) also gave generating functions, derived in a probabilistic 
context. He showed that if one lets c{N, k) denote the number of permutations in Sn 
having exactly k cycles, then one has the generating function 

N 

c{N, k)x^ = x{x + l)(x + 2) • • • (x + - 1). (3.10.1) 

fc=i 

Goncharov computed the mean and variance of the number of cycles n{a), as follows. 

Theorem 3.10.1. (Touchard 1939, Goncharov 1944) Draw a random permutation a 
from the symmetric group Sn on N elements with the uniform distribution. 
(1) The expected value of the number of cycles in a is: 

^ 1 

i=i 

In particular one has the estimate 

i?[n(a)] = log iV + 7 + 0(1). 



(2) The variance of the number of cycles in a is: 

Var\n{a)\ := E[{n{a) - E[n{a)]f] = Hm - Hm,2 = Y,-. - Y.-^- 
In particular one has the estimate 



N , N 



7 — T 

j=i ■' i=i ■' 



^/Var[n{a)] = + (i " J^) ^f=^ + ^ii^og iV)"!). 

Proof In 1939 Touchard [Ml P- 291] obtained the formula E[n{a)] = Hn- In 1944 
Goncharov [132, Sect. 15] derived both (1) and (2), using ()3.10.ip . A similar proof was 
given in 1953 by Greenwood [141^ p. 404]. Note also that for all J > 1 the J-th moment 
of n(a) is given by 

E[n{ay] = ^{x^yMr,{x)\.=,. 

It follows using (j3.10.ip that the J-th moment can be expressed as a polynomial in the 
m-harmonic numbers H^^m for 1 < m < J. 

Goncharov |132l Sect. 16] used the moments to show furthermore that the normalized 
random variables nia) := "(°")~-^["(°")] satisfy a central limit theorem as — t- cxd, i.e. 

they converge in distribution to standard unit normal distribution. □ 

Goncharov [132| also obtained information on the distribution of the maximum cycle 

M((t) as well. Define the scaled random variable := ■^M((t). Then, as — t- 

00, Goncharov showed these random variables have a limiting distribution, with the 
convergence in distribution — t-^ Li, in which Li is the probability distribution 
supported on [0, 1] whose cumulative distribution function 

Fi(q) := Prob[0 < Li < a] 
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is given for < q < 1 by 

fc=l ti,-- ,tf.>i 

with Fi{a) = 1 for a > 1. Much later Knuth and Trabb-Pardo [175^ Sec. 10] in 1976 
showed that this distribution is connected to the Dickman function p{u) in ^3.51 bv 

Fi(q) = p(-), < a < 1, (3.10.3) 

a 

and this yields the previously given formula ()3.5.3p . From this connection we deduce 
that this distribution has a continuous density fi{a) da given by 

= -p'{-)\ = p{- — -)-, 
a a a 

with the last equality deduced using the difference-differential equation (??) for p{u). 

In 1966 Shepp and Lloyd |263] obtained detailed information on the distribution of 
both m{a) and M{a) and also of the r-th longest and r-th shortest cycles. Euler's 
constant appears in the distributions of longest and shortest cycles, as follows. 

Theorem 3.10.2. (Shepp and Lloyd 1966) Pick a random permutation a on N elements 
with the uniform distribution. 

(1) Let Mf.(cr) denote the length of the r-th longest cycle under iteration of the per- 
mutation a € S^. Then the k-th moment E[Mr{a)^] satisfies, for k > 1 and r > 1, as 
N ^ oo, 

for positive limiting constants Gr^k- These constants are explicitly given by 

Gr,, = [ "-^^^^^--'-'-^'^-'dx, (3.10.4) 

in which Ei{—x) = ^T^^ exponential integral. They satisfy, for fixed k, as 

r — )• oo, 

lini{k + lYGr,k = ^e-''\ 

r^oo /e! 

in which 7 is Euler's constant. 

(2) Let rriricr) denote the length of the r-th shortest cycle of the permutation a G S'jy. 
Then the expected value E[mr{cr)] satisfies, as N ^ 00, 

lim 7^/^ = -e-^. 3.10.5 
v^oo (logA^)'' r! ^ ^ 

Furthermore the k-th moment E[mr{a)'^] for k > 2 and r > 1 satisfies, as N ^ 00, 

lim /i^-^-^^^J , = -J— r -^e^^i-^^e-^dx. (3.10.6) 
AT^oo iVfc-i(logiV)'-i (r-l)!yo (^-1)! 

Proof. The exponential integral Ei{x), defined for x < is given in (|2.5.13|) . Here (1) 
appears in Shepp and Lloyd |263[ eqns. (13) and (14) ff., p. 347]. Here (2) appears as 
[2S1 eqn. (22), p. 352]. □ 

Shepp and Lloyd {263j also obtained the limiting distribution of the r-th longest cy- 
cle as — )• 00. In terms of the scaled variables Li^^ := j^Mr{a), they deduced the 
convergence in distribution 

4^) L„ 
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in which is a distribution supported on [0, ^] whose cumulative distribution function 
Fr(a) := Prob[0 < Li < q], < a < J, 



is given by 

Fria) = l + ^ (;~l)!(^37)!i / 



111 

' I- cUidh dtn 

1+ ■■+*fc<sr ti t2 tn '' 

«i.---,tfc>i 



and we set Fr{a) = 1 for ^ < a < oo. This extends to general r the case r = 1 treated 
by Goncharov |132j . In 1976 Knuth and Trabb-Pardo [1751 Sec. 4] defined Pr{u) by 
Pr{u) := Fr{^), so that pi{u) = p{u) by (jS.lO.Sp . The functions pr{u) for r > 1 are 
uniquely determined by the conditions 

(1) (Initial condition) For < u < 1, it satisfies 

Pr{u) = 1. 

(2) (Difference-differential equation) For u > 1, 

Upl{u) = -pr{u - 1) + Pr-l{u - 1), 

with the convention po{u) = 0. 
Knuth and Trabb-Pardo formulated these conditions in the integral equation form 

P dt 

Pr{u) = l-J {pr{t - 1) - pr^i{t - 1)) — . 

One may directly check that these equations imply pr{u) = 1 for < n < r. 
Associated to the longest cycle distribution is another interesting constant 

The formula (j3.10.4p gives 

POO 

X= e-^'-^^^-^^dx = 0.62432 99885 .. . (3.10.7) 
Jo 

This constant A is now named the Golomb-Dickman constant in Finch |1131 Sec. 5.4], 
for the following reasons. It was first encountered in 1959 work of Golomb, Welch and 
Goldstein [126], in a study of the asymptotics of statistics associated to shift register 
sequences. Their paper expressed this constant by an integral involving a solution to a 
differential-difference equation, described later in Golomb |128[ p. 91, equation (33)]. In 
1964 Golomb [127| asked for a closed form for this constant, and the Shepp and Lloyd's 
work answered it with (I3.10.7p above. It was evaluated to 53 places by Mitchell |216j in 
1968. In 1976 Knuth and Trabb-Pardo Sect. 10] showed that 

A = 1 - r ^du, (3.10.8) 
Ji W 

where p{x) is the Dickman function discussed in ^3.51 This associates the name of 
Dickman with this constant. In fact the constant already appears in de Bruijn's 1951 
work [501 (5-2)], where he showed 

A = / , ^^^}.r, du, 
Jo + 1)2 
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as noted Wheeler [3101 p- 516]. According to Golomb (128^ p. 192] another integral 
formula for this constant, found by Nathan Fine, is 

A= /\"(^)dx, (3.10.9) 



in which Li{x) := for < x < 1. In 1990 Wheeler pp. 516-517] established 

the formula 

Jo u + 2 

empirically noted earlier by Knuth and Trabb-Pardo, and also the formula 

/•oo 



in which Ein(t) is given in (I3.5.5p . The Golomb-Dickman constant A is not known to 
be related to either Euler's constant 7 or the Euler-Gompertz constant 5, and it is not 
known whether it is irrational or transcendental. 

The Golomb-Dickman constant A appears as a basic statistic in the limit distribution 
Li for the scaled longest cycle, whose expected value is given by 

rl /"^ 1 / 1 

E[]Li]= adFi{a) = - -p'(-)da. 
Jo Jo a a 

Letting x = ^ we obtain 

E\u] =- r ^dx = r p{x - = A, 

Jl X Jl 

see [175, (9.2)]. 

In 1996 Gourdon [135^ Chap. VII, Theoreme 2] determined a complete asymptotic 
expansion of the expected values E[M{a)] as — )■ 00 in powers of which both the 
Golomb-Dickman constant and Euler's constant appear, with the initial terms being 



£|MMl = AiV+lA-£:i + o(-i,) 



The omitted higher order terms for /c > 2 in Gourdon's asymptotic expansion are of 
form "^^f^^ in which each Pk{N) is a oscillatory function of A^ which is periodic with an 
integer period pk that grows with /c as A; — ?■ 00. 

A result in a different direction relates Euler's constant to the probability of all cycles 
having distinct lengths as A^ — t- 00. This result follows from work of D. H. Lehmer |195j . 
Theorem 3.10.3. (Lehmer 1972) Draw a random permutation a on N elements with 
the uniform distribution. Let Pdi^) denotes the probability that the cycles of a have 
distinct lengths, then 

lim Pd{N) = e-^. 

Proof. The cycle lengths of a permutation a form a partition A = (1^^^, 2*^^, • • • , A^^'^) of 
A^ ( written A h A^), having n(A) = ci + C2 + • • • + cat cycles with cycle length j occurring 
Cj times. The number of permutations A^(A) giving rise to a given partition A is 

A^! 

^^^^ ~ ci!c2!---CAr!l^i2'=2 ••• A^cjv' 

and by definition we have 



XhN 
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In the case that ah cycles are of distinct lengths, then all Cj = or 1, and the formula 
for A^(A) becomes 

W 

N{\) = , 

aia2 ■ ■ ■ ttn 

in which ai < 02 < • • • < a„ with n = n(A) are the lengths of the cycles. (The division by 
Oj reflects the fact that a cyclic shift of a cycle is the same cycle.) Lehmer [195] assigned 
weights (l'^i2'^2 • • • N'^'^)~^ to each partition, and counted different sets of such weighted 
partitions. Theorem 2 counted partitions with distinct summands, in his notation 



X — ^ ' 



1 



X\-N ^ ^ " 

with the prime indicating the sum runs over partitions having distinct parts, and showed 
that — )• e~'^ as — )• 00. (Lehmer 's weights agree with the weights assigned to 
partitions from the uniform distribution on S'at only for those partitions having distinct 
parts.) □ 

3.11. Euler's constant and random permutations: shortest cycle. There are 
striking parallels between the distribution of cycles of random permutations of Sat and 
the distributions of factorizations of random numbers in an interval [l,n], particularly 
concerning the number of these having factorizations of restricted types, allowing either 
factorizations with only small factors, treated in §3.5( or having only large prime factors, 
treated in §3.61 where we take N to be on the order of logn. This relation was first 
observed by Knuth and Trabb-Pardo |175j . whose 1976 work (already mentioned in §3.51 
and §3.10p was done to analyze the performance of an algorithm for factoring integers 
by trial division. They observed in |175l Sect. 10, p. 344]: 

Therefore, if we are factoring the digits of a random m-digit number, the 
distribution of the number of digits in its prime factors is approximately 
the same as the distribution of the cycle lengths in a random permutation 
on m elements! (Note that there are approximately Inm factors, and 
Inm cycles.) 

They uncovered this connection after noticing that the Shepp-Lloyd formula for Golomb's 
constant A in (I3.10.7p . as evaluated by Mitchell [216], agreed to 10 places with their 
calculation of Dickman's constant, as defined by (I3.10.8p . A comparison of the k-th 
longest cycle formulas of Shepp and Lloyd in Theorem 13.10.21 for cycles of length aN 
then revealed matching relatives pk{x) of the Dickman function (defined in §3.5p . These 
parallels extend to the appearance of the Buchstab function in the distribution of random 
permutations having no short cycle, as we describe below. 

Concerning the shortest cycle, recall first that a special case of Theorem 13.10.2( 2) 
shows that the expected length of the shortest cycle is 

E[m{a)] = e'^ log iV(l + o(l)) as TV ^ 00, 

It is possible to get exact combinatorial formulas for the probability 

P{N,m) := Prob[m(o") > m : a £ Sn]- 

that a permutation has no cycle of length shorter than m. There are several limiting 
behaviors of this probability as — t- 00, depending on the way m = m{N) grows as 
N ^ 00. There are three regimes, first, where m is constant, secondly where m{N) = aN 
grows proportionally to A^, with < a < 1 (i.e. only long cycles occur) and thirdly, the 
intermediate regime, where m{N) — )• 00 but "^^^ — >• as — t- 00. 
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First, suppose that m{N) = m is constant. This case has a long history. The derange- 
ment problem (Probleme des menages), concerns the probabihty that a permutation has 
no fixed point, which is the case m = 2. It was raised by Remond de Montmort |251j in 
1708 and solved by Nicholas Bernoulli in 1713, see [2521 pp. 301—303]. The derangement 
problem was also solved by Euler in 1753 [99', E201], who was unaware of the earlier 
work. The well known answer is that 



P(iV,2) 



and this exact formula yields the result 



TV 

E 

j=0 



i-iy 



lim P{N,2) = -. 
Af-5-oo e 



In 1952 Cruder |144j gave a generalization for permutations having no cycle shorter than 
a fixed m, as follows. 

Theorem 3.11.1. (Cruder 1952) For fixed m > 1 there holds 



lim Prob[?Ti(cr) > m : a G Sn] = e 
in which Hm denotes the m-th harmonic number. 



-H„ 



(3.11.1) 



Proof. Let PN{m,k) count the number of permutations of having exactly k cycles, 

each of length at least m, and set PN{m) = Yl,k=i PN{fn^ k), so that P{N, m) = -^-^^^r^. 
These are given in the exponential generating function 



oo TV 
Z 



N=0 



N 



k=l 



exp 



exp 



u(log 



1 



m—1 

E 



1 



Cruder [1441 Sect. 8, eqn. (85)] derives from the case u = 1, 

exp 



m—1 



m—1 



N=0 



m 



1 — X 



Using this fact Cruder [1441 Sect. 9, eqn. (94)] derives limAr^.oo p^^m) 
the result. 



□ 



Cruder also observes that one obtains Euler's constant from these values via the scaling 
limit 

lim lim mP{N,m) = e~^. 

m— >oo N^oo 

An asymptotic expansion concerning convergence to the limit (j3.11.ip for constant m 
is now available in general circumstances. Panario and Richmond [2371 Theorem 3, eq. 
(9)], give such an asymptotic expansion for the r-th largest cycle being > m for a large 
class of probability models. 

Secondly, for the intermediate range where m = m{N) — )• oo with ^ — t- as — t- oo, 
there is the following striking universal limit estimate, in which e~'^ appears. 



58 



JEFFREY C. LAG ARIAS 



Theorem 3.11.2. (Panario and Richmond 2001) Consider permutations on N letters 
whose shortest cycle m{a) > m. Suppose that m = m{N) depends on N in such a way 
that m{N) — )• oo and "^^^ — )• both hold as N ^ oo. Under these conditions we have 

c — 

Fioh[m{a) > m{N) : a e Sn] ^ , as TV oo. 

m 

Proof. This resuh is deducible form a general asymptotic result of Panario and Richmond 
[2371 Theorem 3, eq. (11)], which apphes to distribution of the r-th largest cycle of a 
logarithmic combinatorial structure. □ 

Thirdly, we consider the large range, concerning those permutations having shortest 
cycle of length at least a constant fraction aN of the size of the permutation. Here the 
resulting scaled density involves the Buchstab function uj{u) treated in ^3.6| with a = i, 
with u > 1. 

Theorem 3.11.3. (Panario and Richmond 2001) Consider permutations on N letters 
whose shortest cycle m{a) > m. Suppose that m = m{N) depends on N in such a way 
that "j^^ — )• a holds as N ^ oo, with < a < 1. Under these conditions we have, for 
fixed < a < 1, that 

Fioh[m{a) > aN : a e Sn] ^ oj{l/a)^, as N oo, (3.11.2) 

where uj{u) denotes the Buchstab function. 

Proof. This result follows from a general result of Panario and Richmond [2371 Theorem 
3, eq. (10)], which gives an asymptotic expansion for the r-th largest cycle. The proof 
for this case uses in part methods of Flajolet and Odlyzko [117] . □ 

The matching of the large range estimate ()3.11.2p as q — t- with the intermediate 
range estimate follows from the limiting behavior of the Buchstab function uj{u) — )• e~'^ as 
ti — 7- oo. One may compare Theorem l3. 1 1 . 2l with the sieving result in Theorem l3.6.11 They 
are quite parallel, and in both cases the factor e~'^ arises from the limiting asymptotic 
behavior of the Buchstab function. 

The appearance of the Buchstab function in the sieving models in ^3.61 and the per- 
mutation cycle structure above is precisely accounted for in stochastic models developed 
by Arratia, Barbour, and Tavare [17] . [18| in the late 1990's, which they termed Poisson- 
Dirichlet processes (and which we do not define here). They showed ([17], [19', Sec. 1.1, 
1.2] that (logarithmically scaled) factorizations of random integers, and (scaled) cycle 
structures of random permutations as — t- oo lead to identical limiting processes. To 
describe the scalings, one considers an ordered prime factorization of an random integer 
m = Y\i=iPij with 2 < pi < p2 < ■ ■ ■ < Pk- Let k := il.{m) count the number of prime 
factors of m with multiplicity, one defines the logarithmically scaled quantities 

Oj := , 1 <t < k. 

log m 

The random quantities drawn are (0^,02, a|.), where k itself is also a random variable, 
and these quantities sum to 1 . To a random permutation a of Sjy let ai < 02 < • • • < 
a„(o-) represent the lengths of the cycles arranged in increasing order. Then associate to 
a the normalized cycle lengths 

fli := l<i<n. 

Here n = n{a) is a random variable, and these quantities also sum to 1. Arratia, 
Barbour and Tavare [IB] showed that in both cases, as A^ — )• 00 (resp. n — )• 00) these 
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random variables converge to a limiting stochastic process, a Poisson-Dirichlet process 
of parameter 1. They furthermore noted ( [19f p. 34]) that both models required exactly 
the same normalizing constant: e~'^. 

The model for random cycles applies as well to random factorization of polynomials 
of degree N over the finite field GF{q), as g — )• oo, see Arratia, Barbour and Tavare |16] . 



3.12. Euler's constant and random finite functions. Euler's constant also appears 
in connection with the distribution of cycles in a random finite function F : [1,N] — )• 
[1 , A*"] . Iteration of a finite function on a given initial seed leads to iterates of an element 
having a preperiodic part, followed by arrival at a cycle of the function. Thus not every 
element of N belongs to a cycle. 

In 1968 Purdom and Williams [246] established results for random finite functions that 
are analogous to results for random permutations given in §3.101 They studied the length 
of the longest and shortest cycles for a random function; note however that the expected 
length of the longest cycle grows proportionally to ^/N, rather than proportional to 
as in the case of a random permutation. 

Theorem 3.12.1. (Purdom and Williams 1968) Pick a random function F : [1, A^] — )• 

[1, A^] with the uniform distribution over all such functions. 

(1) The length M{F) of the longest cycle of F under iteration satisfies 

,i„ mm = 

N^oo ^/N V 2 

in which A is the Golomb-Dickman constant. 

(2) Let m{F) denote the length of the shortest cycle under iteration of the function 
F . Then the expected value E[m{F)\ satisfies 

lim S^ffll = i,-.. 

Af^oo log A^ 2 

Proof. For fixed A^, Purdom and Wilhams |246| use the notation E[M{F)] = Ep^^iii) := 
E[M{F)] and Ef^^i{s) := E[m{F)]. The result (1) follows from [MQ p. 550, second 
equation from bottom]. Purdom and Williams state a result related to (2) without 
detailed proof. Namely, the top equation on |246(, p. 551] states 



E[miF)] = 5i,iQ„(l, 1) + o(Q„(l, 1)), 



where 



j=i ^ -^^ 

Here the values Sr^k are taken from Shepp and Lloyd |263j . where for k = 1 they are 
given by the right hand side of (I3.10.5p . and for /c > 2 by the right hand side of (I3.10.6p . 
and in particular ^i^i = e~'^. To estimate Qn(l, 1) given by (I3.12.ip we use the identity 
( [2161 p. 550]) 

«»(»)= ^:^^o=l■ 

We observe that Qni)--, 1) is the expected value of the function logj with respect to the 
probability distribution given by the individual terms in (5n(0). One may deduce that 
(5n(l, 1) = ^ logn + o(logn) by showing that this probability distribution is concentrated 
on values j = ^/n + o[^/n). This gives (2) . □ 
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The study of random finite functions is relevant for understanding the behavior of the 
Pollard [242j "rho" method of factoring large integers. The paper of Knuth and Trabb- 
Pardo |175j analyzed the performance of this algorithm, but their model used random 
permutations rather than random functions. An extensive study of limit theorems for 
statistics for random finite functions was made by V. F. Kolchin [177j . A very nice 
asymptotic analysis of the major statistics for random finite functions is given in Flajolet 
and Odlyzko Theorems 2 and 3]. 

3.13. Euler's constant as a Lyapunov exponent. Certain properties of Euler's con- 
stant have formal resemblance to properties of a dynamical entropy. More precisely, the 
quantity C = might be interpretable as the growth rate of some dynamical quantity, 
with 7 = log C being the associated (topological or metric) entropy. Here we present re- 
sults about the growth rates of the size of products of random matrices having normally 
distributed entires. These growth rates naturally involve 7, together with other known 
constants. 

Before stating results precisely, we remark that random matrix products often have 
well-defined growth rates, which can be stated either in terms of growth rates of matrix 
norms or in terms of Lyapunov exponents. Let {A{j) : i > 1} be a fixed sequence of 
dx d complex matrices., and let M{n) = A{n)A{n — 1) • • • ^(2)74(1) be product of dx d 
complex matrices. We define for an initial vector v E and a vector norm || • ||, the 
logarithmic growth exponent of v to be 



A(v):=limsupilogM!M. (3.13.I) 

n— ^00 n V 



n—>oo n 

This value depends on v and on the particular sequence {A{j) : j > 1} but is independent 
of the choice of the vector norm. For a fixed infinite sequence {A(j) : j > 1} as v varies 
there will at most d distinct values A* < A2 < • • • < for this limit. 

Now let the sequence {A(j) : j > 1} vary, with each A{j) being drawn independently 
from a given probability distribution on the set oldxd matrices. In 1969 Furstenberg and 
Kesten |119j showed that for almost all vectors v (off a set of positive codimension) the 
limiting value A^ takes with probability one a constant value A^, provided the distribution 
satisfied the condition 

/ log+||^(x)||d/x(j;) < 00, (3.13.2) 

where (X, $],;u) is the probability space, || • || is a fixed matrix norm, required to be 
submultiplicative, i..e IIM1M2II < ||Mi||||M2||, and where log"*" \x\ := max(0, log |x|). A 
corollary of the multiplicative ergodic theorem proved by Oseledec in 1968 (see [236], 
|247j ) asserts that with probability one all these logarithmic exponents take constant 
values, provided the probability distribution on random matrices satisfies ()3.13.2p . These 
constant values Ai < A2 < • • • < A^ are called the Lyapunov exponents of the random 
matrix product, with A^ being the maximal Lyapunov exponent. For almost all random 
products and almost all starting vectors vq, one will have A(v) = A^, compare Pollicott 
[243| . The value := exp(Ad) may be thought of as approximating the exponential 
growth rate of matrix norms ||M(n)"^M(n)||-'^/^, as n — t- c«. 

In general the Lyapunov exponents are hard to determine. In 1984 J. Cohen and C. M. 
Newman [601 Theorem 2.5] obtained the following explicit result for normally distributed 
variables. Note that condition (j3.13.2p holds for such variables. 

Theorem 3.13.1. (Cohen and Newman 1984) Let{A{l)ij : 1 < i,j < d} be independent, 
identically distributed (i.i.d.) normal random variables with mean and variance s, i.e. 



EULER'S CONSTANT: EULER'S WORK AND MODERN DEVELOPMENTS 61 

A{l)ij ~ N{0,s'^). Let v(0) he a nonzero vector in W^. Define for independent random 
draws 

v(n) = A{n)A{n -!)••• A(l)v(0), n = 1, 2, 3, ... 

and consideW\ 

A(v(0)) := lim -log||v(n)|| (3.13.3) 

n— >oo ri 

using the Euclidean norm || • || on W^. Then with probability one the limit exists in 
\3.13.3\) . is nonrandom, is independent ofv{0), and is given by 

A = ^[log(s2) + log2 + V(^)], (3.13.4) 

in which ijj{x) is the digamma function. Moreover as n —)• oo the random variables 
log(e~''*"||v(n)||) converge in distribution to N{0,a'^) with 

= i^'(-). 

4^ ^2^ 

It is easy to deduce from Cohen and Newman's result the fohowing consequence con- 
cerning d X d (non-commutative) random matrix products. 

Theorem 3.13.2. Let A{j) be dx d matrices with entries drawn as independent identi- 
cally distributed (i.i.d.) normal random variables with variance 1, i.e. A^(0, 1). Consider 
the random matrix product 

S{n) = A{n)A{n-l)--- A{1), 



and let \\S\\p = y^^^j \Sij\'^ be the Frobenius matrix norm. Then, with probability one, 
the growth rate is given by 

{2e-'r+^d/2-i if d is even, 

(3.13.5) 
l^-,+2H,_,-H,,.,y, ,fd^sodd. 

Here Hn denotes the n-th harmonic number, using the convention that Hq = 0. 

Proof. We use an exponentiated version of the formula (j3.13.4p in Theorem l3.13.1l setting 
the variance s = 1 and multiplying both sides of (j3.13.4|) by 2 before exponentiating. We 
obtain, for a nonzero vector v(0), that with probability one, 

lim (||5(n)v(0)|n" =2e^(i). (3.13.6) 

n— >oo 

where || • || is the Euclidean norm on vectors. To obtain the Frobenius norm bound we 
use the identity ||5'|||. = X]j=i I I'S'^jl P where e^ is the standard orthonormal basis of 
column vectors. On taking vectors Vj(0) = ej we obtain 

d 

\\Sin)\\l = ^\\S{n)^rM\\^ 

and using (I3.13.6p and the fact that limn-s.00 = 1, we obtain with probability one that 

lim (||5(n)|||)" =2e^(i). 



■^"'^Cohen and Newman write the left sides of (|3.13.3() and (|3.13.4p as log A, but here we replace these 
by the rescaled variables A(vo) and A, in order to match the definition of Lyapunov exponent given in 
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On applying the digamma formulas for ^(^) given in Theorem l3.1.H with the convention 
Hq = 0, and noting for d = 2m + 1 that 

2 2 2 ' 

we obtain (13.13.50 . □ 

Rephrased in terms of Lyapunov exponents, Theorem 13 . 13.21 savs that the correspond- 
ing d X d matrix system, for odd dimension d has top Lyapunov exponent 

Xd = l (-7 - log 2 + 2Fd„i - , 

and for even dimension d has top Lyapunov exponent 

Ad = ^ (-7 + log2 + F|„^) . 

The simplest expressions occur in dimension d = 1, giving Ai = ^(—7 — log 2), and in 
d = 2, giving A2 = ^(-7 + log 2). 

Note that any real number can be obtained as a maximal Lyapunov exponent in 
Theorem 13.13.11 simply by adjusting the variance parameter s properly in the normal 
distribution A^(0, s). It follows that the content of the specific formulas in Theorem l3. 13.21 
is in part arithmetic, where one specifies the variance s of the normal distributions to be 
a fixed rational number. Moreover, if we rescale s to choose normal distributions A^(0, s) 
with variance s = 2 in dimension 1 and with variance s = ^ in dimension 2, then we 
obtain the Lyapunov exponents A^^ = ^ (~7) these dimensions. 

Returning to the theme of Euler's constant being (possibly) interpretable as a dynam- 
ical entropy, one approach to the Riemann hypothesis suggests that it may be related 
to an (as yet unknown) arithmetical dynamical system, cf. Deninger |64j . [65], [66], |67j . 
There is a statistical mechanics interpretation of the Riemann zeta function as a par- 
tition function given in Bost and Connes [30], with subsequent developments described 
in Connes and MarcoUi [6IJ. For some other views on arithmetical dynamical systems 
related to zeta functions see Lagarias |185j . |187j and Lapidus fl92', Chap. 5]. Here 
we note the coincidence that in formulations of the Riemann hypothesis given earlier in 
Theorems 13.7.11 and 13.7.31 a "growth rate" e^'^ naturally appears. 

3.14. Euler's constant and periods. In 2001 M. Kontsevich and D. Zagier [181] de- 
fined a period to be a complex number whose real and imaginary parts are values of ab- 
solutely convergent integrals of rational functions with rational coefficents, over domains 
in given by polynomial inequalities with rational coefficients. More conceptually 
these are values of integrals of an algebraic differential form with algebraic coefficients, 
integrated over an algebraic cycle, see Kontsevich jlSOi Sec. 4.3]. We will call such 
an integral a period integral. They observed that for the real part one can reduce such 
integrals to the case where all coefficients are rational numbers, all differentials are top 
degree, integrals are taken over divisors with normal crossings defined over Q. In partic- 
ular they introduced a ring V of effective periods, which includes the field of all algebraic 
numbers Q as a subring, as well as a ring of extended periods V obtained by adjoining 
to V the constant which conjecturally is not a period. 
One has, for n > 2 the representation 
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as a period integral. Thus convergent positive integer zeta values are periods. One more 
generally can represent all (convergent) multiple zeta values 

C(ni,n2,...,nfc) := V m nl nj:- 

mi>m2>...>mfe>0 i ^ /« 

at positive integer arguments as period integrals. Other period integrals are associated 
with special values of L-functions, a direction formulated by Deligne [62]. An overview 
paper of Zagier j314j gives a tour of many places that multiple zeta values appear. 

Some other examples of periods given by Kontsevich and Zagier are the values log a 
for Q a positive real algebraic number. In addition the values r(2)'3 where | is a positive 
rational number are periods, see Kontsevich and Zagier [18H p. 775], Andre Chap. 
24]. This can be deduced using Euler's Beta integral ()2.3.ip . Another well known 
example is (pQl (5.11)]) 

dt (r(i))2 



giving a period of an elliptic curve with complex multiplication. 

Periods arise in many places. They arise in evaluating Feynman integrals in quantum 
field theory calculations. For more work on periods and zeta values in this context, see 
Belkale and Brosnan [28]. Periods appear as multiple zeta values in Drinfeld's associator, 
cf. Drinfeld [73], [73], [75]. Periods appear in Vassiliev's knot invariants, due to the 
connection found by Kontsevich [179j , see Bar-Natan [25] . They appear in expansions of 
Selberg integrals, cf. Terasoma [291] . They appear in the theory of motives as periods of 
mixed Tate motives, cf. Terasoma [292], [293] . and Deligne and Goncharov [63]. Other 
constants that are periods are values of polylogarithms at integers. Algebraic multiples of 
beta values B{^, f^) occur as periods of differentials on abelian varieties, cf. Gross |143[ 
Rohrlich appendix]. In 2002 P. Cartier (53) gave a Bourbaki seminar expose surveying 
multiple zeta values, polylogarithms and related quantities. 

There is a good deal known about the irrationality or transcendence of specific periods. 
The transcendence of C(2ra) for n > 1 is immediate from their expression as powers of tt. 
For odd zeta values, in 1978 Apery ([lO], [11]) established that C(3) is irrational, see van 
der Poorten [244] and §3.151 for more details. In 1979 F. Beukers [31] gave an elegant 
proof of the irrationality of C(3) suggested by the form of Apery's proof, showing that 
certain integer linear combinations of period integrals defined for integer r, s > by 

Ir,s '■= I / -; x^y^dxdy= III jz —dzdxdy 



JO 



'^-xy Jq Jo Jo 1 - {1 - xy)z 



were very small. He evaluated 



'0,0 



C r-i^dxdy = 2C(3), (3.14.1) 
Jo ^-xy 



and for r > 1, 



Ir^r = 2(C(3) - 4- 4 4)= 2(C(3) - Hn,,). 



13 2^ r^' 

As a necessary part of the analysis he also showed for r > s that each I^.s = Is,r is a 
rational number with denominator dividing where is the least common multiple 
[l,2,...,r]. It is now known that infinitely many odd zeta values are irrational, without 
determining which ones (Rivoal }255] . Ball and Rivoal [21]). In particular Zudilin |315] 
showed that at least one of C(5), C(7), C(9), C(ll) is irrational. See Fischler jll4] for a 
survey of recent developments on irrationality of zeta values. The state of the art for 
showing transcendence of periods is surveyed by Waldschmidt [306] . 
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There is currently no effective way known to decide wlietlier a given number is a period. 
One can certify a number is a period by directly exhibiting it as a Q- linear combination of 
values of known period integrals. However no properties of periods are currently known 
which would be useful in distinguishing them from non-periods. Kontsevich and Zagier 
|18H Problem 3] raise the problem of exhibiting a single specific number that is provably 
not a period. 

As mentioned in the introduction it is conjectured: Euler's constant is not a period. 
The conjectures that e and 7 are not periods appears in |18H Sec. 1.1]. We have already 
seen a number of different integral representations for Euler's constant, which however 
involve exponentials and/or logarithms, in addition to rational functions. These integrals 
do not resolve the question whether 7 is a period. 

Kontsevich and Zagier |18H Sec. 4.3] also suggest enlarging the class of periods to 
allow exponential functions in the integrals. They define an exponential period to be an 
absolutely convergent (multiple) integral of the product of an algebraic function with 
the exponential of an algebraic function, taken over a real semi-algebraic set, where all 
polynomials entering the definition of the algebraic functions and the semi-algebraic set 
have algebraic coefficient^. The set of all such values forms a ring SV, which we call 
the ring of exponential periods, and which contains the field Q of all algebraic integers. 
This ring is countable, and is known to include the constant e, all values r(2) at rational 
arguments, namely | G Q \ '^<o, and the constant 



It also includes various values of period determinants for confluent hyper geometric func- 
tions studied by Terasoma [2901 Theorem 2.3.3] and by Bloch and Esnault \35\ Proposi- 
tion 5.4]. Here we observe that: 

(1) Euler's constant 7 S £V. 

(2) The Euler-Gompertz constant 6 G £V. 

M. Kontsevich observes that an integral representation certifying that 7 E £V is obtain- 
able from (|2.2.8|) by substituting — log x = for < a; < 1 and log x = for 
X > 1, yielding 



An integral representation certifying that 6 G £V is given by (I2.5.10p . stating that 



It is not known whether (5 is a period. 

3.15. Diophantine approximations to Euler's constant. Is Euler's constant ratio- 
nal or irrational? This is unknown. 

One can empirically test this possibility by obtaining good rational approximations 
to Euler's constant. The early calculations of Euler's constant were based on Euler- 
Maclaurin summation. In 1872 Glaisher |125j reviewed earlier work on Euler's constant, 
which included his own 1871 calculation accurate to 100 places ( |124j ). which uncovered 
a mistake in an earlier calculation of Shanks. In 1878 the astronomer J. C. Adams [2] 
determined 7 to 263 places, again using Euler-Maclaurin summation, an enormous labor 
by hand. In 1952 J. W. Wrench, Jr [312J . using a desk calculator, obtained 7 to 328 

■^^The semialgebraic set is cut out by a system of (real) polynomial inequalities fj{xi, Xn) > 0, and 
it is supposed that the coefficients of the fj are all real algebraic numbers. 
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places. In 1962 Knuth |174j automated this approach and obtained Euler's constant to 
1272 places with an electronic computer. This record did not last for long, as in 1963 
D. W. Sweeney |287j introduced another method to compute Euler's constant, based on 
the integral representation 

^ = ^ klk "I —du-logn, 
k=i ''^ 

By making a suitable choice of n in this formula, Sweeney obtained 3566 digits of 7. The 
following year Beyer and Waterman [37] used a variant of this method to compute 7 to 
7114 places, however only the initial 4879 digits of these turned out to be correct, as 
determined by later computations. In 1977 Brent |:4i] used a modification of this method 
to compute 7 and to 20700 places. In 1980 Brent and McMillan [33] formulated new 
algorithms based on the identity 

_ C/(n) _ KQ{2n) 
^ ~ ~ /o(2n) ' 

whose right hand side contains for = the modified Bessel functions 



^^(^) - g fc!r(' + fc + i) ■= -tJ^^' 

and 



00 



fc=0 

00 k 2 

Vin) := Io(2n) = E (^) • 



The second term < KQ{2n) / lQ{2n) < vre" can be made small by choosing n large, 
hence 7 is approximated by y^j, choosing n suitably (their Algorithm Bl). With this 
algorithm they computed both 7 and to 30000 places. They also determined the initial 
part of the ordinary continued fraction expansions of both 7 and , and established the 
following result. 

Theorem 3.15.1. (Brent and McMillan 1980) 

(1) If^ = ^ is rational, with mo,no G Z>o then its denominator uq > 10^^*^*^". 

(2) If = ^ is rational, with mi,ni G Z>o then its denominator ni > 10^^'^*''^. 

Proof. We use the argument of Brent (41j . For any real number 6 Theorem 17 of Khinchin 
|170j states: if ^ is an ordinary continued fraction convergent for 6 then \qn9 — Pn\ < 
\q9—p\ for all integers p and q with < |g| < It therefore suffices to find a convergent 
^ with g„ > 10^^'^'^'' such that ^ ^ 9, and the latter inequality holds if qn+i exists. □ 

The continued fraction expansion of 7 begins 

7= [0; 1,1, 2, 1,2, 1,4, 3, 13, 5, 1,1, 8, 1,2, 4, 1,1, 40, 1,11, 3, 7, 1,7, 1,1, 5, 1,49, 4,...], 

see \269\ Sequence A002852]. Brent and McMillan also computed various statistics for 
the initial part of this continued fraction, concerning the number of partial quotients of 
different sizes of the initial continued fractions of 7 and of e"^ , in order to test whether 7 
behaves like a "random" real number. They found good agreement with the predicted 
limiting distribution of partial quotients for a random real number 9 = [ao, oi, 02, • • • ]. 



66 



JEFFREY C. LAG ARIAS 



It is a theorem of Kuz'min that the distribution of the j-th partial quotient Prob[aj = A;] 
rapidly approaches the Gauss-Kuz'min distribution 

p(A;) :=log2(l + 7)-log2(l+ ^ ' 



k' k + l" 

as J — )• oo, see Khinchin |1701 III. 15]. 

Diophantine approximations to 7 can be extracted from various series expansions for 
7 involving rational numbers. In 1910 Vacca |304j found the expansion 



Llog2 



k 

k=l 

1 1\ /I 1 1 1\ /I 1 1 , 

+ 2- -- + --- + 3--- + — +••• 



7 



2 3/ V4 5 6 7/ V8 9 15 
Truncating this expansion at A; = 2m — 1 gives approximations to 7 satisfying 

n / Pm ^ 4(logm + 1) 
U < 7 < . 

Qm m 

In 2010 Sondow |274) gave refinements of this expansion yielding approximations with 
a convergence rate 0( ^°^J^ ) for an arbitrary but fixed r > 2. However approximations 
of these types fall far short of establishing irrationality of 7, since the denominators qm 
grow exponentially with m. 

Much recent work on Diophantine approximation of explicit constants 9 has focused 
on finding a series of rational approximations ^ with the properties: 

(1) Each sequence Un,Vn separately satisfies the same linear recurrence with coeffi- 
cients that are polynomials in the ring Z[n], where n is the recurrence parameter. 

(2) The initial conditions for the recurrences Un,Vn are rational, so that both se- 
quences consist of rational numbers. 

(3) One has ^ — )- as n — )• 00. 

We give examples of such recurrences below, e.g. (I3.15.ip . For such sequences the 
power series U{z) = Yl^=i'^nz'^ resp. V{z) = X^^i ^'^^^ each (formally) satisfy a 
homogeneous linear differential equation DF{z) = in the z-variable, whose coefficients 
are polynomials in Q[z], i.e. the operator has form 

j=0 

We will refer to such a sequence of approximations satisfying (1), (2) as being of H-type, 
and will call any real number obtainable as a limit of such approximations an elementary 
H -period. More generally, the ring T-LV of H -periods will be the ring generated over Q by 
elementary i?-period^. The ring T-LV is clearly countable, and its possible relations to 
rings of periods V or exponential periods £V discussed in ^3.14l remain to be determined. 

It turns out to be useful to single out a subclass of such numbers, which are the set of 
such approximations for which the linear differential operators D have a special property: 

G-operator Property. At some algebraic point zq G Q the linear differential operator 
with coefficients in C[2;] has a full rank set of a G -function solutions. 

The notion of G-function is presented in detail in Bombieri [38], Andre [1], [8], and 
Dwork, Gerotto and Sullivan [77]. We give a definition of a subclass of G- functions 
obtained by specializing to the rational number case. (The general case allows algebraic 
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number coefficients all drawn from a fixed number field K.) It is a function given by a 
power series F{z) = X^^g^""^"' with rational a„ such that 

(1) There is a constant C > such that \an\ < C", so that the power series F{z) 
has a nonzero radius of convergence. 

(2) There is a constant C" > such that the least common multiple of the denomi- 
nators of ai, ..,a„ is at most (C)". 

(3) The an = ^ are solutions of a linear recurrence having coefficients that are 
polynomials in the variable n with integer coefficients. 

The notion of G-operator was formulated in Andre [H IV. 5], [6l Sect. 3]. Andre's 
definition of G-operator is given in another form, in terms of it satisfying the Galochkin 
condition (from Galochkin [120]), but this class of operators coincides with the definition 
above. The minimal differential operator satisfied by a G-function is known to be a 
G-operator in Andr'e's sense [6l Theorem 3.2], by a result of D. and G. Chudnovsky. 
Conversely, any G-operator in Andre's sense has a full rank set of G-function solutions 
at any algebraic point not a singular point of the operator ([Q Theorem 3.5]).. 

Andre ^ p. 719] has shown that G-operators, viewed in the complex domain, are of 
a very restricted type. They must be Fuchsian on the whole Riemann sphere P^(C), i.e. 
they have only regular singular points, including the point oo, see also [8', Theorem 3.4.1]. 
A conjecture formulated in Andre [61, p. 718] is that all G-operators should come 
from (arithmetical algebraic) "geometry" , i.e. they each should be a product of factors 
of Picard-Fuchs operators, controlling the variation of cohomology in a parametrized 
family of algebraic varieties defined over Q. 

We will say that a real number 6 that is a limit of a convergent sequence of rational 
approximations ^ with Un,Vn coefficients in the power series expansions of two (rational) 
G-functions is an elementary G-period. As above we obtain a ring QV of G-periods, 
defined as the ring generated over Q by elementary G-periods. This ring has recently 
been shown in Fischler and Rivoal [1151 Theorem 3], to coincide with the field Frac(G), 
in which Frac(G) is the fraction field of a certain ring G constructed from G-function 
values. The ring G is defined to be the set of complex numbers /(a) where a E Q 
and /(z) is any branch of a (multi-valued) analytic continuation of a G-function with 
coefficients defined over some number field which is nonsingular at z = a. Fischler and 
Rivoal [1151 Theorem 1] show that the ring G coincides with the set of all complex 
numbers whose real and imaginary parts can each be written as /(I), with f{z) being 
some G-function with rational coefficients (as in the definition above) and whose power 
series expansion has radius of convergence exceeding 1. Their result [1151 Theorem 3] 
also establishes that the set of elementary G-periods is itself a ring containing all real 
algebraic numbers, which coincides with Frac(G) H M. Again, possible relations of the 
field GV with the Kontsevich-Zagier rings of periods V and extended periods £V remain 
to be determined 

As already mentioned in §3.141 a famous result of Apery (|10j. [11] ) in 1979 established 
that the period ("(3) is irrational. Apery's result was obtained by finding rational ap- 
proximations ^ to C(3) which were generated by two different solutions {pn, Qn ■ n > 1}, 
to the second order linear recurrence with polynomial coefficients: 

Un = (34n^ - Sln^ + 27n - 5)n„_i - (n - lfun-2 (3.15.1) 



■^^The paper [1151 Sect. 2.2] discusses the possibility that the equahty G = P = ^[^] might hold, 
and sketches an argument due to a reviewer that supports the conjecture that G C "p. 
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with initial conditions Po = 0,pi =6 and qo = l,qi = 5, respectively. Here are now 
called Apery numbers, and are given by 



while pn are given by 



see Fischler [114, Sect. 1.2]). These approximations are closely spaced enough and good 
enough to imply that, for any e > 0, and all q > q{e) one has 

|C(3)-^| > g-^+^for ah p G Z, 

with exponent 6 = 13.41782, which implies that C(3) is irrational. The method permits 
deriving families of approximations given by related recurrences for CI™-) for all m > 2, 
but only for m = 2 and 3 are they are good enough to prove irrationality. For each 
m > 2 the associated power series 



n>0 



n=0 



will be solutions to a linear differential equation with polynomial coefficients, and are of 
-ff-type. For m = 2, 3 these approximations are known to be of G-type. A conceptually 
new proof of irrationality of C(3) using the same approximations, but connecting them 
with modular forms on the congruence subgroup ri(6) of SL(2, Z) was presented in 1987 
by Beukers [33j. 

Apery's discovery was based on series acceleration transformations of the series C(3) = 
Z^^i initial example of Apery of such accelerations is the identity 

see [2441 Sect. 3]. The partial sums of this series converge at a rate 0(4""™) to fC(3), 
but are by themselves insufficient to prove irrationality. Apery's proof used further series 
acceleration transformations, described in [2441 Sect. 4]. It is interesting to note that the 
identity (I3.15.4p is a special case of identities in an 1890 memoir of A. A. Markoff [207] on 
series acceleration, as is explained in Kondratieva and Sadov [ 178] . In 1976 Gosper ( |133j . 
|134j ) noted related identities, also found by series acceleration techniques, including 

_ 30n - 11 

^"^il6(2n-l)n3(V)^- 

A new approach to Apery's irrationality proof of ^(3) was found by Nesterenko [226\ in 
1996, which used a relation of the Apery approximations Pn,qn to Fade approximations 
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to poly logarithms. Nesterenko also found a new continued fraction for 2("(3), 

2C(3)=2+ ^— . (3.15.5) 

2 + 



4 + 



1 



3 + 



2 + 



4 + 



4 

6 + 

9 

5 + 



4 + ... 

ai 

Writing this as 2(^(3) = 6o + t ) the numerators a„ grow quadratically in n while 

&! + ••• 

the denominators bn, grow linearly with n, and are given for k > (excluding ai) by 

a4k+i = k{k + 1), 64fc+i = 2A; + 2, 04^+2 = (^ + 1)(^ + 2), 64^+2 = 2/i; + 4, 
a4k+3 = {k + lf, 64A:+3 = 2A; + 3, 04^+4 = (A: + 2) 2, 64^+4 = 2^ + 2. 

This result was obtained starting from Pade type simultaneous polynomial approxima- 
tions to the logarithm, dilogarithm and trilogarithm found by Gutnik |149j . see also the 
survey of Beukers [32]. When specialized at the point x = 1 these Pade approxima- 
tions yield the sequences qn,—2pn- The continued fraction (I3.15.5p was obtained from 
a Mellin-Barnes type integral for a special case of the Meijer G- function [211], a gen- 
eralized hypergeometric function (see [Ml Chap. V]). A survey of methods using Pade 
approximations to polylogarithms is given in Nesterenko [227] . This approach has led to 
further results, including a new proof of irrationality of C{3) of Nesterenko [228]. 

A natural question concerns whether 7 has rational approximations of H-type. A 
family of such rational approximations to 7 were found in 2007 through the combined 
efforts of A. I. Aptekarev, A. I. Bogolyubskii, D. V. Khristoforov, V. G. Lysov and D. 
N. Tulyakov in the seven papers in the volume |12j . Their result is stated in Aptekarev 
and Tulyakov [15] in the following slightly improved form. 

Theorem 3.15.2. (Aptekarev, Bogolyubskii, Khristoforov, Lysov, Tulyakov 2007) Eu- 
ler's constant is approximated by a ratio of two rational solutions of the third-order 
recurrence relation with polynomial coefficients 

(16n-15)gn+i = (128n2 + 40n2-82n-45)^„-n2(256n^- 240n2 + 64n-7)gn-i 
+ (16n+l)n2(n-l)2g„_2. 

The two solutions {pn : n > 0} and {qn : n > 0} are determined by the initial conditions 

po := 0, pi := 2, p2 := 31, 

and 

qo := 1, qi := 3, q2 := 50. 
They have the following properties: 

(1) (Integrality) For n > 0, 

Pn e Q'n e 2- 

(2) (Denominator growth rate) For n> 1, 



qn = (2n)!^^ + 0{^] 

\/n V\/vr(4e)3/« ^ y/n' 
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(3) (Euler's constant approximation) Fern > 1, 

^ _ ^ = _2^e-2v^ f 1 + 0(^)y 

The original result showed in place of (1) the weaker result that qn & and DnPn £ ^ 
with Dn = l.c.m.[l, 2, n]. The integrality result for p„ follows from a later result 
of Tulyakov [300j . which finds a more complicated system of recurrences that these 
sequences satisfy, in terms of which the integrality of pn is manifest. These approximants 
clearly are of H-type. 

The results in Theorem 13 . 1 5 . 2 1 are based on formulas derived from a family of multiple 
orthogonal polynomials in the sense of Aptekarev, Branquinho and van Assche 
These lead to exact integral formulas 







Qn{x)e ^dx, 



and 



poo 

Pn-iqn = / Qn{x)e~''\ogxdx, (3.15.6) 
Jo 

in which Qn{x) are the family of polynomials (of degree 2n) given by 

Note for Qq{x) = 1 that (j3.15.6p becomes the integral formula (j2.2.8p for r'(l). Also in 
2009 Rivoal |256| found an alternative construction of these approximants that makes 
use of Fade acceleration methods for series similar to Euler's divergent series (I2.5.ip . His 
method applies more generally to numbers 7 + log x, where x > is rational. In 2010 Kh. 
and T. Hessami Pilehrood |159[ Corollary 4] found closed forms for these approximants, 

k=0 ^ ^ k=0 ^ ^ 

where the Hk are harmonic numbers. These rational approximations ^ thus are of a 
kind similar in appearance to Apery's approximation sequence ()3.15.2p . (j3.15.3p to C(3), 
and also satisfy a third order recurrence. 

Recently Kh. and T. Hessami Pilehrood |1601 Theorem 1, Corollary 1] constructed an 
elegant sequence of rational approximations converging to Euler's constant, also analo- 
gous to the Apery approximations to C(3), with slightly better convergence properties 
than those given in Theorem 13.15.21 and having a surprising connection with certain 
approximations to the Euler-Gompertz constant. They set 

«- = E(j^'' Pn = ^il] kl{2H^.k-Hk). 
k=0 ^ ^ fc=0 ^ ^ 

The sequence Qn satisfies the homogeneous second order recurrence 

qn+2 = 2{n + 2)qn+i - (n + ifqn- (3.15.7) 

and initial conditions = 1, (/i = 2. The sequence Pn satisfies the inhomogeneous second 
order recurrence 

Pn+2 = 2(n + 2)pn+l -{n + l) Pn -— , 

n + 1 

with initial conditions pq = 0,pi = 1, and it satisfies a homogeneous third order recur- 
rence. 
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Theorem 3.15.3. (Kh. and T. Hessami Pilehrood 2013) Let {qn)n>o, {Pn)n>o be defined 
as above. Then qn G 7^ and DnPn G where Dn = /.c.m.[l, 2, n], and for all n > 1, 

7-- = -e-^^(27r + 0(^)). 

Here the growth rate of the sequence qn is, for all n > 1, 

qn = nl^^l-^ + 0{^) 

These approximations converge to 7 from one side. The sequence needed to clear 
the denominator of p„ is well known to have growth rate 

r) _pn(l+o(l)) 

an asymptotic result that is equivalent in strength to the prime number theorem. 

The Euler-Gompertz constant b also has a convergent series of rational approximations 
with the same denominator sequence qn as the Euler constant approximations above. 
The new numerators Sn satisfy the same recurrence (|3.15.7p as the denominators qn, i.e. 
Sn+2 = 2(n + 2)sn+i — (n + l)^Sn, but with initial conditions sq = 0, si = 1. The Sn are 
integers, and are also given by the expression 

" (n — 1)\ 



Sn = Y,a{k-l){k + l)(''\ 
k=l ^ ^ 



(^-1)!' 

in which the function values a{m) are 

m 
k=0 

see [2691 Sequence A002793]. Here the a(m) are the partial sums of Euler's divergent 
series discussed in §2.51 These approximations |^ are exactly the partial quotients ^ of 
the Laguerre continued fraction (12.5. lOp for the Euler-Gompertz constant. 

One has the following rate of convergence estimate for this sequence of rational ap- 
proximations ( |160l Theorem 2]). 

Theorem 3.15.4. (Kh. and T. Hessami Pilehrood 2013) Let {qn)n>o, (sra)n>o be defined 
as above. Then qn,Sn G and for all n > 1, the approximations Sn/qn to the Euler- 
Gompertz constant satisfy 



The proof of Theorem 13. 15.41 relates these approximations to values of a Mellin-Barnes 
type integral 

^ r+^-r(s-n)2 

where the vertical line of integration has c > n. The integral is expressed in terms of 
Whittaker function W^^y{z) (a confluent hypergeometric function). The proof displays 
some identities involving the Euler-Gompertz constant and Whittaker function values, 
namely 

^ _ ^-iAo(l) 

using 

W^-l/2,0(l) = ^, W^l/2,0(1) = 
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n 


Qn 











1 





1 





1 


2 


1 


1 


1 


2 


7 


8 


2 


4 


3 


34 


118 


6 


20 


4 


209 


1450 


12 


124 


5 


1546 


53584 


60 


920 


6 


13327 


461718 


60 


7940 


7 


130922 


31744896 


420 


78040 


8 


1441729 


699097494 


840 


859580 


9 


17572114 


25561222652 


2520 


10477880 


10 


234662231 


341343759982 


2520 


139931620 



Table 2. Approximants Pn/ln and Sn/qn to 7 and 5, with Dn = [l,2,...,n]. 



In Tabled we summarize data on the approximations {QmPm 

) to 7 and 6, recalling 
that pn are rational numbers, while Qn, DnPn, Sn are integers, with Dn = [1,2, 

The approximations given in Theorems 13. 15. 3] and 13. 15. 41 certify that both 7 and 5 are 
elementary i7-periods, so that both 7, (5 E HV. These approximations are of insufficient 
quality to imply the irrationality of either Euler's constant or of the Euler-Gompertz 
constant. It remains a challenge to find a sequence of Diophantine approximations to 7 
(resp. 5) that approach them sufficiently fast to certify irrationality. 

One may ask whether 7 has Diophantine approximations of G-type. This is unknown. 
The current expectation is that this is not possible, i.e. that 7 QV. Based on the 
observation of Rivoal and Fischler that GV = Frac(G), and the the belief ( |1151 Sect. 
2.2]) that the ring G may coincide with the Kontsevich-Zagier ring V of extended periods, 
this expectation amounts to a stronger form of Conjecture 11.0.2] 

We conclude this section with recent results concerning G-type and H-type approx- 
imations to other constants considered by Euler. In 2010 Rivoal |258] showed that the 
periods r(^)'^, where ^ > , can be approximated by sequences ^ of rational num- 
bers of the G-type, so are elementary G-periods. In another paper Rivoal [257j finds 
sequences of rational approximations showing that the individual numbers r(^), which 
are not known to be periods for < A; < n, are elementary iJ-periods. 

Finally we note that various authors have suggested other approaches to establishing 
the irrationality of Euler's constant, e.g. Sondow |271j . |273j . and Sondow and Zudilin 
[277]. 

3.16. Transcendence results related to Euler's constant. The first transcendence 
results for numbers involving Euler's constant came from the breakthrough of Shidlovskii 
|264j . |265j in the 1950's and 1960's on transcendental values of functions, which is 
detailed in his book [266]. The class of E -functions was introduced in 1929 by Siegel 
\2<a7\ p. 223], consisting of those analytic functions F{z) = '^Z'^^QCn^ whose power 
series expansions have the following properties: 

(1) The coefficients c„ belong to a fixed algebraic number field K (a finite extension 
of Q). For each e > 0, the maximum of the absolute values of the algebraic 
conjugates of Cn is bounded by 0{n^'^) as n — )• 00. 

(2) For each e > such that there is a sequence of integers qo-,qi-,q2--- such that qnCn 
is an algebraic integer, with qn = 0{n'^^) as n — )• 00. 
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(3) The function F{z) satisfies a (nontrivial) linear differential equation 

j=Q 

whose coefficients are polynomials Rj{z) G where Ki is some algebraic 

number fieldF^ 

Each E'- function is an entire function of the complex variable z. The set of all i?- functions 
forms a ring E (under (+, x)) which is also closed under differentiation, and this ring is 
an algebra over Q. 

We also define a subclass of £'-functions called E* -functions in parallel with the class 
of G-functions, following Shidlovskii [266|, Chap. 13]. An analytic function F{z) = 
Yl'^=i '^«7ir ™ E* -function if it has the following properties: 

(1) F(z) is an E'-function. 

(2) There is some constant C > such that the maximum of the absolute values of 
the algebraic conjugates of c„ is bounded by 0(C") as n — >• oo. 

(3) There is some constant C > such that there is a sequence of integers qo,qi,q2--- 
such that each QnCn is an algebraic integer, and with g„ = 0{{C')'^) as n — )• oo. 

This definition narrows the growth rate conditions on the coefficients compared to Siegel's 
conditions. The set of all ii^* -functions forms a ring E* closed under differentiation which 
is also an algebra over Q. Certainly E* C E, and it is conjectured that the two rings are 
equal ([2661 P- 407]). Results of Andre ([6], [7j) apply to the class E*. 

Siegel \267\ Teil I] originally applied his method to prove transcendence of algebraic 
values of the Bessel function Jo{a) for algebraic a where Jo(q;) 7^ 0. In formalizing his 
method in 1949, Siegel [268J first converted an n-th order differential operator to a linear 
system of first order linear differential equations 

where the functions Rij{x) are rational functions. In order for his proofs to apply, 
he required this system to satisfy an extra condition which he called normal. Siegel 
was able to verify normality in a few cases, including the exponential function and 
certain Bessel functions. A great breakthrough of Shidlovskii [264] in 1959 permitted 
the method to prove transcendence of solution values for all cases where the obvious 
necessary conditions are satisfied. The Siegel-Shidlovskii theorem states that if the Ei{x) 
are algebraically independent functions over the rational function field K(x), then at 
any non-zero algebraic number not a pole of any function Rij{x), the values Ej{a) 
are algebraically independent. For treatments of these results see Shidlovskii |266j and 
Feldman and Nesterenko |llll Chap. 5]. 

More recent progress includes a study in 1988 of the Siegel normality condition by 
Beukers, Brownawell and Heckman [34j which obtains effective measures of algebraic 
independence in some cases. In 2000 Andre [U Theorem 4.3] (see also [8]) established that 
for each £'*-function the minimal order linear differential equation with C{z) coefficients 
that it satisfies has singular points only at z = and z = 00; that is, it has a basis of n 
independent holomorphic solutions at all other points. Furthermore 2; = is necessarily a 
regular singular point and 00 is an irregular singular point of a restricted type. This result 
goes a long way towards explaining why all the known i?*-functions are of hyper geometric 



'This condition on Ki could be weakened to require only that the coefficients are in C[z]. 
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type. Andre ^ applied this characterization to give a new proof of the Siegel-Shidlovskii 
theorem. 

Returning to Euler's constant, in 1968 Mahler |204j applied Shidlovskii's methods to 
certain Bessel functions. For example 

oo 

^o(^) 



(n!)2 4^ 

n=0 ^ ' 



is a Bessel function of the first kind, and 

n=l ^ ' 

is a Bessel function of the second kind, see [H (9.1.13)]. Here J^i^z) is an £'*-function, 
but io(-2)j which contains Euler's constant in its expansion (??), is not, because it has a 
logarithmic singularity at x = 0. However the modified function 

%{z) := ^Yo{z) - (log(|) +7) Jo(^) = E(-1)"''7;^(t)" 

n=l ^ 

is an £'*-function, and Mahler uses this fact. He noted the following special case of his 
general results. 

Theorem 3.16.1. (Mahler 1968) The number 

vr lo(2) 

2 Jo (2) ^ 

is transcendental. 

We have Jo(2) = E^=o 7^ and ^Yo{2) = jJo{2) + J]- ^(-l)-i^, from which 
it follows that Mahler's transcendental number is 

2 Jo(2) y-00 M);^- 

More recent transcendence results apply to Euler's constant along with other con- 
stants. In 2012 Rivoal [259] proved a result which, as a special case, implies the tran- 
scendence of at least one of the Euler constant 7 and the Euler-Gompertz constant S. It 
uses the Shidlovskii method and improves on an observation of Aptekarev [13] that at 
least one of 7 and S must be irrational. This transcendence of at least one of 7 and 6 was 
also established about the same time by Kh. and T. Hessami Pilehrood [16_0| Corollary 
3]. 



5 := dw 



Theorem 3.16.2. (Rivoal 2012) Euler's constant 7 and the Euler-Gompertz constant 

dv 

l + w Jq 1 - log f 

together have the following properties: 

(1) (Simultaneous Diophantine Approximation) One cannot approximate the pair (7, 5) 
very well with rationals (-, -) having the same denominator. For each e > there is a 
constant C{e) > such that for all integers p,q,r with q ^ 0, there holds 

q q H-^+'^ 

where H = max(|p|, \q\, \r\). 
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(2) (Transcendence) The transcendence degree of the field '[^{e^^^ 5) generated by e,j 
and 5 over the rational numbers Q is at least two. 

Proof. This follows from a much more general result of Rivoal's [259, Theorem 1], which 
concerns values of the function 



POO 

g^{z):=z-^ {t + zT-^e-'dt. 
Jo 



Here we specialize to the case a = 0, and on taking z = 1 we have 

r"Oo t 



poo -t 



n ■ nV 



using Hardy's integral (I2.5.10p . The function Qq{z) satisfies a linear differential equation, 
but is not an E'-function in the sense of transcendence theory. Rivoal makes use of the 
identity 

-l + \ogz = -e-'gQ{z)-8{-z), (3.16.1) 
valid for z G C \ M<o, where the entire function 

n=l 

is an i?*-function. 

Result (i) follows from the assertion of Theorem 1 (ii) of Rivoal j259J , specialized to pa- 
rameter values a = and z = 1. This result uses explicit Hermite-Pade approximants to 
the S*-functions l,e'',£a{z), where £a{z) = Y.m=o m!(m+a+i) ' ^'^^^ " -1,-2,-3,.... 
Here £{z) corresponds to a = — 1 with the divergent constant term dropped from the 
power series expansion. 

Result (ii) follows from the assertion of Theorem 2 (ii) of Rivoal |259j . which asserts 
for (complex) algebraic numbers z (— oo,0] that the field generated over Q by the 
three numbers e^, 7 + logz, Qq{z) has transcendence degree at least two. It makes use 
of ()3.16.ip . which can be rewritten as the integral identity 

--f = \ogz + zl e'^"" log t dt + e'"" dt 

Jo Jo t + z 

It uses the specialization z = 1, which essentially gives Hardy's identity (j2.5.1ip divided 
by e, which we may rewrite as 

-f(-l)=7 + -. 

e 

The result then follows using Shidlovskii's Second Fundamental Theorem [266^ p. 123]. 

□ 

Rivoal |259j remarks that the transcendence result Theorem 13.16.21 (2) is implicit in 
the approach that Mahler [204] formulated in 1968. At the end of his paper Mahler |204^ 
p. 173] states without details that his method extends to show that for rational vq > 
and nonzero algebraic numbers a, the integrals 

Ik{u,a) := [ x^«-^(logx)'=e-"^(ix. A: = 0,1,2,... 



are algebraically independent transcendental numbers. Choosing fo = l,a = 1 and 
= 0, 1, we have 



/o(l,l) := [\--dx = l--, 
Jo e 
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and, using (j2.5.9p . 

/i(l,l)= / {logx)e~''dx = / (logx)e~^<ix / log(x + l)e~''(ix = -7 . 

Jo Jo ^ Jo 6 

Mahler's statement asserts that these are algebraically independent transcendental num- 
bers. This assertion implies that Q(7, 5, e) has transcendence degree at least two over 
Q. 

We next present recent results which concern the transcendence of collections of gen- 
eralized Euler constants discussed in ^3.81 These results are applications of Baker's 
results on linear forms in logarithms, and allow the possibility of one exceptional alge- 
braic value. Recall that the Euler- Lehmer constants, for < /i < A; studied by Lehmer 
|196j . are defined by 



7(/i,/c) := lim ( V 



1 logx 
n k 

0<.n<x 
n=h ( mod k) 



In 2010 M. R. Murty and N. Saradha [2231 Theorem 1] obtained the following result. 

Theorem 3.16.3. (Murty and Saradha 2010) In the infinite list of Euler- Lehmer con- 
stants 

{7(/i, k) ■.l<h<k, for all k > 2}, 

at most one value is an algebraic number. If Euler 's constant 7 is an algebraic number, 
then only the number 7(2, 4) = j'j in the above list is algebraic. 

Their proof uses the fact that a large set of linear combinations of j{h, k) are equal to 
logarithms of integers. Then Baker's bounds for linear forms in logarithms of algebraic 
numbers are applied to derive a contradiction from the assumption that the set contains 
at least two algebraic numbers. An immediate consequence of Theorem 13. 16^ is the fact 
that over all rational numbers < x < 1 at least one of r(a;),r'(x) is transcendental, 
with at most one possible exceptional x. ( \223\ Corollary]). 

In a related direction Murty and Zaytseva |2241 Theorem 4] obtained a corresponding 
transcendence result for the generalized Euler constants 7(f^) associated to a finite set 
of primes that were introduced by Diamond and Ford |68j . Recall from ^3.81 that these 

constants are obtained by setting Cfiis) = Y\p\Q (^1 — P~^^C{s) aiid defining 7(il) as the 
constant term in its Laurent expansion around s = 1, writing 



Cc(^) = ^ + 7(f^) + f;7n(f^)0 
s — 1 ^-^ 

n=l 



In particular 7(0) = 7. 

Theorem 3.16.4. (Murty and Zaytseva 2013) In the infinite list of Diamond-Ford gen- 
eralized Euler constants 7(^2) where Q runs over all finite subsets of primes (including 
the empty set), all numbers are transcendental with at most one exception. 

The constants can be expressed as finite linear combinations of Euler-Lehmer 

numbers, cf. (I3.8.ip . so this result is closely related to the previous theorem. This result 
is also proved by contradiction, using Baker's results on linear forms in logarithms. If 7 
were algebraic then 7(0) would be the unique exceptional value. However one expects 
there to be no algebraic value in the list. 
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4. Concluding remarks 

Euler's constant appears sui generis. Despite three centuries of effort, it seems unre- 
lated to other constants in any simple way. 

A first consequence arising from this is: the appearance of Euler's constant in ostensi- 
bly unrelated fields of mathematics may signal some hidden relationship between these 
fields. A striking example is the analogy between results concerning the structure of 
factorizations of random integers and cycle structure of random permutations. These re- 
sults are given for restricted factorizations in ^3.5\ ^3.61 and for restricted cycle structure 
of permutations in ^3.101 ^3.111 respectively. Results in these two fields were initially 
proved separately, and were later unified in limiting probability models developed by Ar- 
ratia, Barbour and Tavare [T7]. These models account not only for the main parts of the 
probability distribution but also for tail distributions, where there are large deviations 
results in which Euler's constant appears. On a technical level an explanation for the 
parallel results in the two subjects was offered by Panario and Richmond \2S7\ at the 
level of asymptotic analysis of generating functions of a particular form. This research 
development is an example of one of Euler's research methods detailed in §2.61 that of 
looking for a general scheme uniting two special problems. 

In a second direction, the study of Euler's constant and multiple zeta values can now 
be viewed in a larger context of rings of periods connected with algebraic integrals. The 
sui generis nature of 7 is sharply formulated in the conjecture that it is not a Kontsevich- 
Zagier period. On the other hand, it is known to belong to the larger ring of exponential 
periods £V, as described in §3.141 This research development illustrates another of 
Euler's research methods, to keep searching for a fuller understanding of any problem 
already solved. 

In a third direction, one may ask for an explanation of the appearance of Euler's 
constant in various sieve method formulas. Here Friedlander and Iwaniec [118^ p. 239] 
comment that these appearances stem from two distinct sources. The first is through 
Mertens's product theorem 13.3.21 The second is via differential-difference equations of 
the types treated in ^3.51 and ^3.61 

In a fourth direction, there remain mysteries about the appearance of Euler's constant 
in number theory in the gamma and zeta function. 

Observation. Euler's constant 7 appears in the Laurent expansion of the Riemann 
zeta function at s = 1; here the zeta function is a function encoding the properties of 
the finite primes. Euler's constant also appears in the Taylor expansion of the gamma 
function at s = 1; here the gamma function encodes properties of the (so-called) "infinite 
prime" or "real prime" (Archimedean place). 

The analogy between the "real prime" and the finite primes traces back to Ostrowski's 
work in the period 1913-1917 which classified all absolute values on Q ("Ostrowski's 
Theorem"). Such analogies have been pursued for a century; a comprehensive treatment 
of such analogies is given in Haran [150j . 

A mild mystery is to give a conceptual explanation for this occurrence of Euler's 
constant in two apparently different contexts, one associated to the finite primes, and 
the other associated to the real prime (Archimedean place) . There is a known connection 
between these two contexts, given by the product formula for rational numbers over all 
places, finite and infinite. It says 

JJ|r|p = l, rGQ. 

p<oo 
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in which | • \p denotes the p-adic valuation, normahzed with \p\p = - and |x|oo being 
the real absolute value. Does this relation explain the common appearance? On closer 
examination there appears to be a mismatch between these two occurrences, visible in 
the functional equation of the Riemann zeta function. The functional equation can be 
written 

in which appears the completed zeta function 

as) :=vr-ir(i)C(.)=vr-fr(£)n(l-^)"'- 

The completed zeta function is given as an Euler product over all "primes", both the 
finite primes and the Archimedean Euler factor 7r~2r(|). In this formula evaluation at 
the point s = 1 for the Riemann zeta function corresponds to evaluation of the gamma 
function at the shifted point = f = ^- The mismatch is that for the gamma function, 
Euler's constant most naturally appears in connection with its derivative at the point 
s' = 1, rather than at the point s' = ^, cf (|3.1.ip . Is there a conceptual explanation of 
this mismatch? Theorem 13.1.11 associates Euler's constant with the digamma function 
at s = ^, but in that case an extra factor of log 2 appears. Perhaps this mismatch can 
be resolved by considering an asymmetric functional equation for the function r(s)(^(s). 
Here one has, valid for Re(,s) > 1, the integral formula 

Jo 

One of Riemann's proofs of the functional equation for (^{s) is based on this integral. 

Another mystery concerns whether there may exist an interesting arithmetical dynam- 
ical interpretation of Euler's constant, different from that given in ^3.121 The logarithmic 
derivative of the archimedean factor 7r~2r(|) at the value s = 2 (that is, s' = 1) produces 
the value —g/2 where 

g ■= 7 + log7r. 

The constant g appears elsewhere in the "explicit formulas" of prime number theory, as 
a contribution at the (special) point x = 1 in the Fourier test function space, cf. [39l p. 
280]. Perhaps one should search for a dynamical system whose entropy is the constant 
-g ■= - logTT « -1.721945. 

In a fifth direction, Euler's constant can be viewed as a value emerging from reg- 
ularizing the zeta function at s = 1, identified with the constant term in its Laurent 
series expansion at s = 1. There are a number of different proposals for "renormalizing" 
multiple zeta values, including the value of "^(1)", with the intent to preserve other 
algebraic structures. This is an active area of current research. Various alternatives are 
discussed in Ihara, Kaneko and Zagier [167], Guo and Zhang [147] , |148| . Guo, Paycha, 
Xie and Zhang [14'6j and Manchon and Paycha [206J . Some of these authors (for example 
|167| ) prefer to assign the "renormalized" value to "C(l)", to be obtained by further 
subtracting off Euler's constant. 

This direction generalizes to algebraic number fields. For background we refer to 
Neukirch |230t Chap. VII. §5]. The Dedekind zeta function C,k{s) of a number field K 
is given by 

where A runs over all the ideals in the ring of integers Ok of the number field, and Nj^^q 
is the norm function. This function analytically continues to a meromorphic function on 
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C which has a simple pole at s = 1, with Laurent expansion 

Ck{s) = + ao{K) + f2a,{K)is - ly . (4.1) 

i=i 

The term a-i{K) encodes arithmetic information about the field K, given in the analytic 
class number formula 

ReSs=i{CKis)) ■■= a-i(i^) = ^^-^^ , 

in which hx is the (narrow) class number of K, Rk is the regulator of K (the log- 
covolume of the unit group of K) and qk is the genus of K, given by 

9K--= log 

in which wk is the number of roots of unity in K, dx is the absolute discriminant of 
ri and r2 denote the number of real (resp. complex) Galois conjugate fields to K, inside 
Q. The constant term aQ{K) is then an analogue of Euler's constant for the number 
field since ao(Q) = 7- A refined notion studies the constant term of a zeta function 
attached to an individual ideal class C of a number field K. Such a partial zeta function 
has the form 

where the sum runs over all integral ideals A of the number field belonging to the ideal 
class C. The Dedekind zeta function is the sum of ideal class zeta functions over all 
hx ideal classes. These partial zeta functions C,{C,s) extend to meromorphic functions 
of s on C with only singularity a simple pole at s = 1 having a residue term that is 
independent of C, i.e. a-i{C,K) = -^a-i{K) (|230l Theorem VII (5.9)]). 

For imaginary quadratic fields the constant term ao(C,i^) in the Laurent expansion 
at s = 1 of an ideal class zeta function can be expressed in terms of a modular form 
using the first Kronecker limit formula. This limit formula concerns the real-analytic 
Eisenstein series, for t = x + iy with y > and s € C 



E{r, s)= Yl 



Imr + np* ' 

(m,n)^(0,0) ' ' 

which for fixed r analytically continues to a meromorphic function of s G C and has a 
simple pole at s = 1 with residue tt and Laurent expansion 

oo 

E{t, ,) = ^ + co(r) + ^ Cn{T){s - ir. 
s — 1 ^-^ 
i=i 

The first Kronecker limit formula expresses the constant term co(t) as 

Co(r) = 27r(7 - log 2) - log(Vy|r/(r)|2), 
in which r/(r) is the Dedekind eta function, a modular form given by 

oo 

n(T) = J](l - g'^), with q = e^^'\ 

n=l 

see j313l Sect. 2]. In 1923 Herglotz |157j found an analogous formula for the constant 
term in the Laurent series expansion at s = 1 of C,{C,s) for an ideal class C of a real 
quadratic field, which involved Dedekind sums and an auxiliary transcendental function. 
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In 1975 by Zagier [313] obtained another formula for this case which is expressed using 
continued fractions and an auxihary function. 

In 2006 Ihara |168] studied another quantity extracted from the Laurent expansion 
(j4.ip of the Dedekind zeta function of a number field at s = 1, given by 



which he termed the Euler-Kronecker constant of the number field K. Here "y{K) rep- 
resents a scaled form of the constant term, and it generalizes Euler's constant since 
7(Q) = 7- Motivated by the function field case, Ihara hopes that ^{K) encodes inter- 
esting information of a different type on the field K, especially as K varies. Under the 
assumption of the Generalized Riemann Hypothesis, Ihara obtained upper and lower 
bounds 



where Dk is the (absolute value of the) discriminant of K and ci,C2 > are absolute 
constants. Further unconditional upper bounds on ^{K) were obtained by Murty [225j 
and lower bounds were given by Badzyan |22j . 

Finally, ongoing attempts to place Euler's constant in a larger context have continued 
in many directions. These include studies of new special functions involving Euler's 
constant, such as the generalized Euler- constant function 7(2;) proposed by Sondow and 
Hadjicostas [ 275] . Two other striking formulas, found in 2005 by Sondow |272j . are the 
double integrals 



Guillera and Sondow |145j then found a common integral generalization with a complex 
parameter that unifies these formulas with double integrals of Beukers [31] for C,{2) and 
C(3), see (|3.14.ip . 

Many other generalizations and variants of Euler's constant have been proposed. 
These include a j?-adic analogue of Euler's constant introduced by J. Diamond [69j in 
1977 (see also Koblitz |176j and Cohen [55^, Sect. 11.5]) and g-analogues of Euler's con- 
stant formulated by Kurokawa and Wakayama [183J in 2004. A topic of perennial interest 
is the search for new expansions and representations for Euler's constant. For example 
1917 work of Ramanujan |249] was extended in 1994 by Brent [42j. Some new expan- 
sions involve accelerated convergence of known expansions, e.g. Eisner [82], K. and T. 
Hessami Pilehrood ^158j . and Prevost |245j. Recent developments of the Dickman func- 
tion led to the study of new constants involving Euler's constant and polylogarithms, cf. 
Broadhurst [47] and Soundararajan [279] . 

All these studies, investigations of irrationality, connections with the Riemann hy- 
pothesis, possible interpretation as a (non)-period, and formulating generalizations of 
Euler's constant, may together lead to a better understanding of the place of Euler's 
constant 7 in mathematics. 
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